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THE ROLE OF NUMERICAL ANALYSIS IN AN 
UNDERGRADUATE PROGRAM* 


GEORGE E. FORSYTHE, Stanford University 


1. Where do mathematics students go? Before considering the aim of under- 
graduate mathematics education we should recall approximately who the stu- 
dents are who take college mathematics courses, and where they are going. I 
assume that everywhere, as at Stanford, most mathematics students major in 
other fields. Usually students of science and engineering form most of our clien- 
tele. The increasing role of mathematics in economics, psychology, etc., is bring- 
ing a number of serious students of social science. And we have an increasing 
number of mathematics majors—students who are naturally dearest to our 
hearts. 

It should be obvious that the engineers and other nonmajors are primarily 
interested in the application of mathematics, and not in its structure. The best 
of the nonmajors will apply mathematical analysis in their professions. A large 
number of the less gifted, particularly among the engineers, will face a con- 
siderable number of computations in their future work, and this may represent 
the main future use of their mathematical background. 

And our majors, where are they bound? One wishes in vain that all the best 
of them would either teach high school or take graduate work. Of those who 
attain the Ph.D., about half go into industry and apply mathematics in various 
areas. Almost all of them will be connected to some degree with automatic 
computation. But the majority of our undergraduate mathematics majors are 
lured at once into the market place, where they are greatly in demand as serv- 
ants of the fast-multiplying family of fast-multiplying computers. 

Why are they in such demand? There seem to be over 3000 automatic digital 
computers now installed in the United States, with more on the way. As a 
rough estimate, each automatic computer needs to have 10 attendants who serve 
it as mathematicians—programmers, coders, analysts, supervisors, etc. The re- 
sulting requirement for 30,000 computer mathematicians should be compared 
with the combined membership of the American Mathematical Society, Mathe- 
matical Association of America, Society for Industrial and Applied Mathe- 
matics, Association for Computing Machinery, Institute of Mathematical Sta- 
tistics, and American Statistical Association—under 20,000 persons. While some 
makeshift arrangements are possible, the disparity in numbers is creating the 
unprecedented demand (and salary) for the new A.B. in mathematics. 

It must be noted as a digression that as a practitioner of mathematics, the 
new A.B. is not in a strong position, despite industry’s demand for him. He 
knows little of what he needs to do a good job—superficially because his teachers 
have not prepared him to be a practitioner of mathematics, but fundamentally 

* Invited address at the Thirty-ninth Summer Meeting of the Mathematical Association of 
America, August 27, 1958. The preparation of the manuscript was sponsored in part by the Office 
of Naval Research under project NR 044 211. 
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because of the brevity of his mathematical education. As a mathematically 
half-educated man he needs to learn both caution and confidence. But, if he is 
confident, wise in human relations, aggressive, and lucky, he will become a 
supervisor of coders—a position in which with his half education he sometimes 
does serious damage to his company and the reputation of scientific computing. 
As far as technical advancement is concerned, coding appears to be a terminal 
occupation, except to those with considerable graduate work in mathematics. 

For discussions of mathematicians in the market place, see Rees [12] and 
Fry [5]. 

In summary, a very large number of undergraduate students of mathematics 
are going to apply mathematics in their future work, and much of this applica- 
tion will involve automatic computation. How should this influence our pro- 
gram? 


2. Aims of undergraduate mathematics education. In the author's opinion, 
we should have the following goals for the mathematical education of our 
undergraduates, whether majors or not. 


a. The student should above all learn as much as possible of the structure» 
significant ideas, and points of view which constitute mathematics. This is the 
sine qua non. 

b. The student should learn to read independently the mathematical litera- 
ture at his level of comprehension, and to use mathematical language to write 
facts, proofs, and ideas precisely and unambiguously, in a style acceptable to 
not-too-strict editors of mathematical journals. With the growth of applied 
mathematics, a good many papers are being written by the A.B. in mathe- 
matics; it is unfair to the authors, to editors, and to readers that so many of 
their authors write as mathematical illiterates. 

c. The student should know the tools of the mathematical profession (books 
and machines), and how and where to find them and to keep up to date on them. 
Within his general limitations, he should be competent in the library. He 
should know roughly what bibliographies exist, what tables of useful functions 
there are, what machines can and cannot do, and so on. Since he is fully com- 
petent to learn it, he should be able to code one automatic computing machine 
quite well. Finally, he should have a general idea of what practitioners of pure 
and applied mathematics there are, and their abilities. For example, we must 
teach the future engineer to consult with statisticians, numerical analysts, and 
other experts at an appropriately early time, instead of trying to go it alone 
until an emergency arises. 

d. The student should have cultivated and practiced the solution of mathe- 
matical problems new to him. This is where research ability is built. It hardly 
matters whether the problems have been solved before by some one else. 

e. The student should have gone fairly deeply into some other field of 
knowledge—preferably some field where mathematics is used, like a physical 
science or the mathematical aspects of economics. He should have acquired as 
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much experience as possible in communicating with practitioners of this other 
field, so as to learn to translate mathematical terms into their language, and vice 
versa. Here he should have struggled with the mathematical idealization of 
problems from the real world. 

f. The student should have learned to enjoy mathematical study. 


It seems to me that these constitute valid goals also for mathematics in 
grade school, secondary school, and graduate school, if the standards of writing 
and item c are modified appropriately. In graduate school one might wish to add: 


g. The student should receive some supervised experience in serving as a 
mathematical consultant. 


3. The abstract versus the concrete. In this post-sputnik era one hears and 
reads a great deal of criticism of mathematical education at all levels, but par- 
ticularly in high schools. Many, perhaps most, mathematical critics deplore the 
lack of rigor and abstraction, and the lack of understanding of the mathematical 
structure of algebra, geometry, and other subjects. I believe that criticism to be 
valid. But I have a different complaint—that much teaching is not intuitive 
enough nor well enough illustrated. I wish mathematics classes would get down 
to earth more often and make actual application of newly learned theorems to 
concrete problems of a realistic nature. 

In these ideas I agree with Felix Klein [9], and quote from his book, pages 
15 and 16: “I have already emphasized the fact that, in the schools, applications 
accompany arithmetic from the beginning, that the pupil learns not only to 
understand the rules, but to do something with them. And it should always be 
so in the teaching of mathematics! Of course, the logical connections, one might 
say the rigid skeleton in the mathematical organism, must remain, in order to give 
it its peculiar trustworthiness. But the living thing in mathematics, its most im- 
portant stimulus, its effectiveness in all directions, depends entirely [“entirely” 
is too strong—G.E.F.] upon the applications, i.e., upon the mutual relations 
between those purely logical things and all other domains. To banish applica- 
tions from mathematics would be comparable to seeking the essence of the liv- 
ing animal in the skeleton alone, without considering muscles, nerves and 
tissues, instincts, in short, the very life of the animal.” Klein goes on to say 
that, while in research differentiation between pure and applied mathematicians 
may be essential, such differentiation is not reasonable in teaching. 

It is not contradictory to criticize in the same breath mathematics teaching 
for not being abstract enough and for not being concrete enough. I think both 
these criticisms are frequently valid. Taken together, they imply a hazy mathe- 
matics course which has neither the beautiful structure of a logical entity nor 
fascinating application to the physical world. To the extent that mathematics 
courses are hazy and shapeless, the author claims there are two directions in 
which reform is needed—to strengthen the skeleton, and to add a variety of 
applications. 

For example, a student with a poor knowledge of high school algebra is weak 
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on structure—for example, he is unaware that the whole subject flows from a 
few postulates and principles, or that proofs are just as relevant to algebra as to 
geometry. But he is also weak in the application of algebra. He is notably unable 
to convert word problems into algebraic form—always a hard part of algebra. 
He is also weak in arithmetic, and lacks good intuition about what happens when 
one gives variables numerical values. (Can — x really be positive?) 

The use of high school algebra is purely for an example. Quite parallel in- 
stances occur everywhere in school and college. 

It must be emphasized that the concrete application of an abstract theory is 
rarely in practice just a special case of the abstract theory! To illustrate this 
paradox, consider the multiplication of real numbers x, y. We know that xy =yx, 
and that x(yz) =(xy)z. However, in the concrete application to digital comput- 
ers, one is forced to curtail the decimal (or binary) expansions of numbers at 
some fixed precision. Thus, x, y, s are represented by digital numbers #, 9, 2. 
Moreover, exact multiplication of digital numbers does not lead to a digital 
number, and so one must be content with a pseudomultiplication X, where the 
pseudoproduct is digital. Difficulties arise from the fact that often #X9#9XzZ 
and £X(jXzZ)#(#X9) Xz. This is perhaps a typical example of essentially 
new problems created by concrete applications of a mathematical theory. The 
problems may in turn lead to a new mathematical theory to solve them. 

Since my personal interest is in the concrete side of mathematics, I shall 
concentrate on it in the following, leaving structure to others more interested 
and better qualified. For other opinions on the teaching of concrete mathematics, 


see Tukey [14] and Weyl [15]. Tukey refers to monastic and secular mathe- 
matics. 


4. What is numerical analysis? Numerical analysis (apparently a postwar 
term) is a branch of applied mathematics, namely the art and science of using 
digital computation to solve scientific problems. Traditionally numerical analy- 
sis included such matters as numerical quadrature and the solution of poly- 
nomial equations and ordinary differential equations. Since automatic digital 
computing machines are newly on the scene, the part of numerical analysis deal- 
ing with their use in scientific computation is the live area for research. And, 
in practice, the large machines are dominating scientific computation in this 
country. Hence the problems of using such machines constitute most of what I 
consider to be numerical analysis today. Since the machines can deal with vast 
quantities of data, the solution of large problems is frequently asked. For exam- 
ple, one thinks today of solving linear algebraic systems with 100 to 1000 un- 
knowns, where, with a desk calculator, one with 10 unknowns requires most of a 
day. But, even more than size, the automatic nature of these machines has 
changed the character of problems. We now consider that numerical analysis 
covers any type of problem which an automatic computer should be able to 
solve—any problem which can be reduced to an algorithm or to a flow chart. 
At this point we consider such things as formal integration of a rational function 
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as numerical analysis. In fact, most exercises other than “word problems” and 
proofs which occur in undergraduate courses involve some algorithm, and may 
be considered as part of numerical analysis! 

Numerical analysis has different branches. Among these are the design of 
algorithms, the analysis of errors, and techniques for coding (mainly the use of 
machines to prepare codes for machines). 


5. What numerical analysis should the undergraduate study? The amount 
of analysis and algebra at the command of the undergraduate limits the depth 
to which he can pursue numerical analysis, and in particular error analysis. This 
tends to force the undergraduate numerical analysis course into a collection of 
recipes (simple algorithms). But this is hardly an excuse for a course, because 
the student should be well enough educated to be able to read and apply 
algorithms for himself. To define an appropriate syllabus for undergraduate 
numerical analysis is therefore a bit delicate. I consider the following to be 
about right: 


a. He should learn to code one automatic digital computing machine. Actual 
experience operating a machine is certainly valuable and stimulating, but it is 
not essential. One can learn to code any machine without seeing it, and indeed 
all machines are coded before they are built. Except for the specialist, it will 
suffice to learn to code in an algebraic language like FORTRAN, which can be 
transformed by special translator codes into machine-language codes. A stand- 
ard international algebraic language for computers, called ALGOL, is now being 
developed. 

b. He should get a feeling for what automatic computers can do, and what 
they cannot do. He should learn that they are no substitute for creative thought, 
and yet that they can do a good deal of what passes for thought in this world. 
The student must learn the distinction. He should also learn that collaboration 
with an automatic computer compels precise formulation of the problem, and 
100 percent accuracy in preparation of the code. In this respect the automatic 
computer really forces that precision of thinking which is alleged to be a product 
of any study of mathematics. 

c. He must learn the tools of the computer—what tables, machines, books, 
codes, bibliographies, and people exist to help him solve future computing 
problems. 

d. He should learn a few typical algorithms, including, for example, New- 
ton’s method for finding zeros, Crout’s method for solving linear systems, and 
Euler’s method for ordinary differential equations. If not too many algorithms 
are undertaken, it should be possible to study each fairly deeply. 

e. Heshould have been exposed to some instances of treachery in computing, 
enough to teach him great caution. For example, he might find out that, whereas 
the polynomial 


F(x) = — — 2) -- (2 — 20) = — + + 201 
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has zeros 1, 2,---, 20, the slightly rounded polynomial F(x)—2-**x!* has 
among its zeros numbers near 20.846 and 13.99 +2.5i (Wilkinson [16]). This 
example should keep any one from using polynomial-solving routines blindly! 
We hope it would teach the student that you still have to think in the era of 
automatic computers. 

f. He should have some exposure to careful error analysis. For example, he 
might have studied Wilkinson’s round-off error analysis of Crout’s method 
[17]. A fruitful point of view here is to see how much one has to change the 
data to make the computed answer correct, rather than to ask how wrong the 
answer is for the given data. 

g. He should have a firsthand acquaintance with arithmetical operations on 
numbers, gained from desk computation and also, if possible, from work with 
an automatic computer. 

Professor Perlis has remarked that, whereas we think we know something 
when we learn it, and are convinced we know it when we can teach it, the fact 
is that we don’t really know it until we can code it for an automatic computer! 


6. How can we teach this numerical analysis? There seem to be two ways of 
teaching numerical analysis within the course structure of the American college. 
One is to teach special courses in it, and the other is to mix it into other relevant 
mathematics courses. Obviously this year’s choice between these methods de- 
pends on the present existence of teaching materials and trained staff. But one 
may also ask which would be the preferable solution, assuming that teaching 
materials and staff will eventually support the prevailing decision. 

Many professional numerical analysts, the author included, agree that 
numerical analysis should be mixed into most undergraduate mathematics 
courses. From our point of view, many undergraduate exercises in such subjects 
as calculus, linear algebra, and differential equations are in fact solvable with 
automatic computers, and hence belong to numerical analysis. We think their 
real importance is the illumination of mathematical theory, and hence that 
such exercises should stay where they are now: in the theoretical courses. How- 
ever, to pick one example from linear algebra, we feel that the conventional 
system of 3 linear equations in 3 unknowns with small integral coefficients does 
not furnish a sufficiently bright illumination of linear algebraic theory. And it 
certainly fails dismally as a fair illustration of linear algebra in industry today. 
We therefore feel that more realistic examples should be considered in the linear 
algebra course. One could easily consider 10 or 20 equations in 10 or 20 un- 
knowns, and expect students to produce a complete flow chart, or a code in 
some algebraic language like ALGOL. 

Turning to ordinary differential equations, we feel that no one understands 
the real nature of a differential equation dy/dx=f(x, y) until he has experi- 
mented with computational or graphical procedures for its solution. Such a 
simple process as Euler’s method, while not accurate, is both useful and worthy 
of study. What normally passes for numerical methods in differential equations 
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courses bears little resemblance to how such equations are solved in practice. 
Why not mix realistic practice right with the theory, where they can influence 
and strengthen each other? 

There is one exception, a sub’ect sufficiently demanding and sufficiently 
remote from other undergraduate mathematics to warrant a special course. 
This is coding. As we have indicated, we think every undergraduate mathe- 
matics student should know how to code some machine fairly well. (I would 
also include all undergraduate students, for I feel that the computer revolution 
will have such a great impact on all our lives that every college graduate should 
understand it intimately. Possibly it will eventually be taught in the ninth 
grade for the same reason.) Since coding presupposes no mathematics beyond 
arithmetic, it can be taught to freshmen. I recommend a two-hour-per-week 
semester course in coding, to be taken as early as possible. There are texts 
already available for such a course, including one by Andree [1]. The coding 
course should also include the material in c of Section 5—an exposure to bibli- 
ographies, tables, codes, and other library material. For, if numerical analysis is 
integrated with other mathematics, the student may nowhere else be taught by 
specialists in numerical analysis. 

Until a student actually has a digital computer to work with, the greatest 
benefit of a coding course is the point of view. The student will learn what con- 
stitutes creative thought, and what is merely routine. The flow chart, almost es- 
sential to coding, is one of the greater tools in scientific and technological think- 
ing. It probably originated in the block diagram of the electrical engineer. It 
would be a useful device in mathematical teaching and thinking. For example, 
for learning the basic concept of functions of one variable and their composition, 
the following block diagram seems far superior to the less graphic functional 
notation y=f(g(x)): 


With coding behind him, the student can analyze many algorithms from the 
point of view of preparing an actual or potential code to solve them. And, to the 
extent that he has. a machine available, he can solve realistic problems in any 
area. 

I should like to point out that numerical analysis was mixed into a theoreti- 
cal course in mathematics by Dickson in his book [3]. Such a course in the 
theory of equations was taught for years in American colleges, and is only now 
beginning to disappear in favor of more linear algebra. Dickson and other 
authors following his lead included Horner’s and Newton’s root-solving pro- 
cedures, evaluation of determinants, solution of linear algebraic systems, and so 
forth. With the evolution of computing equipment and numerical analysis, 
much of Dickson’s explicit and implicit advice on numerical analysis became 
obsolete and misleading, but for many of us it was the only numerical analysis 
we ever saw in college. 

Let us consider the common practice of teaching one or more separate courses 
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in numerical analysis for undergraduates. These courses have not required cod- 
ing as a prerequisite, although this is changing now. Judging from the textbooks, 
some of these courses tend to be filled more with algorithms than with analysis, 
and thus may be shallow for competent students. Milne’s book [11] is excellent, 
but most of it could be read by the student whenever he needed a certain 
technique. The basic mathematical concepts of the book, like orthogonal poly- 
nomial systems and the representation of linear functionals, should really be 
covered in other courses. 

The other extreme in a course is represented by Hildebrand’s beautiful book 
[7]. This includes real analysis of round-off errors, truncation errors, and so on, 
and is a great pleasure to a numerical analyst. But the typical good under- 
graduate is simply not enough of an analyst to do justice to the course, if I can 
judge from Stanford and U.C.L.A. While the instructor can teach the analysis 
as he goes along, it is questionable whether a special numerical analysis course 
is the correct place for this. 

In at least one university (U.C.L.A.) there are four undergraduate numerical 
analysis courses covering the fields of interpolation, linear algebra, ordinary 
differential equations, and approximation, respectively. These are excellent 
courses, and there seems to be no real reason for not merging most of them with 
basic undergraduate courses in corresponding subjects—say advanced calculus, 
linear algebra, ordinary differential equations, and Fourier series. Both sets of 
courses would profit from the merger. If coding were common knowledge of all 
mathematics students, there would seem to be no necessary obstacle to the 
merger. 

As nearly as one can tell from Gnedenko’s article [6], the Russian practice 
seems to be to include numerical analysis problems in a mathematical labora- 
tory. In his article Gnedenko lists several required problems in numerical analy- 
sis. These problems are large ones, and require plenty of both analysis and com- 
putation. How the laboratory is coordinated with the mathematical lectures, 
we have no idea. 

In summary, I favor a special coding course for all students, to be taken 
early. I should like to see material prepared so that numerical analysis could 
otherwise be sifted in with the basic mathematical theory, as illustrations of the 
theory and as motivation for more theory. In the next section are examples of 
where this mixing might occur. 


7. Examples of numerical analysis mixed into other courses. In this section 
I shall outline places where numerical analysis comes into standard courses in 
college mathematics, and in particular where problems should be treated with 


flow charts or codes. Naturally these vary from small problems to quite large 
ones. 


Analytical geometry. Almost any exercise, like finding the bisector of the 


angle between two given lines ax+by+c=0 and aix+biy+c:=0, where the 
constants are arbitrary. 
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Finding areas and volumes included by lines or planes. 

Testing for concurrence or collinearity. 

Determining the type of a conic corresponding to a second-degree equation 
in x and y. 

Calculus. Formal integration, e.g.; making a flow chart for integrating a 
general rational function of x of degree no higher than 6, considering all cases. 

Studying the asymptotic convergence of Newton’s process for getting real 
zeros of a function. What about multiple roots? Bounding the error in a pre- 
sumed root. 

Use of differentials to get the linearized expression for the error in evaluating 
a polynomial. Use of polynomials of degrees like 20 or 100 for examples. 

Study of sequences of difference quotients whose limit is dy/dx. Accelerating 
the convergence of the sequences. Central differences versus one-sided differ- 
ences. 

Numerical quadrature on an automatic computer. Study of sequences of 
quadratures by the trapezoidal rule, as h-0. 

The error in the trapezoidal rule, when f’’(x) is integrable, or not integrable. 

Upper and lower bounds for integrals. Study of Riemann sums as a sequence 
to be accelerated to its limit, the Riemann integral. 

Effect of rounding numbers on the sharp discontinuity in the formula for 
Jxndx when n= —1. See Forsythe [4]. 


Algebra. Number of multiplications required to evaluate a polynomial of 
degree n. See Todd [13]. 

Rapidly computable definition of determinant of A. (The definition in terms 
of n! terms is impossibly slow for large 1.) 

Evaluation of determinants by row elimination. Importance of choosing 
large pivotal elements. Methods of doing so conveniently with matrices of order 
1000 or more, stored on magnetic tape. 

Coding the calculation of the row rank of a set of row vectors. 

Computation of a+b on a computer, with minimum error. 

Actually solving linear algebraic systems. Round-off errors. 

Actually computing eigenvalues of matrices. The great distinction between 
symmetric and nonsymmetric matrices in this problem. 

Determination of whether a quadratic form is definite. 

Solution of linear inequalities. Duality theorems, and their use in computa- 
tion. 

“Finite calculus” of sums and differences. 

Convergence of linear iterative processes. 

Bounding the round-off errors in pseudoarithmetic operations; see Section 3 
above. 


Logic. Preparing a code to test the truth or falsity of a class of propositions. 
All sorts of questions relating computability to the theory of recursive functions. 


Advanced calculus. The nature of the remainder in Taylor’s theorem. 
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Interpolation theory and remainders. Convergence of interpolation poly- 
nomials. 

Direct methods in calculation of minimizing functions in calculus of varia- 
tions. 

Numerical summation of series with a computer. Code for computing 
Yon20f(m)x", where f(n) is a polynomial of low order, and | x| <i. 


Ordinary differential equations. Coding and running Euler’s method for solv- 
ing dy/dx =f(x, y). Bounding its error, and comparing with experience on test 
problems. 

Better methods for numerical integration of same differential equation. 

Two-point boundary-value problems. 

Computing the solution of an initial-value problem for y’’ +ay’+ay=0, 
where the characteristic equation nearly has a double root, and where the 
coefficients are 10-decimal numbers; see [4]. 


Probability and Statistics. Monte Carlo method for integration. Stratified 
sampling methods for greater accuracy. 

Generating normal random processes on a digital computer, with a pre- 
scribed power spectrum density. 

Estimating the power spectral density from an observed instance of a random 
process. 

Use of automatic computer to find mean and variance of given data which 
are very numerous. 


Number theory. Generation and testing of pseudorandom numbers for use 
in Monte Carlo experiments. 

Tests for primality of large numbers. 

All kinds of combinatoria! problems. 


Geometry. Iteration of rational transformations—theory and practice. 
Study of gradient methods on a computer. 


Fourier series. Numerical summation of a given series; acceleration of con- 
vergence. Orthogonal polynomials over finite point sets, and their use in curve 
fitting. 


8. Implementing the recommendations. Since it does not seem practical to 
teach an undergraduate course without following a textbook fairly closely, 
courses are likely to change only when textbooks change. There is no textbook 
for a modern integrated version of a pure mathematical field and its numerical 
application. Moreover, of all the numerical analysis texts published in the past 
10 years, only one seems really concerned with the problems of automatic digital 
computation. While desk calculators are far from dead, existing textbooks 
make it possible for any one to learn desk-calculator technique very well indeed. 
Future textbooks in numerical analysis must emphasize the use of automatic 
machinery. The one book to deal with automatic computation is the recent 
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Modern Computing Methods [2], regrettably anonymous but actually written 
by L. Fox, E. T. Goodwin, F. W. J. Olver, and J. H. Wilkinson. This is very 
clear, too short to be comprehensive, but full of annotated references to major 
sources. 

I therefore seem to be making a recommendation which is almost two gener- 
ations ahead of the current textbooks. The first generation might follow Modern 
Computing Methods and treat various aspects of computation by automatic 
computers. The second generation might integrate this material into basic 
existing mathematics courses. One can hope that some books will leap the whole 
step at once. 

There is one method of accelerating the appearance of textbooks in which 
mathematical theory and numerical applications are suitably intertwined— 
namely to write applied supplements to existing theoretical texts. For example, 
one might prepare a companion text to a textbook in the theory of linear alge- 
bra, in which numerical methods are discussed from the point of view of auto- 
matic computation, with suitable problems. Such a companion text might pre- 
sume that the student knows how to code, but that he is learning the theory of 
linear algebra for the first time from the basic textbook. The companion text 
would refer to the basic text for theory, and would indicate the appropriate 
stages at which the companion book should be used. 

Such companion texts would be useful for at least the following subjects: 
calculus, differential equations, statistics, linear algebra, advanced calculus. It 
would be well to include suitable analytical applications, as well as numerical 
ones. 

My basic recommendation for integrated courses in mathematical theory and 


‘application may have to await the appearance of suitable integrated texts or 


companion texts. What can be done meanwhile, aside from urging one’s friends 
to write expositions? 

For one thing, the coding course can be initiated at once anywhere. One 
needs only a textbook, say Andree [1] or McCracken [10], plus coding manuals 
which the International Business Machines Corporation or other computer 
manufacturers would probably be glad to furnish through their local represen- 
tatives. McCracken’s text deals with a mythical machine. 

No matter what courses are given, I believe all students of mathematics 
should have easy access to semiautomatic or fully automatic desk calculating 
machines, and some problems to solve with their aid. Even such an apparently 
simple problem as using an 8-bank calculator to find close upper and lower 


bounds for 
1.56898x? — 0.45678x? + 9.99938x — 0.42057 for x = 4/2 


will give students a surprising amount of insight into some problems of numeri- 
cal calculation. Some experience with desk calculators is almost an indispensable 
prerequisite to use of an automatic computer. 

If there is an automatic computer, students should be encouraged to use it 
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to find upper and lower bounds for various quantities which they deal with in 
their several theoretical courses. Such exercises should develop a good deal of 
power in numerical analysis. 

Finally, until integrated texts are available, it will probably be profitable to 
offer a specialized course in numerical analysis to undergraduates, using as a 
text such existing books as Modern Computing Methods [2], Hildebrand [7], 
Milne [11], or—better yet—selections from each. There are other good text- 
books also. Such a course should include as exercises, problems which I envision 
as later appearing in several theoretical courses. 


9. Acknowledgments. The author wishes to emphasize that the present 
thoughts are the result of informal conversations with many friends and col- 
leagues, who deserve credit for the good ideas but no blame for the author's 
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A CONNECTED TOPOLOGY FOR THE INTEGERS 
SOLOMON W. GOLOMB, Jet Propulsion Laboratory—California Institute of Technology 


A topology D for the positive integers is obtained when those arithmetic 
progressions {an+b} with (a, 6) =1 are taken as a basis for the open sets. They 
form a basis because the intersection of two such progressions is of the same type, 
or empty, as is easily verified. Note that every nonempty open set, being a 
union of arithmetic progressions, must be infinite. 

This topology furnishes an interesting proof of 


THEOREM 1. The number of primes is infinite. 


Proof. If p is prime, the progression {np} is closed, since its complement is 
{npt+1}U{np+2}U te U{np+(p—1)}, a union of open sets. Consider the 
union X =U, {np} extended over all primes. If this is a finite union of closed 
sets, then X is closed. But the complement of X is { 1}, which is neither empty 
nor infinite. Since the complement of X is not open, X cannot be closed, the 
union is not a finite one, and the number of primes is infinite. 

(A similar proof, in a stronger and very disconnected topology, was given by 
Furstenberg [2].) 


THEOREM 2. The topology D is Hausdorff. 


Proof. Given distinct positive integers s and t, choose a prime p (by Theorem 
1) which exceeds max (s, #). Then {pn+s} and {pn+t} are disjoint open sets 
which separate s and ¢. 


THEOREM 3. The topology D is connected. 


Proof. Suppose the integers could be represented as the union of two disjoint 
nonempty open sets 0; and QO». Let {aim+bi} be a basis set in Q;, and let 
{axn+b2} be a basis set in 02. Let a be a multiple of a;. If a were in O2, we would 
have a=Amno+B, where {An+B}(CO:. Since (A, B) =1, we would have (a, A) 
=1, and hence (a, A) =1. But then {a:n+b,} and {An+B} would intersect 
infinitely often, contradicting disjointness of O, and O,. Thus all multiples of 
a must belong to O,. Similarly the multiples of a2 must belong to O2. But then 
the common multiples of a, and a2 must belong to both O,; and O2, which contra- 
dicts disjointness. 


The author has recently learned that a proof of the connectedness of the 
topology D, without reference to number theory, was presented by Morton 
Brown at the April 1953 meeting of the American Mathematical Society in 
New York [1]. 


THEOREM 4. The topology D is not regular. 


Proof. Suppose that open coverings are given for the closed set {2n} and 
for the point {1} outside it. Any open covering of {1} not intersecting {2n} 
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must include a progression {en+1}, where e is even. That is, e€{2n}. Let 
{an+b} be the member of the open covering of {2n} which contains e, so that 
e=anot+b. Since (a, b)=1, we have (a, e)=1, whereby {an+b} intersects 
{en+1} infinitely often. Thus the closed set {2n} and the point {1} cannot 
have disjoint open neighborhoods. 


THEOREM 5. The topology D is not compact. 


Proof. The union U, {mp—1} extended over all primes is an infinite open 
covering for the positive integers. Since the omission of any progression {ng— 1} 
leaves the number g—1 uncovered, the Heine-Borel property fails. 

Actually, the topology D is not even locally compact, because every locally 
compact Hausdorff space is regular. For a proof of this, as well as for the more 
basic definitions of point-set topology, the reader is referred to [5]. 

Dirichlet’s theorem, which asserts that every progression {an+b} with 
(a, 6) =1 contains infinitely many primes, has an elegant formulation in terms 
of the topology D. 


THEOREM 6. Dirichlet’s theorem is equivalent to the assertion that the primes are 
a dense subset of the integers in the topology D. 


Proof. Assume first the validity of Dirichlet’s theorem. Then every non- 
empty open set contains primes, so that the primes are a dense subset of the 
integers. Conversely, assume that the primes are a dense subset. Then every 
nonempty open set, and in particular all the progressions {an+5} with (a, 6) =1, 
must contain primes. It is well known [4] that if every such progression contains 
at least one prime, then every such progression contains infinitely many primes. 
(In topological terminology: “If the closure of the primes is the integers, then 
the derived set of the primes is the integers.”) 

It appears unlikely that a complete topological proof of Dirichlet’s theorem 
can be given along these lines without the introduction of powerful new ideas 
and methods. 


Another familiar fact capable of topological formulation is 
THEOREM 7. In the topology D, the interior of the set of primes is empty. 


Proof. If there were an open set consisting entirely of primes, there would 
be a progression {an+b} with 1<b<a consisting entirely of-primes. But with 
my =a+b+1, ano +b=(a+b)(a+1), which is composite. 


It is interesting to consider also the topology D’ for the positive integers, 
which has as a basis those progressions {an+b} with (a, 6) =1 for all n>N. 
(Here N is allowed to assume all values.) This topology is clearly stronger than 
D, although Theorems 1 to 7 are still valid in D’. Moreover, certain theorems 
related to Eratosthenes’ sieve are available in D’. In particular, 


THEOREM 8. The set of positive integers m such that 6m—1 and 6m-+1 area 
pair of “prime twins” is closed in D’. 
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Proof. It is known [3] that the numbers m in question are precisely those 
positive integers not expressible in the form 6¢b+a+6 for any a21 and b21. 
Thus the complement of our set is Us21 (66+ 1)a+d}, where each progression is 
restricted to a21, and is open because (60+1, 6) =1. The union is open in D’, 
because it is a union of open sets. Thus the integers m for which 6m—1 and 
6m-+1 are both prime form a closed set. 
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AN APPLICATION OF TURNS AND SLIDES TO 
SPHERICAL GEOMETRY 


J. M. FELD, Queens College 


1. Introduction. The ordinary or Euler spherical triangle, the sides and 
angles of which are assumed to lie between 0 and 7, is determined uniquely 
by three points on a sphere, provided no two of them are diametrically opposite. 
Since each side of an Euler triangle has two poles, the choice of a polar triangle 
is ambiguous; moreover, more than one Euler triangle may exist with parts 
satisfying the law of cosines on a sphere. Moebius [3], consequently, generalized 
the concept of a spherical triangle in order to remove these ambiguities by orient- 
ing the sides and the angles and allowing them to assume all real values between 
0 and 27. Moebius’s device, though adequate for a satisfactory interpretation 
of the formulas involving the sines and cosines of the sides and angles, did not 
go far enough. Such formulas Felix Klein [2] called formulas of the first kind 
to distinguish them from others, essentially different, connecting the sines and 
cosines of the half-sides and half-angles, which he called formulas of the second 
kind. Among the latter there are, for example, the elegant formulas first pub- 
lished by Delambre (1807) and later independently by Mollweide (1808) and 
by Gauss (1809). What is remarkable about the Delambre formulas is that they 
carry a double sign and that there exist two classes of spherical triangle, one for 
which the upper sign is valid and another for which the lower sign is. It was 
Study [4] who grasped the full significance of these formulas; he showed that to 
establish their full import it was necessary to extend the range of values of the 
sides and angles of a spherical triangle to include all real values. We intend in 
this paper to prove Delambre’s formulas, and analogous ones for spherical 
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polygons, by means of simple applications of some elementary transformations of 
oriented lineal elements on a sphere, which, besides being intrinsically interest- 
ing, provide a procedure for obtaining more general results and, in addition, 
serve to make more natural Study’s concept of a spherical triangle. The trans- 
formations we have in mind are special cases of one called the spherical whirl [1]. 


2. Turns and slides. Let the unit sphere S have its center at O, the origin 
of a right-handed orthogonal Cartesian frame §o. Let an oriented lineal element* 
& be tangent to S at the point P; then we call the great circle passing through 
P tangent to & and oriented in the same sense as &, the great cycle of & We as- 
sume all lineal elements are oriented. Let the lineal element & have its point at 
(1, 0, 0) and let it be directed so that its great cycle passes through (0, 1, 0) 
and has for its positive orientation the counterclockwise sense when viewed 
from the point (0, 0, 1). We shall call & the primitive lineal element on S and 
Fo, its associated frame. 

Let €1, é2, és be the Hamilton quaternion units such that = ere, 
= —1, =e;, and let 


= Xoo + + Xee2 + = — — Xee2 — 


Then a rotation of S around an arbitrary diameter is given by the Hamilton- 
Cayley formula 


N(x)u* = £ux, N(x) = 


where u and u* are unit vectors (pure quaternions) with their initial points 
at O. The components x; of x are the homogeneous Euler parameters of the rota- 
tion. Let & be an arbitrary lineal element on S’, a unit sphere concentric with S, 
and let x be the quaternion of the rotation of S such that &)—&; then we let, as 
we may, the components of x serve as the homogeneous coordinates of &. For 
convenience, where no confusion will result, we shall let x designate both the 
rotation &)— 6 and &. Evidently, if x designates &, so does kx, where k is a non- 
zero scalar. If N(x) =1, x is said to be normalized. To any quaternion x, N(x) #0, 
correspond two normalized ones, differing only in sign, namely +x/[N(x) ]"2. 
We assume from here on that all quaternions are normalized. 

If the rotation x revolves S around the unit vector v through the angle 28, 
we can set x= —cos 0+ sin 9, so that its associated normalized quaternion 


—x = — cos (x — 6) + dsin (x — 8), i= 


designates a rotation of angle 27 — 20 around @. The rotation x also rotates F into 
another Cartesian frame § situated relative to & as §o is situated relative to &0; 
we shall call $ the frame associated with &. The point (0, 0, 1) of ¥ will be called 
the pole of the great cycle of &. 

A lineal element transformation x—x* given by 


* A lineal element on a sphere S is the combination of a point P on S and a line ¢ tangent to S 
at P. The lineal element is said to be oriented if t is oriented. 
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(2.1) = x4, 


where a= —cos a+u sin a, u?=—1, represents a lineal element rotation by 
means of which all lineal elements on S are rotated through an angle 2a around 
an oriented diameter having the direction of the vector u. Now, let us suppose 
that every lineal element x is rotated through the same angle 2a, but in such a 
way that each lineal element revolves about an oriented axis situated relative 
to its associated frame § as the axis u is situated relative to the primitive frame 
§. Such a transformation we have called a spherical whirl [1]. Evidently, if the 
whirl transports x—>x*, x*%=(x#)a, so that the whirl has the equation 


(2.2) x* = ax. 


The whirls (2.2) constitute a three-parameter group %®,, skew-isomorphic to 
the group of lineal element rotations (2.1). We shall be particularly concerned 
with certain one-parameter subgroups of Qs, namely, the turns, the slides and 
the dilatations. 

A turn T, is a whirl such that every lineal element is rotated around its point 
through an angle r. The group of turns is given by 


(2.3) T,: x* = [—cos (1/2) + sin (r/2)]x. 

A slide S, is a whirl such that every lineal element slides tangentially along 
its great cycle through the arc-length ¢. The equation of the group of slides is 
(2.4) Se: x* = [—cos (¢/2) + es sin (0/2) ]x. 

It will be observed that the identity turn and identity slide are Tx. and 
Sur, respectively, where k is an integer, and that the identity transformation is 
(2.5) I: 2* = px, p = +1, 

p taking the upper sign when k is odd and the lower one when k is even. 

A dilatation D, is a whirl such that every lineal element x is transported to a 
position orthogonal to a great cycle c which is itself orthogonal to the great 
cycle of x, the point of x traversing an arc of length 6 on c. The dilatations form 


a one-parameter group. The sense of c is chosen so that the equation of this 
group is given by 


(2.6) Ds: x* = [—cos (6/2) + ee sin (6/2) ]x. 


The dilatations, unlike the turns and slides, are contact transformations. 
Evidently, 


Ds = = Srj2T 
and 


= SS. = = I. 


We define a lineal element polarity ® by means of the equation 
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(2.7) x* = 2%; 

clearly, we can express @ also in the equivalent forms: 

(2.8) @ = = = = 

With the aid of (2.3) and (2.6) we obtain for ® the equation 
@: x* = + 


from this, or from (2.8), we get ®?x = —x; that is, ®@?=J7. Also by means of these 
equations we see that 


(2.9) S, = = 


As the lineal element x slides along its great cycle c a distance ¢, its polar, 
@x, turns around its point Pi, which is the pole of a, through an angle o; and, 
as x rotates around its point P; through an angle a, its polar, @x, slides on its 
great cycle, of which P is the pole, a distance ¢. 

Let x and y be two lineal elements on S. If x and y have a common point 
the whirl x—y is a turn; if x and y have a common great cycle the whirl xy 
is a slide. If they have neither a common point nor a common great cycle, let 
P, and cs be the point and great cycle respectively of x; P3 and @ the point and 
great cycle of y. Let P: be a point of intersection of cs; and c, and let c, be a great 
cycle that passes through P; and P;. Let us now subject x to a sequence of slides 
on the great cycles c; and turns at the points Pj, as follows: 


(1) a slide S, on ¢s, starting from P; and terminating at Ps, the slide includ- 
ing possibly any number of circuits of cs; 

(2) a turn Ts at P2, starting from the direction of cz; and terminating in a 
direction properly tangent to 4, permitting the mobile lineal element to com- 
plete any number of revolutions during the turn; 

(3) aslide S, from P: to P3; 

(4) aturn T, from to ¢2; 

(5) a slide S, from P; to Pi; and 

(6) a turn T, from c to ¢c3. With T,, the lineal element is brought back to 
its initial position. The complete circuit is given by 


(2.10) TST = ps p= 1. 


3. Spherical polygons. We define a Study spherical triangle by means of a 
circuit (2.10); the slides determine its sides and the turns its angles. Degenerate 
triangles are not excluded. The first three steps in the circuit (2.10) produce the 
whirl x—y. Multiplying the two members of (2.10) on the right by # and on the 
left by T_,S.4T_., we obtain two equivalent forms for the whirl, namely, 


(3.1) SeTpSe = pT-~S-sT_a- 


By means of great cycles passing through P; and P; and orthogonal to q and ¢, 
respectively, the whirl (3.1) can be resolved also into products such as TDS, 
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SDT, etc., where T, D and S represent suitable turns, dilatations and slides. 

Since any quaternion defines a whirl, it can be expressed in any one of these 
factored forms. Consequently, if we replace x in the rotation formula u* =Zux 
by one of the sets of three factors into which it can be resolved, we obtain for 
each set a parametric representation of the rotation, the parameters represent- 
ing the magnitudes of the slides, turns and dilatations that occur in the product. 
In particular, we get Euler’s parametric form when we let x = S,TsS-. 

Let x be a rotation of angle @ around an axis carrying the unit vector 
CoS Cos +e3 Cos v; then we can let 


x = — cos (6/2) + vsin (6/2) = TaSpS.. 


Substituting for the turns and slides from (2.3) and (2.4), we obtain, by equat- 
ing the corresponding components: 


a+c 
cos — = cos — cos , 
2 2 2 
cos A sin — = sin —cos , 
3.2) 2 2 
cos wp sin — = sin — sin ) 
2 2 2 
B . atc 
cos vy sin — = cos — sin . 
2 2 


If we multiply both members of (2.10) on the left by @ and insert PG =I 
between the factors in the left-hand member, we obtain: 


(PT 2) (PSP) = pPx; 


hence, by virtue of (2.9), we see that as x makes a circuit of the Study triangle 
(2.10), its polar, ®x, makes a circuit of its polar triangle. 

Inserting in (3.1) the formulas for the turns and slides given by (2.3) and 
(2.4), we obtain the formulas of Delambre: 


B at+ec b at+y 
cos — cos = — pcos— cos , 
2 2 2 2 
B a-c¢ aty 
sin — cos = pcos— sin ) 

(3.3) ( + 1) 
p= 
sin — sin = — psin— sin , 

2 2 2 2 
B . e+e 
cos — sin = psin— cos , 
2 2 
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It will be observed that the left members of the equations (3.3) appear also 
in equations (3.2); consequently, with every Study triangle there is associated 
an orthogonal transformation, the axis and angle of which are given by equations 
(3.2). By cyclic permutation we get two more sets of equations like (3.3). 

There are Study triangles for which only those Delambre formulas are valid 
in which p= +1, and also others for which the Delambre formulas with p= —1 
are valid. The full significance of the double sign was made clear by Study 
[4, 5]. In order to formulate Study’s principal theorem bearing on this circum- 
stance, let us, with Study, call the triangles belonging to the class in which 
p=+1 proper, and those belonging to the class in which p= —1 improper. 
Furthermore, let us consider the six-space Sg, the points of which have for their 
coordinates the quantities a, b, c, a, B, y, all of which are unrestricted real num- 
bers. Then Study’s theorem can be stated as follows: 


The set of points in Sg which are the images of the aggregate of Study triangles 
separates into just two continua (compact, connected point sets), one corresponding 
to the set of proper triangles and the other to the set of improper triangles. 


The Euler triangles belong to the set of proper triangles. 

In the same manner we can, of course, by means of a circuit of slides and 
turns, define a spherical polygon of any number of sides; the slides determine 
the sides in sense and magnitude while the turns, in like manner, determine the 
angles. The polarity (2.7) assigns to every polygon a unique polar polygon. The 
sides and angles are unrestricted in magnitude, as in the case of Study triangles. 
Starting from the equation for a circuit of slides and turns for such a polygon, 
like equation (2.10) for triangles, we can derive sets of formulas connecting its 
sides and angles, analogous to those of Delambre. Let it suffice to indicate the 
results obtainable for a spherical quadrilateral determined by the following se- 
quence of slides and turns: 


(3.4) TSaT = pt, p=tl. 


Multiplying both of the members of (3.4) on the right by #, and on the left by 
S_.T_,S_aT_s, we obtain 


Substituting in (3.5) the formulas for the slides and turns given by (2.3) and 
(2.4), and simplifying the results, we obtain the following four formulas: 


[ a+b [ a+ 8B a-—b 
cos + cos cos ——— + | cos — cos cos 
2 2 2 2 


y¥+6 y¥—6 c+d 
= {| cos + cos 00s 


y¥+6 
+ | cos — cos cos F 
2 2 2 


2 


4 


by 


nd 
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[sin + sin [00s sin — sin | 00s 
2 2 2 2 2 2 


= p4| sin sin cos 
E y¥+é6 c+d 
+ | sin — sin cos , 
2 2 
E + sin ] sin + [sin — sin | sin 
2 2 2 2 2 2 
6 —6 —d 
+ | sin — sin sin ‘ 
2 
cos + cos sin + | cos — cos sin 
2 z= 2 2 2 2 
{| 
=— c cos sin 
2 


6 
+ [cos 


Cyclic permutations of the sides and angles provide three more sets of for- 
mulas. It will be observed that these formulas are self-polar; moreover, as in the 
case of Study triangles, there are two kinds of quadrilaterals: those whose sides 
and angles satisfy the formulas with p= +1, and those that satisfy the set with 

=—1. 

In [1] we defined a spherical turbine as a series of oriented lineal elements, 
the points of which lie on a cycle c, and the great cycles of which make a con- 
stant angle a with c. When a=0 we have a tangential turbine, and, when c is a 
point cycle, a point turbine. The equation in lineal element coordinates of the 
point turbine at the point e is clearly ex =e. By means of a lineal element ro- 
tation (2.1), this is transformed into a point turbine at an arbitrary point having 
for its cordinates the components of the unit vector r, namely, Ze.x =r. The dila- 
tation D,; transforms this point turbine into a tangential one with center at r 
and radius equal to 6, given by the equation 


(3.6) &(e, cos 6 — sin 6)x = ¢. 


When 6=7/2, (3.6) is the equation of a great cycle. 
If we let the four components x; of the quaternion x be the homogeneous 
coordinates of a point in projective three-space S;, we obtain a continuous one- 
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to-one mapping of the lineal elements x on S upon the points x of S;. By virtue 
of this mapping, the ® lineal element whirl-rotations x—axb, are mapped on 
the «*® displacements in elliptic three-space E;; indeed, to the whirls x—ax 
correspond the left translations in E3, and to the rotations x—>xb correspond the 
right translations [1]. To the turbines on S correspond one-to-one the lines in 
E;. To a slide of the lineal elements on S corresponds a left translation in which 
the points of E; glide on lines belonging to a congruence of Clifford parallels 
C,, and to a turn on S corresponds a left translation on the lines of another 
congruence of Clifford parallels C,. The sequence of slides and turns in (2.10), 
defining a triangle on S, corresponds to a sequence of left translations that a 
point in E; executes in describing a complete circuit of a skew hexagon, the 
sides of which belong alternately to C, and C;. Similarly, to a spherical polygon 
of m sides defined by slides and turns, there corresponds in E; a polygon of 2n 
sides that lie alternately on lines belonging to Ci and C2. If, however, we let a 
point in EZ; make a circuit, by means of left translations, of a polygon, say a 
hexagon, the sides of which are not restricted to C, and C2, the corresponding 
lineal element on S executes a sequence of six whirls in each of which an “arc” 
of a turbine is traversed; the figure on S that corresponds to a general hexagon 
in Es is consequently a closed chain of six “arcs” of turbines. For such a figure 
we could also obtain formulas analogous to those of Delambre; however, the 
lengths of the “arcs” would be given in an elliptic metric and the formulas would 
involve additional parameters, namely, the six radii of the turbines and the 
angles of inclination of the lineal elements in the turbines. 

An interesting special case arises when the six turbines are alternately tan- 
gential turbines with radii between 0 and 7/2, and point turbines. Such a figure 
is a “spherical triangle” t whose sides are oriented arcs of small cycles. Let us 
assume that the three small cycles have equal radii, and that the sides and angles 
of t are measured in the elliptic metric of their image in E;. It can be easily 
verified that the tangential whirl that makes all of the lineal elements of the 
cycle (3.6) advance along its periphery for the distance a/2 is given by 


(3.7) Wa: x— qx, q = — cos (a/2) + (e: cos 6 — es sin 8) sin (a/2). 
It should be observed that 
W, = D_sT.Ds, W..= 8 = 2/2 8. 
We can now define the sides and angles of t by means of the following circuit: 
(3.8) = px, p= +1, 


where a/2, b/2, c/2, are the lengths of the sides of t, and a/2, 8/2, y/2, are the 
magnitudes of its angles. From (3.8) we obtain 


Substituting for the turns and tangential whirls in (3.9) the expressions given 
by (2.3) and (3.7), we obtain the following formulas: 
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cos 6 sin — sin — cos — cos 
2 2 
c 
= p| cos sinc sin — cos — cos 
cos 6 cos — sin — cos 26 sin — sin — sin — + sin — cos — cos — 
2 2 2 2 2 2 2 
= — p| cos 6 sin — cos + cos— sin , 
2 2 
sin 26 sin sin sin — sin sin cos ry = psin 6 cos sin 


When 6=7/2, t becomes a Study triangle, and these three formulas reduce 
to Delambre’s. However, since the sides and angles are measured in an elliptic 
metric in which 7 is the magnitude of one complete revolution of x around its 
point, as well as of a complete circuit by x of its great cycle, we accordingly let 
a, b, c be the lengths of the sides and a, 8, y the magnitudes of the angles, to re- 
store the spherical metric. 
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AN INEQUALITY ARISING IN GENETICAL THEORY 


H. P. MULHOLLAND anp C. A. B. SMITH, University of Exeter and 
University College, London, England 


1. The study of genetics provides simple examples of the use of such mathe- 
matical tools as probability theory, calculus, and matrices: and it may well be 
of interest to teachers of mathematics and statistics. Recently one of us 
(C. A. B. S.) came across an unusual inequality in a genetical problem, and the 
other author (H. P. M.) succeeded in proving its correctness. The inequality is 
as follows: 
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Let v, u be arbitrary nonzero real kX1 and kXk matrices respectively, with all 
their elements nonnegative, and with u symmetrical, 1.e., u" =u, where u™ denotes 
the transpose of u. Then, for any positive integer n, 


(1) (v?u"v)(v7v)"" (v7) uv" 
with equality if and only if v is a latent vector of u. 


We have received alternative proofs of the first part of this result from 
M. M. Crum and R. Scheuer, and a weaker result has been published by 
S. P. M. Mandel and I. M. Hughes [4] since this paper was first written. We 
explain here why this inequality is of interest, as well as giving a proof. 


2. One of the most striking of biological phenomena is the great variability 
of many species of plants and animals, even when raised in a nearly constant 
environment. For example, some men are short, others tall, some are clever, 
others stupid, and some dark, others fair. It is far from being immediately obvi- 
ous why this should be so. For if natural selection leads to the survival of the 
fittest, why has humanity not been reduced to a single uniform type of maximum 
fitness? It would be very interesting to know the reason; but it would also be 
very presumptuous at present to attempt a detailed explanation for qualities 
such as height, intelligence, or hair color, since our knowledge of how these are 
inherited is very incomplete. Instead, we here consider a rather simplified 
model. In fact, in what follows both the biological facts and the theoretical 
model have been, strictly speaking, oversimplified, but with the essential features 
retained. Those who wish to pursue the matter further can read C. C. Li [3] 
or O. Kempthorne [2], who give references to other work on the subject. 

Many properties of living organisms are determined by objects called genes. 
For example, there are three types of genes, usually called A, B and O, that 
determine the ordinary blood group in man; each individual carries just two 
such genes, so that there are 6 distinct blood groups, or pairs of genes, namely 
OO, AO, AA, BO, BB, AB. Their practical importance is that it is dangerous to 
transfuse blood from an individual carrying an A gene into one without an A 
gene, e.g., AO or AA blood may not be transfused into a BO or BB patient, and 
a similar rule holds for the B gene. (It can be deduced that when only transfu- 
sion is considered it is unnecessary to distinguish between the types AO and AA, 
which are lumped together as group A; similarly BO+BB=group B.) Mendel’s 
laws state that each child gets at random one of the two genes carried by its 
father, and independently and at random one of the two carried by its mother. 
Thus, if the father is AB and the mother OO the child is equally likely to get A 
or B but must get an O; it is therefore equally likely to be AO or BO. A child of 
an AB father and AO mother is equally likely to be AA, AO, AB, or BO. 

In general we consider a set of k types of genes or alleles, Gi, - - - , Gx, with 
the property that any individual carries just two such genes, possibly of the 
same type. In a large population of N individuals there will be a certain number 


$2 
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Ny, of type G,G,: the ratio N,,/N=p,, is the frequency of type G,G, in the 
population. The frequencies »,, describe the population as far as the genes G, 
are concerned. Thus in London approximately 7 per cent of persons have blood 
group AA, 35 per cent AO, 0.4 per cent BB, 8 per cent BO, 3 per cent AB, and 
47 per cent OO (blood group O). The population of N individuals can also be con- 
sidered as one of 2N genes. Thus, since every GiG; individual carries two Gi 
genes, and every GiG, (r~1) carries one, the number of G; genes is 


Ni = Nu. 
A division by 2N, the total number of genes, shows that the frequency of G; is 
(2) fi = Ni/2N = put dpi t+ dpu, 


and similarly for fe,- +--+, fz. Thus in London the A gene has a frequency 
7+3:-35+43-3=26 per cent. In a large population it is generally found that the 
frequencies are the same in males and females, and so equality is assumed in 
the following argument. 

Evolutionary genetics is the study of how a population gradually changes 
from generation to generation. For simplicity we here assume the generations 
to be separate, as with annual plants, which have one generation each year. 
bre(g), f(g) represent the frequencies in the gth generation. Now there is a saying 
that “marriage is a lottery.” If this is taken to mean that marriage partnerships 
are formed at random, it seems to be true as far as blood groups and similar 
qualities are concerned. For while one may choose a partner for charm, intel- 
ligence, beauty or common interests, it isn’t likely that his or her blood group 
will seem particularly adorable. We make the assumption that the genes under 
consideration do not affect the chances of having children: this means that a 
baby chosen at random gets at random one of the genes carried by a random 
father, and similarly and independently a random gene via the mother. Thus 
the chance of a baby of generation (g+1) being of type G,G, is f,(g) -f-(g), while 
the chance of being G,G, (r#s) is 


f(g) -fe(g) + = 2f-(g) -fa(2) 


since either the G, gene comes from the father and the G, from the mother or 
vice versa. But babies with different genes will have different liabilities to dis- 
ease or accident: we let w,, be the chance of a baby of type G,G, surviving to 
adult life. The matrix w of the w,, accordingly represents the effect of natural 
selection. The frequencies p,,(g+1) of the various types G,G, among adults of 
generation (g+1) will be in the ratios of the frequencies at birth multiplied by 
the survival probabilities, i.e., for the types GiGi, GiG2, GiGs, - + + they will be 
in the ratios 


= pulg 1): + 1): pis(g + 1): ---. 
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To simplify the notation we introduce here certain conventions. Unless other- 
wise stated p,,, f, stand for the frequencies p,,(g), f-(g) in the gth generation, and 
P,,, F, for the corresponding values in the following, (g-+1)th, generation. After 
a summation sign Greek letters will be used to denote the suffixes to be summed 
over, and these will be allowed to range over all possible values unless the con- 
trary is indicated: for example, }-f2=1, and (2) becomes f,=3prr+4 > Pra- 

Now it may happen trivially that all the products w,,f,f, are zero, so that 
there are no survivors in generation g+1; we suppose this case excluded. Since 
the frequencies P,,=,,(g+1) must add up to unity, and their ratios are given 
by (3), their actual values are 


Pre = Wrefe/W(f), Pra = ¥ 5), 


where W(f) stands for the quadratic form f’wf= )-weafafs>0. From (2) the 
gene frequencies in generation g+1 are 


(4) F, = fe Wrafa/W(f) = 3f,[W(f) 


which we write for brevity as F=¢(f). This is a recurrence relation connecting 
the gene frequencies in generations g and g+-1. It is easily seen that if w,.f,f,>0 
then w,,F,F,>0 so that, by assuming that the products w,,f,f, do not all vanish, 
we have ensured that the population will continue to exist in generation g+2, 
and the frequencies in that generation are obtained from those of generation 
g+1 by equations similar to (4). By induction, the population continues to exist 
in every subsequent generation g*>g, and W(f(g*))>0, and a relation like (4) 
connects every two successive generations. 

A real numerical example is provided by Allison [1]. In central Africa there 
are genes affecting the hemoglobin (officially called Hb4, Hb’, but here we use 
the letters G:, G2, for consistency) such that GiG; type individuals (the normal 
European type) are susceptible to malaria, while the G2G2 tend to die in child- 
hood from a severe sickle-cell type of anemia. Their survival chances are, ac- 
cording to Allison, approximately in the ratio wu! wi2!wW2=0.8:1:0.2, so that 
the recurrence relation is, by (4), 


(5) Fy = fi(wufi + wief2)/ + 2wifife + woof), 


where f.=1—f,. Thus if we take, say, fi(1) =.1 as initial value, the frequency 
of G; in successive generations will be fi(2) =.280, fi(3) =.464, fi(4) =.579, - +> 
tending eventually to .8, whereas if initially f,(1) =.9 the successive values are 
.889, .878, .868, etc., again tending to the limit .8. From this it seems plausible 
that whatever frequencies (other than 0 or 1) the genes G; and G, have initially, 
they will approximate to fi=.8, fe=.2 respectively after a few generations; 
fi=.8 is a stable equilibrium. Thus all three types GiGi, GiG2, G2G2 will continue 
to exist indefinitely in the population as a result of natural selection. Of course, 
if initially f; =0, it will always remain so, for in our model if a gene is not present 
in any one generation it cannot appear in any subsequent generation. (In nature 
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it could be introduced by migration or mutation.) We call an equilibrium de- 
generate if any of the f, are there zero, otherwise nondegenerate. 

An example of a stable equilibrium with three genes Gi, Gz, G3 is shown in 
Figure 1. Here the selection matrix is taken to be 


4 .7 .6 
w=|.7 .1 


The small circles show the values of the three gene frequencies fi, fe, fs plotted 
as trilinear coordinates, which is possible since fitfe+/:;=1, and the initial 
values are various points near the perimeter of the triangle. The sequences of 
points clearly appear to tend to an equilibrium at E = (3, 4, 4). 


(0,1,0) 


(0,0,1) (1,0,0) 
Fic. 1 


3. It is therefore an important genetical problem to know when such stable 
equilibria exist, and whether a population not initially in equilibrium must 
tend to one under the action of selection. 

When only two genes G; and G; are concerned, a complete analysis can be 
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made by elementary methods. The cases in which wiw22=0= wi are trivial in 
the sense that at least one gene is completely lost after the first generation. We 
therefore exclude these, as well as the cases for which wu fj = = Wafe=0, 
already excluded in Section 2. For other cases the analysis is most easily done 
in terms of the ratio r=fi/fe=fi/(1—fi), which has the range of valuesOXr<S ~; 
if R= F\/F, the recurrence relation (5) becomes 


(6) R= r(wir + W12)/(wier + Wee) 


(and R= « when r= ~). An equilibrium point is one at which R=r. From (6) 
we see that there are always two degenerate equilibria r=0 and © (at which 
only one gene is present); there is also a third, nondegenerate equilibrium r=q 
where 


(7) (wie — wee) / (wie — 


provided that 0<qg<~, unless we have wi =wi2=We, and g is indeterminate: 
but then R=r for all r, and any population is in a state of (neutral) equilibrium. 
Apart from this, the following cases may occur: 


(I) wiu<wi2>w2. A nondegenerate equilibrium g exists. Furthermore, it is 
stable, and whatever initial value other than 0 or ~ we take for 7, the popula- 
tion will eventually tend to the equilibrium g. For from (6) and (7) 


wWurt wee (we — wWy)r 
r—q Wir + We Wier + Wee 


Thus if we choose arbitrarily two positive numbers r’, r’’ such that r’<r<r’”, 
r'<q<r'’, we have 


(8) 0<(R-g)/(r-—g <1-¢, 


where ¢ = / >0. The distance |r—g| of r from the equi- 
librium q is diminished by a factor smaller than 1—c in one generation. It follows 
that r’<R<r’, the argument repeats in the next and every subsequent genera- 
tion, and the distance tends to zero with (at least) geometric rapidity. 


The other cases which follow can be covered by a similar type of argument, 
so we omit details: 


(II) wu >wi2<wz. There is an unstable equilibrium at g; from smaller values 
of 7 the sequence tends to 0, and from larger ones, to ~. 


(III) wuSw2Swe (but not wi=w2=we). For any initial r, except ~, 
successive values tend to 0, the only stable equilibrium. 


(IV) wu (but not wi=wi2=wWe). As case (III), with 0 and 
interchanged. 


4. In general, if f(1) is the initial gene-frequency vector, the vectors in sub- 
sequent generations are f(2) =¢(f(1)), (3) =@(f(2)), and so on, where the func- 
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tion ¢ is defined explicitly by (4); but it is difficult to get any explicit expression 
for f(z) for general 7, or study the properties of the sequence directly. However, 
much information can be got from the following two theorems. (From now on, 
it is to be always understood that a frequency vector, such as f, is to be re- 
stricted to lie in the simplex [all f;20, >°fa=1] and in the region where 
W(f)>0, unless the contrary is explicitly stated.) 


THEOREM 1. E is an equilibrium point if and only if either (trivially) only one 
E;+0 or the function W(f) has a stationary point at f = E when variation is allowed 
only in the components f; for those i for which E;~0; if E;=0, f; is held fixed at 0. 
In particular, a nondegenerate equilibrium E (i.e., one with all E;>0), is a station- 
< point for W(f) with the f; varying freely (apart from the necessary restriction 

fa=1). 


Proof. Suppose, for example, that E,=E,=0, E;>0 (¢>2). The necessary 
and sufficient conditions for W(f) to be stationary at E subject to the conditions 
fi=fe=0, >> fa=1 are that for some 


(9) OW (F)/dfs = = = = d 


when f=E. By Euler’s Theorem }>f.0W(f)/0fa=2W(f), i.e., \=}W(E). Thus 
(9) is equivalent to 


(10) fe = 
for all 14, when f=E: by (4) this means that $(E£) = E, i.e., E is by definition an 
equilibrium. 


COROLLARY. An equilibrium point E is either isolated or is on at least one 
straight segment on which every point f is an equilibrium point with W(f) = W(E). 


Let us consider, for example, those points f for which 
Dhe=1, 


where, for the moment, we do not require every f; to be nonnegative. Since these 
equations are linear in the components f; they are satisfied either at no point, 
, or at one point, or at every point of some r-dimensional linear subspace with 
1<r<k-3. As in the proof of the theorem, on such a subspace W(f) is station- 
ary at every point for variation of f within the subspace, and thus W(f) must 
clearly be constant on the subspace. On imposing the restrictions f; 20 and using 
the theorem we see that the set of equilibrium points E of the special type with 
fi=fe=0 is either empty, or contains only one point, or consists of an r-dimen- 
sional convex body (with 1Sr<k—3) on which W(f) is constant. Similarly for 
equilibrium points of the other special types. Since the number of possible types 
is finite (actually 2*—1), it is clear that an equilibrium point must either be 
isolated or belong to a convex body (corresponding to one of the special types) 
in which every point is an equilibrium point and W(f) is constant. This suffices 
to establish the assertion in the corollary. 


680 AN INEQUALITY ARISING IN GENETICAL THEORY [October 


THEOREM 2. W(F) = W(9(f)) = W(f), with equality if and only if f is an equi- 
librium. 


(We recall that values of f at which W(f) =0 have been explicitly excluded: 
they make F indeterminate.) 


Proof. (We assume the correctness of (1), which will be proved later.) Write 
Vf; and uy Then W(f) = > and, from (4), 


= uv)? = 


Since v’v= =) fa=1, the inequality W(F)2W/(f) is equivalent to 
(v7u*v)(v7v)?= (v7 uv)*. Furthermore from (9) and (10) the condition for equi- 
librium is that for each i either f;=0? =0 or 2 )-wiefa=X, that is, for every i, 


= AD, OF Uv = Fry, 


1.e., that v is a latent vector of u. This reduces Theorem 2 to the case n =3 of 


(1). 
5. Now Theorem 2 implies that 
(11) W(f(1)) S W(4(2)) S W(F(3)) S---. 


If the function W(f) is represented by its graph, a “surface” in k dimensions, 
then the sequence { f(r) } is represented by a sequence of points climbing up the 
surface. It is natural to suppose that they will eventually tend to a “peak,” 
which will be a point of stable equilibrium. This picture of the genetic behavior 
of populations is due to Sewall Wright [5], but his formal proofs do not seem to 
be complete. In any case, unless we impose further restrictions upon the func- 
tion W(f) the situation may be complicated by the presence on the surface of 
saddle points, and also of level ridges or plateaus: the sequence of points may 
converge to a saddle point, which will be a point of unstable equilibrium, and 
(although we conjecture the contrary) it might have several limit points instead 
of converging to a unique limit. However, with the aid of Theorems 1 and 2, 
we are able to demonstrate the following theorem, and we shall leave till Section 
6 the discussion of various genetically appropriate restrictions. 


THEOREM 3. (1) Any limit point of the sequence { f(r) } is an equilibrium point. 
(II) Let W(f) have a strict maximum at E (for admissible values of f): then E is 
a stable equilibrium point, i.e., if amy point £(i) is near enough to E then f(r) —-E 
as (III) If Eis an isolated equilibrium point, then either f(r) as r> 
or E is not among the limit points of { f(r) }. 


Proof. Since the sequence { W(f(r))} is bounded and, by (11), is monotone, 


it must have a unique limit W., say, as ro, and W,,= W(f(1))>0. Also, 
f(r) lies in the bounded closed region 
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R: Dfa=1, W(f) = W(F(1)). 


Hence the sequence {f(r)} must have at least one limit point L in R: for any 
such L we have W(L) = W,,>0, since W(f) is a continuous function of f. Since 
W(L)>0 we have ¢(f) continuous at L. Hence {d(f(r))} must have $(L) as a 
limit point. But this sequence is {f(r+1)}, and so $(L) is a limit point of 
{ #(r)}. Hence, as before, W(¢(L)) = W... Thus W(@(L)) = W(L); whence, by 
Theorem 2, L is an equilibrium point. We have now proved (I), and also that 
W(L) has the same value for every limit point L of { f(r) I. 

For Part II, since W(E£) must clearly be positive, we can surround E by two 
regions S;, S2 consisting of the points f in R at which the distance d of f from E 
satisfies the respective inequalities d<p,, d<p2, where p; and p2 are such that 
pi>p2>0 but otherwise are at our disposal. Now W(f) has a strict maximum at 
E, and so, by Theorem 1, Corollary, E is an isolated equilibrium point. Hence we 
can choose p; to b° such that 

(i) if f is in S, but not at E, then W(f) < W(E); 

(ii) S, contains no equilibrium point other than E. 

Further, since W(E) >0 the function $(f) is continuous at E: hence we can now 
choose p2 to be such that 

(iii) if f is in S, then (Ff) is in Si. 

Now S,—S: is a bounded closed region and does not contain E: hence the great- 
est value m taken by W(f) in S,—S: must be such that m< W(E). 

Suppose now that f(z) is in S; and near enough to E to make W(f(1)) >m. 
Then, by (11), W(f(r)) >m for r=i+1, 1+2, +--+, and so f(r) can never lie in 
Si—S: if r2i. However, f(z) is in S; and, by (iii), as r runs through 7+1, 
i+2,+-- the first point f(r) (if any) to fall outside S; would be preceded by a 
point in S,—S2. But this cannot occur: hence f(r) cannot escape from S;. Since 
S; is bounded and closed, every limit point of { £(r)} is thus in S;. By Part I 
and condition (ii) it follows that E is the sole limit point. Altogether, if f(7) is 
near enough to E, then f(r) —>E as r->~, as required. 

For Part III let us suppose that E is a limit point of { f(r) } : we can then show 
that E is the unique limit of f(r) by modifying the argument used for Part II. 
We have to drop the requirement (i), and consequently cannot claim that 
f(r) is never in S;—S: for r=%. However, we observe instead that it can be in 
S:—S: at most finitely often, since {f(r)} has no limit point in S;—S». Since E 
is a limit point of {f(r)} there are infinitely many points f(i) in S: and (thus 
also in S:), and, as before, for every departure of f(r) from S; there is a point 
of the sequence in S,—S:. Hence there can only be finitely many departures, 
and f(r) must be confined to S; when r is large enough. It follows as before that 
f(r) E as 


6. We remark that better results than those in Theorem 3 can be obtained 
under various restrictions. Firstly, when there is a nondegenerate equilibrium 
E at which W(f) has a strict maximum it is not difficult to show that this is the 
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only local maximum point: for any section of the “W-surface” by a plane 
passing through E will be parabolic with the vertex of the parabola at E. Thus 
E must be the only stable equilibrium. In other cases there can be more than 
one stable equilibrium, necessarily degenerate: a simple case is 


W(fay fas fa) = — 2fs)? + + fo + = [3 + (fr — 2f2)*]/10, 


with local maxima at (1, 0, 0) and (0, 1, 0). 

Secondly, by Part III of Theorem 3, { f(r)} must converge to a unique limit in 
all cases in which there are only isolated equilibria (or, equivalently, in which there 
is only a finite number of equilibria). From the corollary to Theorem 1 and from 
the conditions for an equilibrium point [see (9) and (10) above] it is clear 
that every equilibrium point will be isolated unless one or more of certain 
determinants, formed from the numbers wes, happen to vanish. In this sense 
any cases in which { f(r)} might not converge are exceptional, and we should be 
unlucky if we encountered such a case in a real situation. 

Further, although in our deterministic model convergence to a “saddle 
point” equilibrium E is possible, yet such an equilibrium is unstable, and in 
practice random fluctuations would almost certainly lead the system to some 
point f’, say, with W(f’) appreciably greater than W(E), and after f’, since 
W(f) must increase, f would be kept away from the neighborhood of E. In view 
of all this Sewall Wright’s picture of the process [see Section 5] is “practically” 
correct. 


7. It only remains to prove the inequality (1). If we write 1,=v7u"v, and 
exclude the trivial case v=0, the inequality can be written 


(12) ¥(v) = Mn/Mo — (Mi/Mo)" > 0 


under the conditions (qa) all v;20, (8) all u;;=u;;20, (vy) v is not a latent vector 
of u. We know that there is a real orthogonal matrix, R, say, which reduces u 
to principal axes, 7.e., R’uR=t where t;;=0 when 7+j. Let a= Rv; 
then 


(13) M, = Vvuveata= 
where by condition (y) a is not a latent vector of f, #.e., there exist 7, 7, such that 
(14) a; 40, big ty. 


Now M, is the mth moment of a distribution of masses aj at respective dis- 
tances #;; from the origin. When all ¢;;2=0, (12) is a standard result, true even if 
(a) and/or (8) do not hold. For on allowing to take all real values, a straight- 
forward differentiation and use of Cauchy’s inequality show that (d/dm)? In M, 
> 0, t.e. the graph of In M, as a function of m is concave upwards. Hence of two 
adjacent chords of the graph the one on the right must have the greater slope; 
if h<m<n, 
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n—m m—h 
or 
(15) 


and (12) is the particular case m=1, h=0. (A fortiori (15) is true if m and h are 
even, even if not all ¢;;20.) However, if some ¢;; are negative and m is odd (15) 
need not necessarily hold, even if (a) and (8) are satisfied. For example, take 
u= HAE then M,= M;=4, and M3M, <M? and 
< 2: 

We now establish the truth of (12) under conditions (a), (8) and (vy) by 
reductio ad absurdum. Suppose it untrue; choose u, v of the smallest possible 
dimensions k to satisfy (a), (8) and (vy) but not (12). Clearly k>1; and we can 
suppose all v;>0, for if any v;=0, then on suppressing the 7th row and column 
of u we have a fortiori a contradiction of (12) in a smaller number of dimensions. 
Thus (v) depends only on the direction of the vector v, is a continuous function 
of v, is strictly positive on the boundary of the set of v’s that satisfy (a), where 
some v;=0, and is negative or zero at some point v* that lies inside this set and 
also satisfies (vy). Thus among points satisfying (a) there is at least one point 
V at which the function ¥(v) takes its minimum value m <0 and is stationary. 
If m<0 such a point V must satisfy (y), since it is easy to see that y(v) =0 if v 
is a latent vector of u. If m=0 we may clearly take V to be v*, and again V 
will satisfy (vy). Write ¥(v) =w(a), where a= R’v, and write A= RV: then at 
a= A the condition (14) is satisfied and w(a) is stationary, so that dw(a)/da;=0 
at that point. But, by (13), 0M,/da;=2a,f%, so that for each 7 either A;=0 or 
Moly — Ma—n(Mi/ Mo)*"(Motii— Mi) =0 at a=A. In the latter case tig is a 
root of the equation 


(16) — ni(Mi/Mo)""* + — 1)(Mi/Mo)" — ¥(V) = 0. 


The left-hand side of (16), considered as a function of ¢ only, has a unique 
minimum for nonnegative ¢ when t= M,/M, and its value is then —¥(V)20: 
thus the only possible nonnegative root of (16) is ¢= Mi/ Mo, and so t;;= Mi/Mo 
whenever A;¥~0. But A?’tA/ATA=Mi/Mo, and so 
which can accordingly only be satisfied if t;;= Mi/Mo whenever A;+0. This con- 
tradicts (14) for a= A and so the theorem is established. 
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SEPARABILITY CONDITIONS FOR SOME SELF-ADJOINT 
PARTIAL DIFFERENTIAL EQUATIONS 


C. P. WELLS, Michigan State University 


1. Introduction. Let L., be a linear partial differential operator with inde- 
pendent variables x and y. Denote by T a real transformation x=x(u, v), 
y=y(u, v) with finite, nonvanishing Jacobian. Let the transformed operator be 
given by TL,,=Lue. If Lu» can be written in the form 


(1) F(u, v)(Lu + Ly), 


where L, and L, involve u alone and » alone, respectively, the operator Lz, and 
the equation L.,[y]=0, are said to be separable in the coordinates u, v. We 
shall call u, v a separable system of coordinates if Lzy is separable in u, v. 
Separability conditions and separable systems of coordinates for Laplace’s 
equation, the Helmholtz equation and the Schrédinger equation are well 
known. (See, for example, [1] to [4].) While these constitute the most important 
cases from the point of view of applications, it is of interest that results can be 


obtained for more general types using only elementary methods. In this paper 
we shall discuss separability of the equations 


(2a) Lily] = 0, (2b) Lily] + ky = 0, 
where 
(3) = + + + (Chr) 


Here A, B and C are functions of x and y, k is a constant, and subscripts refer 
to differentiation. We shall determine necessary and sufficient conditions that 
(2a) and (2b) be separable in coordinates u, v and for the special case A =C, 
B=0 we shall determine the separable systems of coordinates. 

We assume throughout that all functions are twice differentiable in a region 
R of the xy-plane and the corresponding region S of the uv-plane where S is 
defined by the transformation T:x=x(u, v), y=y(u, v). We further assume 
that (2a) and (2b) are elliptic, i.e, AC—B?>0 in R. The equations remain 
elliptic under the transformation T but not necessarily self-adjoint. However it 
is easy to show that the transformed equation becomes self-adjoint when multi- 
plied by the Jacobian J=0(x, y)/0(u, v). 


2. Separability of (2a). We define 
M|y, 6] = + BWeby + + 
In order that (2a) be separable in u, v the transformation T must be such that 
(5) N[u] = F(x, y)U:,  M[u,x]=0, = F(x, 9)Vi, 
L{u] = F(x, y)U2, = F(x, y)V2, 
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where F is an arbitrary function of x, y and where here and in what follows, 
U; are functions of u alone, V; are functions of v alone, i=1, 2, 3, 4. We assume 
that U;, V; are different from zero for all 7. 

The transformation T when applied to (2a) gives, after multiplying by J, 


(6) (APu)u + + (B’u)s + = 0, 


where A’=JN[u], B’=JM[u, v], C’=JN[v]. Hence conditions for separability 
are also 


(7) JIN{u] = M[u,v] =0, JN[v] = Us. 


This latter set of conditions is equivalent to (5) in the sense that (7) implies (5). 
This equivalence can be established without difficulty by making use of the self- 
adjoint property of (6), namely 


Let a, be defined as follows: @= i= We now prove the 
following theorem: 


THEOREM 1. The necessary and sufficient conditions that (2a) be separable in 
u,vare 1) & and 6 satisfy the Beltrami equations 


(8) Ai, + Ba, = Wi,, Ba, + Ca, = — Wi, 
with W=(AC—B?*)"?, and 2) functions U(u), V(v) exist such that 
(9) W = U(u)V(2). 


We first prove the conditions are necessary. From (5) we have N[a]=N[é] 
=F, M[a, 6]=0, and from these the Beltrami equations follow at once. Then 
from the Beltrami equations we get 


N[a] = W(a.8, — 4,0.) = 


with H=1/J. Combining this result with equations (7) gives JN [a] = U;V3UP 
= W(U,V:)-"? and similarly JN [6] = UsVi Vi = Hence, (9) fol- 
lows immediately. 

To show that the conditions are sufficient, we note first that the equations 
N[a]=N[s]=F, M[a, 4] =0, follow from the Beltrami equations. Moreover, 
as above, we have = HU(u) V(v) (Ui Vi)? and a similar 
equation for N[é]. These equations are obviously equivalent to (7) which com- 
pletes the proof of the theorem. 


3. Separability of (2b). If we transform (2b) by means of T, multiply by J 
and make use of (6) and (7), we get 


ky 


* 


1 1 
(10) + (Vibe) + 
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Thus we have the theorem: 


THEOREM 2. The necessary and sufficient conditions that (2b) be separable in 
u, v are that a, satisfy (8), U, V exist which satisfy (9) and N[a] is such that 


(11) 1/N[a] = f(u) + g(u). 

We omit the proof which follows readily from Theorem 1 and proceed to 
discuss the special case A = C, B=0. This is a special case that occurs frequently 
in physical applications, e.g., in the study of vibrations of nonhomogeneous 


membranes or flow of heat in substances whose thermal conductivity is a func- 
tion of position. 


4. Differential equations for f, z, F=log U,G=log V. With A=C, B=0 the 
separability conditions of Theorems 1 and 2 become 


(12) Us = Vy, A(x, y) = UV, J= UV(f(u) + g(v)), 


where we have omitted the bars from u and v. Thus we shall be concerned with 


the transformations T: x=x(a, 6), y=4(#, 3), and this involves no loss in gen- 
erality. 


Let z=x-+7y be an analytic function of w=u-+iv, and write z=¢(w). Two 
transformations Ji, are then defined by 2:=¢:(w), 2.—=¢2(w). We shall say 
that two such transformations are of the same type if ¢: and ¢» are related in 
the following manner: 

1) =az,+5, 

2) m=auw.+B8, «0, a real or purely imaginary. 


We have J=23+2? =|¢'(w)|? and 
2 log | 6’(w)| = F(u) + G(x) + log (f(u) + 
with F(u) =log U(u), G(v) =log V(v). Since log |’ (w)| is a harmonic function 
we find that 
(13) + + + + 2) — = 9, 


where here and in what follows, primes indicate differentiation. We omit the 


details, which amount to straightforward differentiation and integration, and 
find from (13) that 


(14) F” = mf? — m, G" = — mg? — — m, 


(15) f'? + mf — nf? + pf — of +r = 0, 

g'* — mgt — ng* — p*g? — gg —7 = 0, 
where m, n, m, p, g, 7 are real constants. These equations determine, to within 
arbitrary constants, the F, G, f and g. We note that putting f= —g, du =idv, 


changes the first of (15) into the second and vice versa. Hence we need consider 
only one of (15) and we choose the one involving g. 
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It is convenient to define h(g) = m*g*+-ng*+ p*g?+q¢g+r and write the second 
of (15) as 


(16) g? — h(g) = 0. 


5. Classification of solutions of (14) and (15). We shall discuss only those 
solutions which are real in both f and g. Constants of integration will be desig- 
nated by aj, bj, ¢:, u;, v;, fi, gi With i=1, 2, - - - , and are to be real. The function 
m(u?—v?)/2+au+bv+ca will be denoted by R(u, v). We assume m0; the 
extension to m =0 is obvious. We shall not discuss all possible cases that arise 
from the solutions of (14) and (15) but only those which seem of most interest 
and simply state the results for the others. 


Case 1. m=n=0, p?=—m,/4. Define complex constants y and 6 by 


= (a1 —1b;)/2, 5=any complex constant whose real part is ¢,/2. From (14) we 
find 


F(u) = — + au + a2, G(v) = + Biv + be, a2 + = 
(a) p=q=r=0. Here 2 log |b’ (w) | = R(u, v)+log (fitgi) and 
= Gat) exp + at + dat. 


Hence we say that all transformations here are of the type 


= 
(b) p=0, gO. In this case we find 
where C and ¢; are complex constants. The transformations are of the type 


(c) h=p?(g—g:)?. Transformations are of type (a). 


(d) p¥0, h= p?(g—gi)(g—ge), g: and ge real and distinct. All transformations 
here can be shown to be of the type 


ow) = f sin (of + $1) exp 


where a; and { are complex constants. 
It is of interest to notice that the cases where F and G are both constant, 
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which implies m, =0, are the cases in which the separable coordinate systems are 
the four cylindrical coordinate systems found by Eisenhart [2]. These are the 
only separable systems for A(x, y) =const. 


Case II. m=0, n¥0. 
(a) h=n(g—g.)*. Transformations are of type I(b). 
(b) h=n(g—g:)*(g—g2). Transformations are of type I(d). 


(c) h=n(g—g:)(g—g2)(g—gs), g1, Ze, gs real and gi >ge>gs. The solution here 
can be found in terms of elliptic functions.* We have 


g = {gi — gz sn? [M(v + 01), ki] } + 01), hil, 
where kj = (g2—gs)/(gi—gs) and M=4n(gi—gs), and 
f = (g1 — g2) [M(u + ko] — gu 


with R=1—£?}. In the remainder the modulus & is to be associated with v and 
ke with u. Explicit reference to these moduli will be omitted. Also for brevity 
we write v for v-+v, and u for u+m. Thus we have 


f+¢= (ge: — g2)| sn Mw|%(1 — sn? Mv dn® Mu)/cn? Mo, 


where w=u-+iv and the modulus of sn Mw is ko. 


The functions F and G can be expressed in terms of the theta functions, 
thus: 


F = p?u? + au + a2 + 2 log 0(Mu, ks), 
G = — po? + + b2 + 2 log kz), 


where p?=m/2+(n/4) [gs+(g:—gs)E/K], and K and E are complete elliptic 
integrals. 
Finally we find that 


2 log | ¢’(w) | = — 0?) 
+ log [(g: — g2)- | sn (Mw) -O(Mw) |20(0)?], 


and hence the transformations for this case are of the type 
= + £9) sn (out + £0) exp 


(d) h=n(g—g:)[(g—a)?+6?], a and b real constants. This case can also be 
expressed in terms of the elliptic functions. We omit the details and give only 
the result that the transformations are of the type 


=f + — en Gout + exp ar. 


* The necessary details conccrning elliptic functions can be found in [5]. 
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Case III. m¥0. 

(a) h=m?(g—g,)*. Transformations are of type I(a). 

(b) h=m?(g—g:)*(g—ge2). Transformations are of type I(b). 

(c) h=m?(g—g:)*(g—ge)?. Transformations are of type I(a). 

(d) h=m*(g—g:)?[(g—a)?+b?]. Transformations are of type I(a). 
(e) h=m?(g—g:)(g—ge)(g—gs)*. Transformations are of type I(d). 


(f) h=m?(g—g:)(g—g2)(g—gs)(g—gs). This case can be shown to yield trans- 
formations that are of type II(c). 


(g) h=m*(g—gi)(g—gz) [(g—a@)?+6?]. This case can be shown to yield trans- 
formations that are of type II(d). 


The results of this section can be summarized as follows: The transforma- 
tions which separate (2b) with A =C, B=0 and which are defined by the real 
and imaginary parts of z=¢(w) with z=x+iy and w=u-+72, are of the type 


(17) ow) = f exp 


where x(¢) assumes one of the values, 1) 1, 2) (+, 3) sin(ig+f), 4) O(ag+h) 
5) {@(ois +1) [1+en(oig ]} 

The coefficient A = A (x, y) = U(u(x, y), v(x, y)) are those for which separable 
transformations exist. Although all transformations which are separable are 
given by (17) it is obvious that an explicit expression for A cannot be written 
in most cases. 
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AN OPTIMAL SEARCH PROCEDURE* 
SETH ZIMMERMAN, Dartmouth College 


The problem which we are concerned with is the following. We are givena 
set of objects, one of which has been previously selected without our knowledge, 
as the one which we must determine. In addition, we have assigned to each ob- 
ject a probability of having been selected, the sum of all such probabilities being 
unity. We must determine the selected object by a series of yes-no questions, 
arbitrarily dividing the set of objects into two groups and asking which it is 
contained in. We continue this yes-no process until the selected object is known. 
For example, we might use the yes-no process to determine which letter of the 
alphabet was randomly selected from a book, given the probabilities (the fre- 
quencies) of each letter. We would now like to know what the best strategy of 
division is, where by a strategy we mean a plan which indicates exactly what 
the next division will be, given that the selected object is contained in any par- 
ticular group of the original objects. For example, given four objects with as- 
signed probabilities (.4), (.3), (.2), (.1), our strategy might be, first divide in 
the following way: (.4, .1), (.3, .2). Then, if the object is learned to be in the left- 
hand group, further divide, (.4), (.1), and if it is in the right-hand group, divide, 
(.3), (.2). The expected value, that is the probable number of divisions for this 
strategy, would be two. 

The questions which will now be answered are: Given any set of objects, 
what strategy should be used to obtain the smallest expected value (the optimal 
strategy)? What is the smallest expected value? 

It is important to observe that although this seems like a method of repeated 
division, it can be viewed rather as one of combination. We could start with the 
entire set as in the above strategy, and describe where the first division, and 
subsequent divisions dependent upon the outcome of the preceding ones, are 
going to be made. Instead, we can proceed in precisely the reverse manner, be- 
ginning with the several isolated objects and indicating a series of repeated 
combinations, the first few of which correspond to the last few divisions, and the 
final one to the first division. Since in our strategy divisions are always into two 
groups, each of the combinations will similarly involve two groups. 

In the above example then, the method of combination corresponding to 
the strategy used would be one of the following two. First, combine the two 
objects with the probabilities (.4) and (.1); next combine those with probabilities 
(.3) and (.2). Finally combine the two groups corresponding to (.4, .1) and 
(.3, .2). The second method of combination would simply contain a reversal of 
the first two steps of the preceding one. It is seen from this example that a 
particular strategy need not have a unique method of combination. However, 
it should be fairly evident that any method of combination will always have 
only one division strategy corresponding to it. Therefore, in order to determine 


* This research was supported by the National Science Foundation through a grant to the 
Dartmouth Mathematics Project. 
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the optimal strategy of division of a set of objects, it is only necessary to deter- 
mine a corresponding method of combination. 

A concept of importance is that of the m-set. If we are given a set of n 
objects, then we call the array of the corresponding probabilities, arranged from 
left to right in descending order of magnitude, an n-set. This left to right ar- 
rangement has no significance except that it will be convenient later on. When 
we combine two of the objects of the m-set, and assign the sum of the two 
probabilities as the probability of the new “object”, we have essentially reduced 
the n-set to an (n—1)-set. Any method of combination is then a way of making 
a combination of two objects n—1 times until the 1-set is reached. 

Another concept which we will use is that of the partial sum. At any stage 
of combination, the sum of probabilities of all the pairs of combined objects we 
call the partial sum. The probability of one of the original objects can be in- 
cluded more than once in a partial sum, since after it has been combined to form 
a new object, this latter object may again be combined. 

We can now state the main theorem, although it will not be proved until 
later, following a series of lemmas. 


THEOREM. The optimal rule of combination is: At each stage combine the ob- 
jects with the two smallest probabilities. 


Lemma 1. The method of combination of any n-set is, after the first step, pre- 
cisely the same as that of some (n—1)-set. 


Proof. In the method of combination of the n-set, there must be a first com- 
bination. Therefore, the (7—1)-set equivalent to the m-set, except for the re- 
placement of these two objects by the combined one, has, from then on, the 
same method of combination as the m-set. We have, in fact, referred to the n-set 
as being reduced to the (n—1)-set by the combination of the two objects. It 
follows that the expected value of the m-set should be that of the (m—1)-set 
plus the probability that we will have to ask another question if the object is 
contained in the two combined ones. Since the latter is exactly the sum of the 
two original probabilities, we can write the equations 


E, = E,-1+ (Pm + Pn), Ey = 0, 


where E, is the expected value of a given r-set, and », and , are the probabili- 
ties of the combined objects. If we proceed by induction, the result is that E, 
is the sum of »—1 combinations of probabilities and the problem is really how 
to combine probabilities two at a time so that the total of these probabilities 


is a minimum. Each original probability is going to be combined one or more 
times. 


LemMA 2. In an optimal method of combination, a larger probability cannot be 
combined more times than a smaller one. 


Proof. lf we had a method whereby a larger probability were combined more 
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often than a smaller one, we could interchange the two probabilities in our oper- 
ations and get a smaller expected value. 


Lema 3. We are given the optimal value of an n-set, and any two probabilities 
pr and p, where p,> p,. If we increase p, by any amount k, and decrease p, by the 
same amount, the optimal value of the transformed set cannot be greater than that 
of the original. 


Proof. By Lemma 2, p, cannot be combined a greater number of times than 
p. and so if we use the same method of combination for the transformed set we 
can add plus & no more than minus &. It is impossible to get a greater expected 
value in this way. In addition, this may not even be the optimal method for the 
transformed set. 


We are now in a position to give an inductive proof of the main theorem. 


Proof of the main theorem. The rule obviously holds for all the 1-sets, of 
which there is only one. We shall now assume that it holds for every r-set, where 
r assumes the values 1, - - - , »—1. If it can be shown that the rule is true for 
any -set, then the theorem is proved, by induction. 

We are given an arbitrary n-set, arranged in descending order of size, 
Pn, * pr. There are three possible first combinations: 


(1) A combination of objectsr and s, that is, an addition of p, and p,, where 
neither r nor s is equal to 1 or 2. 


(2) A combination of r and s where one of r and s is equal to either 1 or 2. 


(3) A combination of 1 and 2. 
It must be shown that the third method is at least as good as the other two. 


(1) Combining r and s gives a partial sum of (p,+,) and the (m—1)-set 
The next combination must be 1 and 2 according to 
the inductive assumption, giving a partial sum of (,+,+ 1+ 2) and the 
(n —2)-set (Pr+Ps), (pith), 

If, instead we first combine 1 and 2, then we can next combine r and s, 
which may not even be optimal. We arrive at the same partial sum and the 
same (n—2)-set. This method is at least as good as combining r and s in the n- 
set. 

(2a) If 1 and r (r>3) are combined, our partial sum is (f:+),) and we have 
the (n—1)-set The next combination must be 
(p2+ ps); this gives a partial sum of ),) and an (n—2)-set pa, 
(pit?r), (p2+Ps), where if r 4. 

If, instead, 1 and 2 are combined, we have a partial sum of (f:+ 2) and an 
(n—1)-set pn, +--+, (Pitpe), +++, ps. The next combination cen be 3 and r, 
which may not even be optimal. This gives the same total partia| sum as before 
and an (u—2)-set, Pa, , ps By Lemma 3, the 
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new (m—2)-set has an expected value no greater than the original, and since 
the partial sums are equal, this method must be at least as good as the former 
combination. 


(2b) If 1 and r are combined, where r is 3, then one of two situations can 
arise. If ($1+3) >, the proof is analogous to (2a). If however, (f:1+p3s) < pu, 
then, after we have combined 1 and 3 we would combine 2 and the combination 
of 1 and 3. In the second method, after we had combined 1 and 2, we would com- 
bine 3 with the combination of 1 and 2. The resulting (m—2)-sets would be the 
same, and the partial sum of the 1 and 2 combination would be lower than that 
of the 1 and 3 combination. Therefore the latter would be at least as good as the 
former method. 

In the proof of (2) it is possible to interchange the values for p; and 2 so 
as to prove the theorem for the combination of 2 and r. All possible first com- 
binations in the n-set have been treated and hence the theorem is proved. 

A rather interesting situation is that of nonuniqueness. If, at any stage of 
the method, there exist two or more probabilities which are equal, we can ar- 
range them in any order we like, as long as they are adjacent. Further, if these 
probabilities represent not individual objects but groups, then the order will 
affect the final combination procedure. The result will be two or more arrange- 
ments which are equally acceptable. 

As final remarks, let me indicate the outcome of the example presented at 
the outset, that of the determination of a letter of the alphabet. We are given 
a set of probabilities corresponding to each letter and the blank space. They 
range in size from .001 for the letter Q to .200 for the blank. Using the rule of 
combination, the expected number of divisions or questions is 4.154. 

The author would like to thank Professors John G. Kemeny and J. Laurie 
Snell for their assistance in the solution and preparation of the problem. 


SOLUTION OF A SET OF GAMES 
BENJAMIN L. SCHWARTZ, Battelle Memorial Institute, Columbus, Ohio 
(present affiliation: Technical Operations, Inc.) 


1. Introduction. This paper deals with a family of zero-sum two-person 
games in the Von Neumann sense [1]. The rules of the games are simple, but 
their implications are sufficiently obscure that the methods of game-theoretic 
analysis are necessary to determine the solutions. In this respect, the games fall 
into a category which, at present, has too few representatives known; namely, 
those games which (unlike tic-tac-toe) are not so simple that a theoretical analy- 
sis is superfluous, and (unlike chess) not so complicated that the analysis is in- 
adequate to provide the solution. 
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2. Description of the games. The games considered are known by the desig- 
nation of the “n-coin game,” where 7 is a positive integer. The particular case 
of the three-coin game will be used to illustrate the rules of play; the general case 
is completely analogous. 

In the three-coin game, each player A and B is supplied initially with three 
identical coins or other counters. 


First Move. Each player chooses, unknown to the other, an integer between 0 
to 3 (inclusive) and places that number of his coins to his right hand, which is 
thereafter kept closed. The number of coins selected by A and B, respectively, 
in this move will be denoted A: and B,. 

Second Move. Player A attempts to guess the total number of coins now held 
in both player’s right hands. He announces his guess, As, an integer between 0 
and 6. 

Last Move. Player B, having heard A’s second move, now attempts himself to 
guess the total. His guess, B2, must dfffer from Az. 


The hands are now opened and the correct total ascertained. If either player 
has guessed correctly, he wins one arbitrary unit from the other. Otherwise, the 
play is a draw. In practice, when a draw occurs, it is customary to repeat the 
game with A and B exchanging roles, until one or the other wins. However, this 
convention need not affect our theoretical studies, which are concerned simply 
with optimal strategies for A and B in a single play of the game. 


3. Heuristic discussion. Each player has some benefits under the rules, and 
it is not immediately obvious which has the greater. Player A has the choice of 
seven numbers (0, 1, 2, 3, 4, 5, or 6) for his second move, while B has only six. 
Moreover, if A guesses the sum correctly, then B’s last move is completely 
futile, since he is defeated already. However, B has the advantage of hearing 
A; before naming Bz, and Az may give B some clue as to Ai, A’s choice at the 
first move. For example, if A2=6, then B can be reasonably sure that A’s right 
hand holds three coins. This information will be useless in the case when B is 
also holding 3, but in any other case, B can be assured of a win. 

Note that A can cancel this advantage to B by making a guess of A2=3 at 
his second move. This reveals nothing to B about A’s first move. On the other 
hand, in order to utilize this method of annulling his opponent’s advantage, 
A must sacrifice his own advantage of complete freedom of choice at Move 2. 

The relative importance of each of these features of the game will be re- 
vealed in the theoretical analysis below. 


4. A simple special case. Before undertaking the analysis of the general 
n-coin game, we shall consider some special cases. Our first example is the one- 
coin game. Observe first, that among the possible strategies for either player, 
there are some which offer no possibility of winning. For example, for A to play 
A,=1, A2:=0 is manifestly absurd. It is possible to exclude such strategies by 
employment of the dominance principle [1], but it is hardly necessary to use 
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so sophisticated a technique to effect so elementary a result. In the remainder 
of this paper, only strategies which give the player a possibility of success 
(“feasible” strategies), will be considered. This does not destroy the generality 
of the results, since, as remarked above, all others may be excluded by dom- 
inance. Henceforth, we shall use the term strategy to denote feasible strategy, 
unless otherwise specified. With this convention, A has exactly four strategies: 


A-I : A, =0, A2=0, A-III: Ai =1, A. =1, 
A-II : A, =0, Az = 1, A-IV: A, = 1, Az = 2; 


and B has only two strategies 


x 0 

B-I: B,=0, as A2=}1, 
1 2 

0 

B-II: By} = 1, as 
x 2 


In case B,=0 and A:=0, B presumably has a lost game, and his second move 
is immaterial; a similar comment holds for the case B} =1, A2=2. These cases 
are denoted above by the (purely arbitrary) notation Bi=x. As a matter of 
practical expediency, B may be best advised to take B.=1 in these two cases 
since this will give him an opportunity to win in the unlikely event that A is 
playing one of the nonfeasible strategies. 

The payoff matrix for the game is given in Table 1. 


TABLE 1. PAYorr FOR THE ONE-Corn GAME (PAYOFF TO A) 


B-I 

A-I 1 —1 
- —1 1 
A-Ill 1 —1 


A-IV -1 1 


From this table, it can be seen that basic optimal strategy mixes [3] for A 
include (and are limited to) 


I and II in equal proportions, 
I and IV in equal proportions, 
II and III in equal proportions, 
III and IV in equal proportions; 


and B has only one optimal strategy, which is to mix B-I and B-II in equal 
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proportions. The value of the game is zero; no play will end in a draw, and 
neither player can protect himself against loss if he attempts to take advantage 
of an opponent’s error. In several of these respects, we shall see that the one- 
coin game is atypical. 


5. The two-coin game. When n=2, a game is obtained which is typical in 
all important respects of all higher values of m, and still has few enough pos- 
sibilities that detailed examination is practical. 

It will be convenient to extend our notation as follows: Let B, be player A’s 
guess of B’s first move, Bi, i.e., B, =A2—A1. Then A’s strategies in the two-coin 
game are nine in number and may be tabulated as in Table 2. 


TABLE 2. STRATEGIES FOR PLAYER A IN THE Two-Corn GAME 


Strategy number B, Ai A: 
A-I 0 0 0 
A-Il 0 1 1 
A-IIl 0 2 2 
A-IV 1 0 1 
A-V 1 1 2 
A-VI 1 2 3 
A-VII 2 0 2 
A-VIII 2 1 3 
A-IX 2 2 4 


The strategies for Player B may be described by the following convention. 
Any strategy for B determines a column in the payoff matrix, provided the A- 
strategies are listed down the side as in Table 2. Conversely, any column of 
nine entries of +1, —1, or 0 determines a B-strategy provided certain conditions 
are observed. The restrictions are stated below for the two-coin game. The form 
these take in the general case (with parameter 7) is given in brackets. 


Condition 1. Let the A-strategies be classed into the following three [n+1] 
mutually exclusive sets: 


All strategies for which either B, =0, or B,=1, - - - , or Bi =n. Then there 
must be +1 entries in the column defining the B-strategy against all A-strategies 
in exactly one of these sets, and no other +1 entries may appear. Since any B- 
strategy whatsoever must have a value of B;, one of the three [7 +1] estimates 
B, must be correct, i.e., lead to a payoff to A of +1. We shall denote a B- 
strategy in which B, =7 as a strategy of Type 7. 


Condition 2. The remaining 6 [2(n+1) ] entries in the column may be filled 
with 0 and —1 entries, but one and only one —1 entry* may be inserted against 


* If nonfeasible B-strategies are considered, this condition should be modified to read “ . . . not 
more than one —1 entry... .” 
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any single value of Az. Thus, strategy A-VI and A-VIII both display A: =3. 
In the column defining the B-strategy (say of Type 0), player B must select 
his call Bz corresponding to A;=3. He may select it so that he will win against 
A-VI, 1.e., Bz =2, or against A-VIII, i.e., Bg =1, but he cannot choose a strategy 
1 which will win against both. 

. Brief consideration will establish that any column conforming to conditions 
(1) and (2) above will define a B-strategy, and conversely. It is not difficult to 
S determine all B-strategies in the two-coin game, of which there are exactly ten. 
L These strategies and the payoff matrix for the game appear in Table 3. (In this 
table +1 and —1 are abbreviated to + and — respectively, for brevity.) The 
reader is urged to construct Table 3 for himself. By so doing, he will greatly 
facilitate his understanding of the proofs given in subsequent sections. 


TABLE 3. Payorr MATRIX, Two-Coin GAME (PAYOFF TO A) 


A-Strategies B-Strategies 
Stratesy 4, B-I B-II B-III B-IV B-V B-VI B-VII B-VIII B-IX B-X 
— Type 0 Type 1 Type 2 

A-I + + + + - 

A-II + + + + - - - 0 

A-III 0 2 + + + + - 0 - - 0 0 

A-IV - - - + + 0 o - 

A-V - - 0 + + 0 

A-VI ins o - + + 

A-VII 2 2 o - - + + + + 
A-IX 2 4 - + + + + 


A somewhat tedious but straightforward examination of this payoff matrix 
suffices to establish the following results. 


Proposition 1. The game is fair. [2] 


PROPOSITION 2. Player B has several optimal mixed strategies. Every one of the 
ten B-strategies is active in at least one of the optimal strategy mixes available to B. 


PRoposITION 3. Player A has only one optimal mixed strategy; namely, A-ITI, 
| A-V, and A-VII in equal proportions. 


6. The general theorems. Statements analogous to Propositions 1, 2, and 3 
hold in the general case. We shall present proofs of these general theorems. Be- 
cause of the awkward nature of the operations, it will prove convenient to ~ 
present the proofs in the language of a specific game, and we shall use » =3, the 
three-coin game, for this purpose. However, it will be observed that the methods 
we employ may immediately be applied to any value of n. 
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THEOREM 1. Jf n>1, the n-coin game 1s fair. 


Proof for the three-coin game. Using notation like that introduced in the pre- 
ceding section, we first show a mixed strategy for B in which the sum of the 
entries in any row is nonpositive. This implies that if the strategies so selected 
are played in equal proportions against any A-strategy, pure or mixed, the ex- 
pected payoff to A will not exceed zero [2]. 

Select any nonnegative integers, 1#j (7, 7<3). We now choose B-strategies 
of Type zand j. In the column defining the Type 7 strategy, we insert —1 against 
every A-strategy which has Bi=j; similarly, the Type j strategy has —1’s 
inserted in those positions where B,=i. The other —1 entries are immaterial 
as long as Condition 2 is observed. It is obvious that such strategies satisfy the 
Conditions 1 and 2 of the preceding section and assure that the row sums are 
nonpositive. B can now play these two pure strategies with equal weight in a 
mixed strategy. This completes the proof that B can force an expected payoff of 
0 or less. 

Player A can force an expected payoff of 0 or more by playing an equal 
weight mix of the four stragegies: 


B, = 0, Ar = 3; B, = 1, As = 3; B, = 2, Az = 3; B, = 3, As = 3. 


By Condition 1, every B-strategy must lose against exactly one of these 
four. Also, according to Condition 2, each B-strategy will win against one and 
only one of them. Therefore, any column sum is zero.* This completes the proof. 
The restriction m>1 is, of course, unnecessary since we have already proved the 
one-coin game to be fair. However, the method of proof used here requires m> 1. 


THEOREM 2. In the n-coin game, any pure B-strategy is active in some mixed 
optimal strategy for B. 


Proof for the three-coin game. Let S be a given pure B-strategy. We shall 
produce a four-strategy mix of B-strategies including S which has all row sums 
nonpositive. 

Strategy S is of some type, say 7; rename S=S;, and begin construction of 
the other three strategies in the mix by choosing them to be one of each of the 
three types among Types 0, 1, 2, and 3 obtained by omitting Type 7. Name each 
strategy by a subscript defining its type, 7.e., So, Si, S2, and S3, of which one is 
S and is completely defined, and the others still have —1 entries to be inserted 
in the columns of the payoff matrix defining them. 

We shall insert these —1 entries such that each row has at least one such —1 
entry. Since each row (defined by B, and A:) has exactly one +1 entry (by con- 
struction) the four strategies So, S:, S2, and Ss; may be played with equal weight, 
giving nonpositive payoff to A. This will prove the theorem. We will assign these 
—1 entries systematically. For each value of Az, there are one or more A- 


* If nonfeasible B-strategies are considered, any column sum is nonnegative. 
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strategies, 7.e., one or more rows in the matrix. Table 4 shows the number of 
occurrences of each A: value. 


TABLE 4, OCCURRENCES OF VALUES OF A? IN FEASIBLE 
A-STRATEGIES IN THE THREE-COIN GAME 


Values of As 0 1 2 3 4 
Number of occurrences 1 2 3 4 3 


Now 


We begin by assigning the four —1’s which fall on rows where A2=3. One of 
these already is entered, in the S; column (since it is a feasible strategy). This 
lies on a row for which B, =j (say), ji. In column S; insert a —1 against an 
A-strategy for which A2=3, and B, +i, j, say at B,=k. This is always possible 
since for each B,, an A-strategy exists with A2=3. Next, in column S;, insert 
—1 against an A-strategy for which A2=3, B, Xi, j. k, say at Bi=m. Lastly, in 
column S», insert a —1 against strategy A.=3, Bi =i. It is obvious that this 
construction is always possible and complies with Conditions 1 and 2. 

In similar manner, three —1 entries are made against A-strategies with 
A:=4. One of these is already determined in column S;. Suppose it occurs 
against the A-strategy with Bi=j’. In column Sy, insert a —1 against an 
A-strategy with A.=4, B, +i, say at B,; =k’. This is always possible. In column 
Sy, insert a —1 against an A-strategy with A,.=4. This may or may not occur 
where B,=i, but it, too, is always possible, since three distinct values of 
B, (Bi =1, 2, and 3) permit A-strategies with A,=4. 

Similarly, the three —1 entries for Az=2, the two entries for A2=1, and the 
two entries for A2=5 can be assigned. The one entry for A2=0 is inserted in any 
one (or more) of columns Si, S2, and S3; the entry for Az =6 is inserted in any one 
(or more) of columns So, Si, and S:. This completes the construction. It can 
readily be seen that the required conditions are satisfied and the theorem proved. 


THEOREM 3. Jf n>1, in the n-coin game, no optimal strategy for Player A has 
any active strategies except: 


Proof for the three-coin game. These strategies have been shown to be active 
in one optimal strategy mix. We now prove that no other pure strategy can be 
active in an optimal mixed strategy for A in the three-coin game. 

Let T be an A-strategy with A2+*3. First consider the case A2<2. We choose 
the following two-strategy mix for B: (1) S2, a Type 2 strategy with —1 entries 
against T and every A-strategy for which B,=3, (2) S:3, a Type 3 strategy with 
—1 entries against every A-strategy for which B;=2. Other —1 entries are 
immaterial as long as condition 2 is satisfied. If these strategies are played with 
equal weights, a nonpositive expected payoff to A is assured for every pure 
A-strategy, and a negative payoff of —1/2 is obtained when strategy T is 
played in particular. Similarly, if A2>4, a mix of a Type 0 and Type 1 strategy 
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suffices to provide a negative expected payoff. 

It remains to consider the case A2=2 or 4 (in general, Aa=u+1). Suppose 
A:;=2. For this case, we construct a four strategy optimal mix for B in a manner 
similar to that employed in the proof of Theorem 2. We begin by defining four 
strategies of Types 0, 1, 2, and 3, denoted So, S:, S2, and S3, respectively. S3 will 
be called the key strategy. Suppose strategy T is defined by B,=i, A2=2. 
Choose an integer 7 $3, 7#i. In column S,, insert a —1 against strategy T. This 
is possible since Bi(T) =i, and 7#7. Now in column Sj, find an A-strategy with 
A:=2, Bi Xi, j. Such a strategy exists, say Bi=k, since three occurrences of 
A:=2 are shown in Table 4. Furthermore, k#3, since the A-strategy Bi =3, 
A, =2 is not feasible. In column S;, assign a —1 against the A-strategy with 
A:=2, B,=j. Also, in column S;, assign a —1 against the A-strategy with 
Az B, = k. 

Now enter a —1 in column 5; (the key strategy) against T. This is the fourth 
—1 entry against A-strategies having A,:=2, and therefore, the second such 
entry against 7. The remaining entries in column S; are entered in any manner 
to create a feasible strategy and satisfy Condition 2. The remaining entries in 
columns So, S:, and S: are filled in the manner described in the proof of Theorem 
2 (except for the —1’s on lines where A2=2 which have already been defined). 
The resulting strategy mix has one +1 and one or more —1 on every line, and 
in particular on line T, it has two —1’s. Hence, when these strategies are played 
in equal proportion against any A pure strategy, the payoff to A is nonpositive, 
and against T in particular, it is negative and equal to —1/4. 

This proves the theorem for the case A2=2. For A2=4, the proof is similar, 
except that So is the key strategy. This completes the proof of the theorem. 


Remark. It is clear from Theorem 3 and the proof of Theorem 1 that A must 
play these strategies with equal weights. 


7. Concluding remarks. Despite the content of Theorem 1, it seems apparent 
that the advantage in the game lies with Player B wherever n>1. He has a wide 
selection of optimal strategy mixes, and in view of Theorems 2 and 3, he can 
select them in such a way that he can penalize A for any departure from the 
strategies in A’s single optimal mix, while simultaneously protecting himself 
against any loss. Player A, on the other hand, is limited to his single mixed 
strategy, and may not deviate from it without risking loss. Unfortunately for A, 
while his strategy assures a nonnegative expected payoff, it also assures a non- 
positive payoff, even if B blunders. While a quantitative measure of the value 
of such a property of a game is not generally agreed on [3], it is difficult to 
imagine any interpretation which would not make this circumstance a plus value 
for B. 

It has been pointed out to the author that in view of this, a wiley A-player 
may choose a nonfeasible strategy such as A1 = 3, A2=2 in order to force a drawn 
play. In accordance with the convention mentioned in Section 2, this player 
would then attain on the next play the B-role, which is preferable. The inter- 
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pretation of this strategem will be left to the reader. 

It will be noted that in the three-coin game, one half of the plays end in 
draws (provided both players play perfectly). This may be seen by examination 
of the proof of Theorems 1 and 3. The single feasible strategy mix which A will 
employ will win one play out of four against any B-strategy, and lose one out 
of four against any feasible B-strategy, and draw in the remaining two plays. 
In general, the ratio of draws to total plays is readily seen to be (n—1)/(m+1). 

An interesting consequence of the theory is that B may play the game with 
only one coin for any m. In the proof of Theorem 1, a construction is given 
which provides an optimal mixed strategy for B comprising Type 0 and 1 
strategies only. The content of Theorems 2 and 3 must be modified, of course, 
if B has only one coin. However, the methods employed in the proof of the 
theorems may readily be applied to determine the new relationships. 
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MATHEMATICAL NOTES 
By Roy Dvusiscu, Fresno State College 


Material for this department should be sent to Roy Dubisch, Department of Mathematics, 
University of California, Berkeley 4, California. 


NOTE ON STIRLING’S FORMULA 
T. S. Nanjunp1an, Mysore University, India 
Stirling’s asymptotic formula, namely, 
~~ +/(2xn) -(n/e)*, 
is usually proved by showing that 
(1) n\| = +/(2nn) -(n/e)*-e™, 0< yn <1/(12m), m=1,2,---. 


Herbert Robbins (this MonTHLY, vol. 62, 1955, pp. 26-29) has shown by an 
elementary method that the estimate for y, in (1) can be replaced by the im- 
proved estimate 


(2) 1/(12m + 1) < yn < 1/(12n), m=1,2,---. 


In this note, we shall prove, by a simple refinement of the argument of Robbins, 
the stronger and more precise result 
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1 On 
3, = 2 —— 0< 4, =1,2,---, 
( = +/(2xn) -(n/e)"- exp (— < n 


as envisaged by Euler’s summation formula. 
For the proof, let a, =n!n—/2(e/n)", 1, =log dn, m=1, 2, We have 


an 1 1 1 
-—(1+-) = (n+ =) log (14+ —) 1, 
e n 2 n 


Using the familiar series 


1 1 1 
og (1+ =) + ) 


we have 


1 1 
3(2n + 5(2n + 


First we derive, as usual, from (*): 


<l n+1 < + (2n + 12n(n + 1) 


On rewriting this in the form 0 </, —In41 < (1/12) { (1/m) —(1/(n+1)) }, we get on 
the one hand 


(4) In > n=1,2.--- 

and, on the other, 

12n 12(n + 1) 


From (4) and (5) we obtain the known result that the sequence of intervals 
Ini ly —(1/(12n)) <A<1,, n=1, 2, - ++, whose lengths evidently tend to zero, 
forms a nest assuring existence of the limit 


(6) = lim 


This, in turn, assures the existence of the limit 


(7) a=lima,=eée> 0. 


We use the known value 
(8) a = (27)1/?, 
Indeed, (8) follows, after (7), on letting n—>© in the easily verified relation 


| 
| 
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Gon 1-3---+-(2n —1) 
and using Wallis’ well-known product for 7. 

We now proceed to improve (4), so as to lead to (3), by finding a suitable 
series majorized by (*). The series obtained by replacing the pth term in (*) by 
3-?(2n+1)-*? for each p was employed by Robbins and led him to (2). It fails, 
however, to yield the more precise estimate we have in mind. On the other hand, 
if we retain the second term there and allow the others to be replaced as above, 
we still get a series majorized by (*) and it is interesting to find that it just 
yields our result. In fact, we have 


3(2n + 1)2 3(2n + 1)?  32/(2n+1)4 


12 n+ 1) 6n(n + 1) 180n?(n + 1)? 4n(n + 1) 
1 1 1 1 
12n(m + 1) 6n(n + 1) 180n?(n + 1)? 2n(n + 1) 
1 3n(m + 1) +1 


12n(n +1)  360n%(n + 1)? 
As before this gives 


1/1 1 1 /1 1 
which shows that 
1 1 1 1 


It follows from (5) and (9) that the sequence of intervals 
12n 12 360n* 


whose lengths obviously tend to zero, forms a nest assuring the existence of the 
limit (6), since the endpoints of J, differ from those of J, by o(1) as n>. 
The result just proved yields what we desire. In fact, we may now write 


(10) 


<1, m=1,2,--- 


Recalling the definitions of a, and J,, (10) together with (7) and (8) clearly 
establishes (3) and the proof is complete. 


r 
1 
| 
’ 
1 On 
n 360n* 
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ESTIMATION IN A CERTAIN PROBABILITY PROBLEM 


GERALD A. HEveEr, Concordia College and Remington Rand Univac 

An example sometimes quoted as a “nonintuitive” probability is the follow- 
ing: Given a collection of K people, find the probability that some two of them 
were born on the same day of the year (assuming that no one was born on 
February 29, and that people are born with equal probability on other days). 
If K is 23, the probability slightly exceeds 1/2, and if K is 50, the odds are 
about 33 to 1 in favor of such a pair. 

This is of course a special case of the following problem. Let an experiment 
with N equally likely outcomes be performed K times. What is the probability, 
P(N, K), that at least one of the outcomes occurs twice? P(N, K) =1—Q(N, K) 
where Q(N, K)=N(N-1) is the 
probability that all K outcomes are distinct (K <<). For given N and K this 
may be evaluated directly, or approximated with the aid of Stirling’s formula. 

The problem becomes more difficult, however, if one assigns values to 
P(N, K) and N, and attempts to solve for K. We are then confronted with an 
equation 


N} 
— = |} 
(N — 
or, if Stirling’s formula is used, 


( N ) N-K+1/2g-K = 1 — 
N-K 


to be solved for K, where t= P(N, K). Such transcendental equations can be 
solved approximately, of course, when WN and # are specified, but to exhibit 
K=f(N, t) explicitly seems rather difficult. It is the purpose of this note to show 
that for 0 <#<1, K is given asymptotically by L(t)/N, where L(é) 
=4/—2 log (1—#). Furthermore, except for extreme values of ¢, this approxima- 
tion is very good even for small N. 

We first prove the 


Lema. Q(N, K) 2a, where a is constant, 0<a<1, then K/N-0as No 


Proof.* Since 1—x <e~* for all x, we have 


K-1 


so that K(K—1) $2N log (1/a), which implies the conclusion of the lemma. 


* The author is indebted to the referee for suggesting this simpler proof of the lemma. 


4 
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THEOREM. Let K, N, and t be related by the equation 
N! 
1 — 
(N — K)!N*® 
(We may assume all three are continuous real variables, O<K<N.) Then 
K? 
—2N log (1 — 2) 
Thus K is given asymptotically by K=L(t)/N, where L(t) =./—2 log (1—#). 


Proof. By Stirling’s inequality we have 
( N — K)- *) 
N-K. 12(N — K) 


N! N N-K+1/2 12N 
(N — K)!N* N-K (12N — 1) 
Since the middle member above is 1—?#, we have, by inverting, 
(1 — = 1/[((1 — (1+ 6], 


where e—0 as N+ (since when N-~), 
Taking logarithms gives 


—log [(1 — )(1+ 0] = K+ (N — K + 4) log (1 — K/N) 


r 


= K — (N — K) — —+ — log (1 — K/N) 
r=1 i N° 


Ke @ 4 
— — log (1 — K/N 


Km +3 (2 1 


The lemma shows that K/N-0 as N-, so we have —log [(1—¢)(1+6)] 
=e'+[K?/(2N)](1+e’), where e’ =} log (1—K/N)-0 as No and 
2 Ky Ks 
0 


as Thus 
K? —log [(1— — 


2N i+é 
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so that 


K? —log [((1 — + — 
—(1 + log (1 — 4) 
It is of some interest to see how good the approximation is for several values 
of N. The author has done this with the aid of the Univac 1103A and some of 


the results are tabulated below. In this table, K is the least integer such that 
N!/|(N—K)!N*]<1-t and K;, is the integer nearest L(t)-/N. 


t=% t= 3 t=} t = 31/32 

N K Ki K Ki K Ki K Ki K Ki 

4 2 1 2 2 3 2 4 3 5 5 
4? 2 1 4 3 5 5 7 7 10 il 
4° 3 2 7 6 10 9 14 13 21 21 
4‘ 5 4 13 12 20 19 27 27 42 42 
45 9 8 25 24 38 38 54 53 84 84 
48 17 16 49 49 76 75 107 107 168 168 
47 33 32 98 97 151 151 214 213 337 337 
48 65 65 195 194 302 301 427 426 674 674 
4° 130 129 389 388 604 603 853 853 1348 1348 
4to 259 258 778 777 1206 1206 1706 1705 2696 2696 
4u 517 516} 1554 1553 2412 2411 3411 3410 5392 5392 
42 1033. 1032} 3108 3107 4823 4823 6821 6820 | 10784 10784 
43 2065 2064] 6215 6214 9646 9645 | 13641 13641 | 21567 21568 
44 4130 4129 | 12429 12428} 19291 19291 | 27282 27281 | 43135 43135 
4% 8258 8257 | 24856 24855 | 38582 38581 | 54563 54562 | 86270 86271 
46 | 16515 16514 | 49712 49711 | 77164 77163 | 109125 109125 | 172541 172541 
417 | 33029 33028 | 99423 99422 | 154326 154325 | 218250 218249 | 345082 345082 


AN ANALYTICAL EXPRESSION FOR [X] 
ArtTHuR Porcgs, Los Angeles City College 
THEOREM. Given the functions 


_ Arcsin | sin | 


F(X) 


and G(X) = lim {1+ | sin x(X — F(X))| }%, 


then for all finite X, |X | (the greatest integer not greater than X) may be written as 
[x] = — | F(X) 1]. 
To prove the theorem it is helpful to establish two lemmas. 


Lemma 1. For K an integer and F(K+5) for }<b<1, F(K+0) 
=1-5. 


LemMaA 2. For K an integer and <3, G(K +6) =1; for }<b<1, G(K+0) 


= 0, 
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To prove Lemma 1, compute F(K+6), where 0S6S}. Then 


Arcsin | sin (Kx + bz) | Arcsin | sin (57) | 


Now consider F(K +5), where }<b<1. We have 


F(K + 6) = 


b. 


Arcsin| sin (Kx + bn) | Arcsin | sin (57) | 


F(K + 6) = 


a — br 
= =1-d). 


The proof of Lemma 2 is also direct. Consider G(K+5), where 0563}. 
Then G(K+b)=limy.. {1+|sin {1+0}%=1, where 
Lemma 1 is used to prove that in this case X — F(X) is the integer K. Then consider 
G(K+b) for }<b<1. We have G(K+B)=limy.. {1+|sin (K—1+2b)n| }¥ 
= limy {1+ | sin (2b7) | ©, since sin (267) for the assumed range of b. 

It is now possible, with the aid of the lemmas, to prove the theorem. We 
assume that X =K+6, where 0<b<1, so that [X]=K. Two cases are consid- 
ered. 

Case 1. 0Sb<}. Substituting X=K-+ 5 in the formula of the theorem, we 
find [X]=K+6—|b+2'"1—1| =K, where each step is justified by one of the 
lemmas. 

Case 2. }<b<1. Letting X =K-+6 in the formula of the theorem, we have 
=K. 

These two cases prove the theorem. 


Since the cost, in cents, for a letter of N ounces may be written C= —A[—N], 
where A is the cost per ounce, it is possible to write an analytical formula for 
this “postal function.” 


The writer wishes to express his appreciation for some very helpful sugges- 
tions on the part of the referee. 


ANOTHER PROOF FOR TWO PROPERTIES OF THE GENERATORS 
OF A RULED SURFACE 


E. K. Hstao, University of Chengchow, China 


C. E. Weatherburn [2] has proved the following theorem for ruled surfaces, 
referred to the generators v=constant and their orthogonal trajectories. 

When a is the unit tangent vector of the generators, so that div a is equal to 
the geodesic curvature of the orthogonal trajectories u=constant, y is the dis- 
tance function for the family of generators and K the Gaussian curvature, then 


a/1 
(1) = = My, (5) 2/(-—K) div a. 
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This theorem is obtained by using the Mainardi-Codazzi equations. Ram 
Behari [1] has proved it without reference to these equations. An alternative 
proof, which only depends on the material in [2], Sections 29, 30 is as follows: 

Suppose u denotes the actual distance along a generator; then we have 
E=1, F=0. Since the normal curvature in the direction of a generator is zero, 
it follows that L=0, so that the first and second curvature are given by: 


(2) J = N/G, 
(3) K = — M°/G. 


The distance function y for the family of generators v=constant is the 
reciprocal of the magnitude of Vv. Consequently, 


(4) ¥ =1/| Vol = 1/V(|r2|2/@) = VG, 
and 
(5) div a = G,/(2G) = yi/y. 


Since Wu+tyj=1, integration with respect to u gives WW,=u—a, and therefore 
y?=(u—a)?+6?, where a and £ are functions of v only. Now 


(6) K = — B*/y. 
On the basis of (5) and (6), we obtain immediately 
— 
and therefore that 
V/(-—K) = = — 2 diva. 


That is, 
Ou 


Ram Behari [1] also proved, now with the aid of the Mainardi-Codazzi 
equations, the following theorem for a ruled surface using an orthogonal co- 
ordinate system with v=constant as asymptotic lines, so that L=0, F=0. 


If the surface is ruled and the mean curvature varies inversely as the distance 
function along the generator, then the total curvature 1s constant along any orthogonal 
trajectory of the generators. 


Here the derivation is somewhat long and complicated. We point out that 
the result can be derived directly: Since ({2], Sec. 29) 


E 
a 
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for 
av au\My/ My du 


Consequently we have 
av 


from which the theorem follows. 
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CLASSROOM NOTES 
Eprtep By C. O. OaKLEy, Haverford College 


All material for this department should be sent to C. O. Oakley, Department of Mathe- 
matics, Haverford College, Haverford, Pa. 


THE HIGHER DERIVATIVE TEST FOR EXTREMA AND POINTS OF INFLECTION 
STANLEY W. Nasu, University of British Columbia 


1. Introduction. Consider a real function f of one real variable x. When the 
second derivative as well as the first vanishes at a point P: (xo, f(xo)), the second 
derivative test for determining whether P is an extremum of f fails to apply. 
Under these circumstances more advanced books on calculus give a higher 
derivative test for determining whether P is an extremum [1-9]. The proof al- 
ways seems to use a higher order mean value theorem, that is, Taylor’s theorem. 
Usually the first of the higher derivatives which does not vanish at xo, say f, 
is required to be continuous in an open interval (a, b) which includes x» [1-7]. 
However, the hypotheses of the test can be weakened: The existence and be- 


= 
a 
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havior of f® at any point but x» is irrelevant; f“-”) must exist throughout some 
open interval (a, b) including xo, but need not be continuous there. Occasionally 
the test is proved in this more general form by applying Taylor’s theorem to 
f*—, then using some extra argument [8, 9]. The present note proves the more 
general form of the test without using Taylor’s theorem, which the author feels 
is too sophisticated for the average first year student of calculus. 


2. Examples. Before beginning the proof we show that the extension of the 
higher derivative test is not trivial by finding functions whose sth derivative 
exists at xo but is discontinuous there. To these functions the extended test can 
still be applied at xo, although the usual test can not. Consider 


= x cos (1/x), x 0; = 0; 
where r=p/(20+1)>0 and p, o are integers. It is not too hard to show that 


(e) 


r— 1 r— 1 s 
(x) * cos (— — =) — sin (—- *) + O(« for x 0, 
x 2 x 2 


=0 #£forr>2s—1, 


where A,,>0 depends only on r and s, and O denotes the order of the other 
terms as x goes to 0. It follows that ¢® is defined at 0 but discontinuous there 
whenever 2s—1<r<2s. When r<2s, the oscillations of ¢ tend to infinity as x 
tends to 0. 

Consider now the functions 


f(x) = + cos (1/2), x0; f(0) = 0; 
g(x) = + cos (1/2), x #0; g(0) = 0; 


where 055=a/(28+1)<1 and a, B are integers. Then f’(0)=f’(0)= --- 
=f-)(0)=0 but f@(0)=(2n)!, and g’(0)=g’’(0)= - - -=g@"(0)=0 but 
g°"t+))(0) =(2n+1)!, while neither f@™ nor g"*» is continuous at 0. By Theo- 
rem 1 below f has a minimum at 0 and g has a point of inflection there. Thus 
Theorem 1 is a nontrivially stronger form of the higher derivative test for ex- 
trema and points of inflection. 

Two other examples may be of interest. Let Yo be an everywhere continuous, 
nowhere differentiable function, such as that of Weierstrass. Let y,(x) 


= for n=1, 2, - - - . Consider h(x) =xf,1(x). Then = 
+jy._;(x) for 7=1, +--+, s—1, while 4 (0) =s¥(0) but h(x) does not exist 
for x0. 
As a last example let u be defined by 
0 forx = 0 
u(x) =} for < |x| for integer n, 
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where r=p/(2o+1)>0 and p, o are integers. Then u’(x) is defined except when 
x=2-", and u’(0) =0. It follows that neither wu nor w’ is continuous in any open 
interval (a, b) which includes 0. The existence of its derivative at a point ensures 
the continuity of a function at the point, but not its continuity in any open inter- 
val including that point. : 


3. The higher derivative test. The higher derivative tests can be stated as 
follows: 


THEOREM 1. Let n denote a positive whole number, f a real function of one real 
variable. 
A. If f’ (xo) =f" (x0) = (x9) =0, then: 
(1) f&" (xo) >0 implies that f(xo) is a minimum value of f and that f is concave 
upwards at the point (xo, f(x0)); 
(2) f°" (xo) <0 implies that f(xo) is a maximum value of f and that f is concave 
downwards at the point (xo, f(x0)). 
B. If f’’ (x0) (x0) = (xo) =0, then: (xo) implies that 
f (x0) is not an extreme value of f, even if f'(xo) =0, but that (xo, f(xo)) is a point of 
inflection of f. 
(1) If f@"*» (xo) >0, then f is concave downwards to the left, concave upwards 
to the right of x». 
(2) If f@**) (xo) <0, then f is concave upwards to the left, concave downwards to 
the right of xo. 
(3) If in addition f'(xo)=0, then f°"*)(x9)>0O implies that f is increasing, 
and f'"+1)(x9) <0 implies that f is decreasing at the point (xo, f(x0)). 
First we prove Theorem 1-A for »=1. We assume as already known the 
following theorems [10]. 


THEOREM 2. If g is a function and g’ (xo) ¥0 exists, then there is an open inter- 
val (a, b) including xo such that g(x) exists for every x within the interval and 
[g(x) —g(x0) |/(x—xo) maintains the same sign as g'(xo) for every such x but xo. 

If g' (xo) >0, then g(u) <g(xo) <g(v) whenever 

If g' (x0) <0, then g(u) >g(x0) >g(v) whenever 


In applying Theorem 2 we shall consistently take (a, b) to be the most in- 
clusive open interval satisfying the conditions of that theorem. 


THEOREM 3. If g’(x) exists for every x in an open interval (a, b) and maintains 
the same sign throughout, then: 
(1) g’(x)>0 for every x in (a, b) implies that g is strictly increasing there and 
that g(a+) <g(u) <g(v) <g(b—) whenever a<u<v<b; 
(2) g(x) <0 for every x in (a, b) implies that g is strictly decreasing there and 
that g(a+-) >g(u) >g(v) >g(b—) whenever a<u<v<b. 


Suppose then that f’(xo) =0 and f(x) >0. Then, applying Theorem 2 to the 
function f’, f’(u) <0<f’(v) whenever a <u <xo9<v<b. It follows from Theorem 3 


= 
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that f(a+)>f(u)>f(xo) and f(xo)<f(v)<f(b—), that is, f(x)>f(x0), where 
xx is any other point in (a, b), and f(xo) is a local minimum value of f. For 
every x in (a, b) but xo, f(x) lies above the line y=f(xo), which is the tangent line 
to the graph of f at the point (xo, f(xo)). Hence f is concave upwards at xp. 
Similarly, if f’ (xo) =0 and f’’ (xo) <0, one can prove the second half of Theorem 
1-A for n»=1. Theorem 1-A with n=1 is commonly called the second derivative 
test for extrema. 

We now prove Theorem 1-A for general m. Suppose that f@"-")(x») =0 and 
f°” (xo) >0. Then, applying Theorem 2 to the function f@*-”, one finds that 
f°") (x) exists for every x in (a, b), the open interval, including xo, where 
(2) (x9) ]/(x—x0) >0 for every x but xo. The second derivative 
test is applied to the function f@*-” and one finds that f("-?) (x) > f("-®) (x9) =0, 
where xx» is any other point in (a, b). The same line of argument is used when 
f°" (xo) =0 and f@")(x9) <0. The existence of f@*-")(x) for every point x in 
(a, b) ensures the continuity there of f and its derivatives of order less than 
2n—1. 

The second derivative test can not be used further, for by hypothesis 


(xo) =0 for r=1, - - - ,2—1. To continue the proof we establish the following 
lemma. 


Lemna 4. Let h be a function for which h''(x) exists for every x in an open inter- 
val (a, b), and let h’ (xo) =h'' (xo) =0 for xo some point of (a, b). 

(1) If h’'(x)>0 for every x in (a, b) except xo, then, within (a, b), h’ is strictly 
increasing, h(x) is the minimum value of h, and h is concave upwards. 

(2) If h(x) <0 for every x in (a, b) except xo, then, within (a, b), h’ is strictly 
decreasing, h(x) is the maximum value of h, and h is concave downwards. 


Suppose that h’ (x9) =h’’ (xo) =0 and that h’’(x) >0 for every other x in (a, b) 
except x9. Then, applying Theorem 3 to the function h’, h’ is strictly increasing 
in the open intervals (a, xo) and (xo, 6) and h’(a+) <h'(u) <h'(x.—) =0 
=h'(xot+) <h'(v) <h'(b—) whenever a<u<xo<v<b. Thus h’ is strictly in- 
creasing throughout the open interval (a, 6). Since h’(u) <0<h'(v), it follows 
from Theorem 3 that h(a+)>h(u)>h(xo) and h(x) <h(v)<h(b—), that is, 
h(x) >h(xo), where xxo is any other point in (a, 6), and h(xo) is a local mini- 
mum value of h. For every x in (a, b) but xo, h(x) lies above the line y=h(xo), 
which is the tangent line to the graph of & at the point (xo, A(xo)). Hence h is 
concave upwards at xo. One can readily prove that h is in fact concave upwards 
for every x in (a, b). However, we do not need this to establish our main result 
and omit the proof of this detail. Similarly if h’ (xo) =h’’ (xo) =0 and if h’’(x) <0 
for every other x in (a, b) except xo, one can prove the second half of Lemma 4. 

We return now to the proof of Theorem 1-A for general n. Suppose we have 
established that f is continuous throughout the open interval (a, 6) and that 
f(x) >f@(xo) whenever By hypothesis (x9) (x9) =0. Then 
by Lemma 4, isstrictly increasing throughout (a, b), <0 (v) 
whenever a<u<x9<v<b, (x9) is the minimum value of within 
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(a, 6); and f°*-® is concave upwards throughout (a, b). Again, f("-® is continu- 
ous throughout (a, >f@*-®)(xo) there whenever x#xo, and by hy- 
pothesis f"—*) (x9) =f"-®)(x9) =0. If this argument can be started for some 
integer r=s, it can be continued for smaller integers until we reach r=1. But 
we earlier established the proper conditions for this argument for r=n—1. At 
the (n—1)st step, with r=1, we find that f’ is strictly increasing throughout 
(a, b), f’(u) <0<f’(v) whenever a<u<xo<v<b, and conclude that f(x») is the 
minimum value of f within (a, b), and f is concave upwards throughout (a, 3). 
Thus the first half of Theorem 1-A is established for general m. The same line 
of argument is used to prove the second half of Theorem 1-A for general n. 

To prove Theorem 1-B apply Theorem 1-A to the function f’. Suppose that 
the other proper hypotheses hold and that f@"*) (x9) >0. Then f”’(u) <0<f’’(v) 
whenever a<u<x9<v<b, as we found in the proof of Theorem 1-A, both 
when m= 1 and in the general case. Thus f has a point of inflection at the point 
(xo, f(xo)), and f is concave downwards in the open interval (a, xo), concave up- 
wards in the open interval (xo, 6). Furthermore f’(xo) is the minimum value 
of f’ in (a, b), and f’ is continuous throughout this interval. If f’(xo) 20, then 
f'(x) >0 for every other xx in (a, b). Then f’ does not change sign in (a, b), 
so f(xo) is not an extreme of f, but rather f is strictly increasing throughout 
(a, b). In particular this is true for f’ (x9) =0, proving the first half of 1-B-3. On 
the other hand, if f’(x9) <0, then by Theorem 2 there is an open interval (c, d), 
possibly a subinterval of (a, b), where [f(x) —f(xo) ]/(x—xo) <0 for every x in 
(c, d) but xo. Then f(u) >f(xo) >f(v) whenever c<u<xo<u<d, and again f(xo) 
is not an extreme value of f. Thus 1-B-1 is established. Similarly, with the other 
proper hypotheses and with f@"*) (xo) <0, one can prove 1-B-2 and the second 
half of 1-B-3. 
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A PROPERTY OF THE BERNOULLI NUMBERS 
L. Carxitz, Duke University 


The Bernoulli numbers in the even suffix notation satisfy the symbolic rela- 
tion 


(1) (B+ 1)"— (n> 1), 


where, after expansion, exponents are degraded to subscripts. By the Staudt- 
Clausen theorem, the denominator of B2, contains the prime 2 to the first power. 
Thus both numerator and denominator of 2B:, are odd after cancellation of 2. 
The more precise statement 

(2) 2Bon = 1 (mod 4) (n = 2) 


seems not to be very well known. It is implicit in a stronger result proved in 
Nielsen’s book ({2] p. 256). 


The congruence (2) can be proved rapidly as follows. From (1) we have 


2n 2 1 
(n > 0), 
r=0 
so that 
2 1 ae 1 
2— (2n+1)+— ) + (2n + 1)2Ba = 0, 
r=2 r 


Assuming that (2) holds up to and including the value »—1, we have 
2— (2n +1) — w + + (2n + 1)2Bon = 0 (mod 4). 
But 
bg ‘) {(1 + 1)24 4 (1 — 1)%#1} = 0 (mod 4) 
for n=1, so that 


(2m + 1)2Bon + 1 — 2(2n + 1) — 


Ill 


2n(2n + 2) = 0 (mod 4). 


(2n + 1)(2Bon — 1) = 2n + ; ‘) 


This evidently proves (2). 
By the same method we can prove that 


(3) 2Be, = 4n + 1 (mod 8) (n = 2), 


: 
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which is also implicit in a result given by Nielsen that is ascribed to Staudt. In- 
deed, since 4n-++1=(—1)"2—1 (mod 8), we have as above 


2—(Qn+1)+3 + ‘) + (2n + 


2 
= 0 (mod 8). 
But 
2n+1 2n n {2n+1 
ar )- >> )+ ar ) 
= (2m + 1){(1 + 1) — (1 — + (1 + + (1 — 


= 0 (mod 8) 


for n=2, so that 


x (4r + 1) (4n + + 1), 


2n+1 
1)2Bu +1 — (n+ 1) 2/ )- (4n + 1)(2n + 1) = 0, 


(2m + 1)(2Bon — 4n — 1) = 0 (mod 8). 
We recall that the numbers D2, defined by 


x 


sinh x (2n)! 
satisfy ({3] p. 28) Do, =2(1—22"-")B,,. Thus (2) implies 
(4) Den = 1 (mod 4), 
while (3) implies 


(5) Do, = 4n + 1 (mod 8); 


both (4) and (5) hold for all »=0. 
We remark that if 2" is the highest power of 2 dividing 2, then the writer 
has proved ([1] formula (10)) that 


(6) 2Bon = 1 (mod 2+) (n 2 2); 


this result is evidently in agreement with (3). 


References 
1. L. Carlitz, A note on the Staudt-Clausen theorem, this MONTHLY, vol. 64, 1957, pp. 19-21. 
2. N. Nielsen, Traité élémentaire des Nombres de Bernoulli, Paris, 1923. 
3. N. E. Nérlund, Vorlesungen iiber Differenzenrechnung, Berlin, 1924. 
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ON TWO FUNCTIONAL EQUATIONS* 
R. D. BoswELL, JR., Mississippi State University 


It is well known that the only real-valued continuous solutions of the func- 
tional equation 


(1) f(x + y) = f() + fO) 


are of the form f(x) = kx where & is a real number. In this note we shall solve two 
functional equations which may be given to a beginning analysis class as ap- 
plications of (1). Both of these equations arise naturally in describing the one- 
parameter groups of affine transformations in the euclidean plane. 


THEOREM 1. The only continuous solutions of the equation 


(2) + 9) = f(x) +f) + — A*)(1 AY), 
where a and A are real numbers and A>O, are of the form 
(3) f(x) = kx — a(i — A?) 


where k is a real number. 


Proof. Let us first find differentiable solutions of (2). Holding y constant and 
differentiating, we get 


+ 9) = f'(%) + al — A¥%)(— A? In A). 


Now holding x constant and differentiating, we obtain f” (x+y) =a(In A)?A*+», 
Let u=x+y; then f’(u)=a(In A)*A*. f’(u) =f9a(In A)2A ‘dt=a(In A)[A*]§ 
=a(In A)(A*—1) and f(u)=af9(In A) ]§. Therefore, 
f(u) =aA*—au(In A)—a. Let f; be any other continuous solution of (2); then 
g=f-—fi is also continuous. We notice that 


g(x + = + ») = f@) — A) + f0) — fi). 


This relation shows that g is of form (1) and, therefore, there exists a real num- 
ber ¢c such that g(x) =cx. Hence, any continuous solution of (2) is of form (3). 


THEOREM 2. The only continuous solutions of f(x+y) = A*f(y) +A"f(x), where 
A is a positive real number, are of the form f(x) =kxA* where k is a real number. 


Proof. Let g(x) = A-*f(x) and observe that g satisfies (1). Hence, there exists 
a real number k such that g(x) =kx. Therefore, f(x) =kxA?*. 


* This paper was presented to the Louisiana-Mississippi Section of the Mathematical Associa- 
tion of America on February 21, 1958. 
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ON THE SYNTHETIC DIVISION PROCESS 


W. W. Comrort, University of Washington anD Tuomas E. Mort, 
The Pennsylvania State University 


The authors present here a pair of proofs of a familiar result from college 
algebra. Stated without proof in-many reputable texts, the theorem, with 
either of the proofs here given, has proved generally palatable to those freshmen 
conversant with the axiom of mathematical induction. The content of the theo- 
rem is, very roughly, that if nonintegral entries arise in the standard synthetic 
division process, then the rational number being tested is not a root of the equa- 
tion in question. The result may be stated precisely as follows: 


THEOREM. Let p/q be a nonzero rational root of the polynomial equation f(x) 
- +--+ +a,=0, where each a, is an integer and (p, q)=1. Let 
bo=do, and, for each positive integer k not exceeding n, let 


be = a0(p/q)* + + + ae. 
Then each b; is an integer. 


First proof. Assume that some b,=1/s is nonintegral, where (r, s)=1. From 
the fact that 


by = [aop* + a,p*—'q axg*]/q*, 


it follows that s divides qg*, so that (p, s) =1. Then (pr, s) =1, so that (pr) /(qgs) 
=(p/q)b; is nonintegral. Then is nonintegral. 

(n—k)-fold repetition of this argument shows that 6, is nonintegral; in par- 
ticular, b, is nonzero. 


Second proof. In this proof we show that each of the numbers 0, is an integer 
divisible by q. 

That q divides the integer bp =a» is well known. Now let 0Sk <n, and sup- 
pose that }; is an integer divisible by g. Then bey: =(p/g)bi+Gey: is an integer 
which, in case k+1=n, is zero and hence is divisible by g. To show that g 
divides by41 in case R+1<n, set x=p/q in the equation f(x) =0 and multiply by 
q”*—! to get 


Since q divides the right-hand side, it divides the left-hand side, which is pre- 
cisely 


Since (p, q) =1, it follows that g divides by4:. 


y 

0 

> 

n 
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AN APPLICATION OF STIRLING’S FORMULA 
Davin ZEITLIN, Remington Rand Univac, St. Paul, Minnesota 


Stirling’s formula states that m!~(2rm)'/2n"e-*. Varying only in certain re- 
finements, the majority of proofs for Stirling’s formula use the relation 


b 
(1) f log x dx = blogs} —aloga—b+4a. 


For our purpose (given below), it suffices to know that there exists elementary 
proofs which, without using (1), establish Stirling’s formula in a weaker form, 
1.€., 


(2) ni ~ 


where C is an unspecified positive constant. Feller [1], Aissen [2], (there may 
be others) have given proofs of (2). The purpose of this note is to show how 
Stirling’s formula, (2), may be used to establish the validity of (1) for real num- 
bers, b>a>0. 

We recall that a continuous function is Riemann-integrable over a closed 
interval. For a>0, log x is integrable on aSxSb< &. Thus, if we divide the 
interval [a, b] into equal parts such that Ax,=(b—a)/n, x,=a+kAx,, 
k=0, 1, +--+, m, we are assured that 


b n 
(3) log adx = lim S, = lim {log (a + Ame. 


a k=l] 


Let a=r/s, b—a=p/q, where p, g, 7, and s are fixed positive integers. Let m 
be an arbitrary positive integer and define n=mqs(b—a) =mps. We note that 
p and s are fixed here, so that if m— ©, then n— ©. Now, Ax, = (b—a)/n=1/mgs 
and 


S, = (log (a + kAx)} Ax 


(4) mpe mns 
MS k=i = mgs 
(5) Tl (“= + 1 + mps) 
k=l mqs (mqs)?!4 (mqr) ! 


At this point, one introduces (proves) Stirling’s formula, (2). Thus, if we define 
Q(n) by the relation, m!=Cn"+!/2e-"Q(n), then lim,.. Q() =1. Using (2), we 
have 


(mar + mps)! _ Qlmagr mps)(mgr + (1 + 


| 
as 
© 
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r+mps 


ps 1/2 
Q(mgr) (r/s)™9" gr 
Substitution of (6) into (5) yields 
r (r/s)+(pla) / \ rls) 
= + (<) F(p, q, 7, m) 
(7) @ 5 
= e**b'a-*F (p, g, r, Ss, m), 
where 
mor + mbs 1/mqe 1/(2mqa) 
Q(mar) 
Let m— ©, Then F(p, g, 7, s, m)—1, n>, and S,—/? log xdx. Thus, from (7), 
b=a+(p/@) 
(9) f log xdx = blogb — aloga— b+ <a. 


Thus far, we have established the validity of (1) for rational numbers, 
b>a>0. Because of the continuity property of the function, log x, and because 
one may choose a= {r/ s} and b= {a+(p/q) } as two rational sequences of num- 
bers approaching two distinct, positive, irrational numbers as limits, respec- 
tively, the extension of the validity of (1) for irrational numbers is immediate. 
Thus, (1) is valid for all real numbers b>a>0. 


References 
1. W. Feller, An Introduction to Probability Theory and Applications, New York, 1950, 
pp. 41-42. 


2. M. I. Aissen, Some remarks on Stirling’s formula, this MONTHLY, vol. 61, 1954, pp. 687- 
691. 
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EpITtEp By Joun A. Brown, University of Delaware, AND 
Joun R. Mayor, AAAS and University of Maryland 


All material for this department should be sent to John R. Mayor, 1515 Massachusetts 
Avenue, N.W., Washington 5, D. C. 


GUIDE ON PURCHASE OF EQUIPMENT 


The Council of Chief State School Officers with the assistance of repre- 
sentatives from science, including staff members of the National Bureau of 
Standards, published in August, 1959 a very comprehensive guide for the pur- 
chase of equipment under the provisions of National Defense Education Act of 
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1958. The guide includes specifications on some 1000 items of equipment recom- 
mended for teaching science, mathematics, and modern foreign languages in 
elementary and secondary schools. There also is included an extensive bibliog- 
raphy for elementary and secondary school libraries. The project was in charge 
of Dr. Edgar Fuller, executive secretary of the Council. Advice on the equipment 
and bibliography recommended for use in the teaching of mathematics was given 
by a group of mathematicians and teachers whose meeting was made possible 
by the School Mathematics Study Group. 

Copies of the publication may be purchased from the state superintendent 
of public instruction in each of the states, or directly from Ginn and Company. 


AAAS STUDY ON THE USE OF SPECIAL TEACHERS OF SCIENCE 
AND MATHEMATICS IN GRADES 5 AND 6 


One of the major activities in the extended Science Teaching Improvement 
Program of the American Association for the Advancement of Science, under 
the new three-year grant to AAAS by the Carnegie Corporation of New York, 
is to be known as the Study on the Use of Special Teachers of Science and 
Mathematics in Grades 5 and 6. The cooperation of the public school systems 
in Cedar Rapids, lowa; Lansing, Michigan; Versailles, Kentucky; and Washing- 
ton, D. C. have been obtained for this project. In each of these school systems, 
at least one special teacher of science and at least one special teacher of mathe- . 
matics will teach only science or only mathematics in grades 5 and 6 to at least 
four different classes. The teachers have been selected by the school systems as 
those having special training, competence, and successful teaching experience 
in the areas to which they are assigned. 

While the use of special teachers of art, music, and physical education is 
fairly widespread in American elementary schools, the question of special teach- 
ers in the academic disciplines such as science and mathematics remains highly 
controversial. Majority opinion of experts in elementary education, both ad- 
ministrators and teachers, seems to favor strongly the retention of the “self- 
contained” classroom at this level. On the other hand, many scientists favor the 
use of specially prepared teachers particularly because of the new and antici- 
pated demands on the elementary teacher. It is quite impossible under the pres- 
ent methods of training elementary teachers for the elementary teacher to take 
what appears to be a necessary minimum of course work in science and mathe- 
matics, or in the other academic disciplines as well. The National Science 
Foundation sponsored curriculum studies, which are giving attention to the 
science curriculum at the elementary level, will stress how serious is the prob- 
lem. Both the School Mathematics Study Group and the Biological Sciences 
Curriculum Study are concerned with grades from kindergarten to grade 8, as 
well as with the high schools, and there are also other privately sponsored ac- 
tivities which are working on materials at the elementary level. 

An attempt will be made to provide a careful evaluation of the Study so 


ee 


| 
| 
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that interested persons will have available objective data to support their cases 
or to persuade them of weaknesses in their present point of view. It is antici- 
pated that the Study will be conducted during a two-year period beginning in 
the fall of 1959 and it is to be hoped that the Study might be extended to more 
schools during the second year. (Abstracted from Science Education News, Sep- 
tember 1959.) 


A NEW KIND OF COURSE FOR THE PREPARATION OF TEACHERS* 
E. P. Norturop, College of the University of Chicago 


Three main streams can be discerned in current efforts to improve secondary 
school mathematics curricula in this country. First, various groups, working 
at both national and local levels, are formulating proposals for curricular change. 
Second, various groups, also working at both national and local levels, are pre- 
paring new classroom materials. Third, numerous programs, supported by in- 
dustry and by private and public foundations, are being offered to improve the 
mathematical and professional competence of teachers. 

These efforts have come into being and have grown too large too fast to be 
well coordinated. Yet their interrelatedness is readily apparent. If the proposals 
for curricular change are to be understood and not rejected out of hand, then 
teachers need to be trained in the relevant mathematical subject matter. If the 
proposals are to be acted upon and not just read and filed away, then teachers 
need to be supplied with relevant text materials. 

These efforts, however well conceived and executed, will not in themselves 
suffice. What is lacking? The more thoughtful of the groups proposing curricular 
change agree that no single curriculum will serve all the students in all the 
schools. Accordingly, when they produce a particular curriculum by way of il- 
lustration of what might be done, they are quick to emphasize that this is only 
an example, and that the teachers of the country must adapt the suggested 
changes to their local needs. Again, the more thoughtful of the groups preparing 
new text materials similarly reject the notion of a universal panacea. Accord- 
ingly, when they produce texts by way of illustration of what might be written, 
they too are quick to emphasize that these are only examples, and that the 
teachers of the country must adapt the suggested materials to their local needs. 

What is lacking now begins to emerge in the picture, as does the place where 
the lack can be remedied. Overburdened classroom teachers, however well 
trained or retrained, are about to be faced with a bewildering and often con- 
flicting variety of curricular proposals and text materials. They need help in 
selecting and organizing materials for use in their own local classrooms. They 
need, in short, help in learning to think and to write for themselves. The help 
they need can and should be provided in the third stream of our current efforts, 
namely, in the programs of re-education being offered to inservice teachers. 


* Presented at a panel discussion, “The Training of Secondary School Mathematics Teachers,” 
Annual Meeting of the Mathematical Association of America, Philadelphia, Pa., January 23, 1959. 
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Let us therefore consider a possible design for a seminar of nine months’ 
duration, intended to provide the individual teacher with a working understand- 
ing of curriculum construction, together with working experience in the prepara- 
tion of materials for use by his own students. It is assumed that the seminar is 
to be attended by twenty-five or thirty teachers, and is to be conducted con- 
currently with, or subsequent to, appropriate courses in mathematics proper. 

During the first three months of the year, and while the teacher is learning 
some of the mathematics he will need, the seminar would be devoted to a down- 
to-earth consideration of factors that have an important bearing on curriculum 
construction. Examination of a few specific curriculum proposals, including the 
bases offered in support of them, would lead naturally and quickly to a discus- 
sion of such factors as the needs or demands of society, the philosophy of educa- 
tion of the school, the characteristics of its students, the nature of the subject- 
matter, the formulation of aims or objectives, the selection and organization of 
learning activities, and testing procedures. Such a discussion, properly con- 
ducted, can serve three purposes. First, it can bring order into possibly dis- 
ordered and unrelated ideas gleaned from earlier courses in education. Second, 
it can provide the teacher with a framework for the writing efforts to follow. 
Third, it can provide him with a rationale by which to justify his proposals and 
sell his product to an audience too often ignored, namely, the folks back 
home. 

The next two months would be devoted to group writing efforts. The class 
would be divided into groups of four or five teachers each, sorted according to 
expressed interest in this, that, or the other general topic. Each group, having 
selected a modest and manageable unit, would settle upon the kind of student 
for which it is intended and what he is expected to achieve from it, would outline 
it, and would write a piece of text complete with problem material. It would be 
made clear at the outset that these initial efforts were not expected to terminate 
in a polished product—that they were simply practice for a more extensive 
effort to come. The seminar would continue to meet, though perhaps less fre- 
quently, for open discussion of progress, major difficulties encountered here or 
there, and so on. The final products would be reproduced and distributed to all 
members of the seminar for analysis and criticism, in the course of which all 
would learn something about the importance of correctness, clarity, and preci- 
sion in mathematical exposition. But there is an additional and very real value 
in this preliminary writing. Here the teacher finds himself, perhaps for the first 
time, facing up to issues and committing himself to paper, and by doing so in 
the company of his fellows he gains as much in confidence as he gains in simple 
experience. 

The last four months of the year would be devoted to the major effort and 
product of the seminar. Each teacher would embark on an individual project 
in the planning of a year course for use in his own high school by his own stu- 
dents. The project would issue in a substantial report consisting of the following 


parts: 


- 
- 
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1. A brief general description of the course. 

2. Its justification, addressed to school board and principal, parents and 
fellow teachers, and including general statements of aims, of shortcom- 
ings of the present program, and of strengths of the program being pro- 
posed. 

3. A full outline of the year’s work. 

4. A complete treatment of a substantial unit of the course (for example, a 
month’s work, including classroom materials to be used by the student 
for further study), an examination over the unit, and a teachers’ manual 
for use by other teachers. 


During this period the seminar would be disbanded as such, though it could 
conceivably meet for discussions with visiting representatives of outside groups 
concerned with curricular proposals or preparation of text materials. Members 
of the staff would of course make themselves available for individual conferences. 
It would be desirable, if it could be arranged, to use the fourth and last month 
of the period for criticism and revision of the reports, especially if they were to 
be used—as well they might—in partial fulfillment of certain degree require- 
ments. 

A seminar of the kind described is perhaps most suitable for use with teach- 
ers who have had several years of teaching experience.{ It could nevertheless 
be adapted—no doubt with lowered sights—for use with new teachers in the 
course of their practice teaching. Indeed, a foreshortened version of the seminar 
might be imagined for use in a summer program of, say, eight weeks’ duration. 
Here the work would take the form of group projects, with a skilled leader as- 
signed to each of the groups of four or five teachers for two purposes: to direct 
the writing efforts and, in the course of them, to call attention to the curricular 
issues involved. In this way the first two parts of the nine-month seminar de- 
scribed above might be successfully telescoped, especially if the several groups 
concerned were to meet together once a week, under a director of the entire 
operation, for discussion of current proposals for curricular change. 

It is to be anticipated that the individual or group reports issuing from semi- 
nars of the kind described in this paper would be of mixed quality. But it can 
be safely asserted that the author or coauthor of even the weakest report would 
have grown enormously in critical understanding of his task. It is often said 
that one does not really know his subject until he tries to teach it. In the same 
vein, one does not really know his subject or his teaching until he tries to write 
what he teaches. Even if what he writes is not worth teaching, the experience 
will have made him a far better and wiser judge of courses and texts that others 
write for him to use. 


+ Evidence that it can be done is furnished by the fact that in the 1958-59 Academic Year 
Institute for High School Mathematics Teachers at the University of Chicago, Mr. Maurice 
Hartung and the author are in midstream with a second such seminar. The author's indebtedness 
to Mr. Hartung, and to other colleagues in education and mathematics who assisted in the design 
and conduct of the seminar, is hereby gratefully acknowledged. 
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ELEMENTARY PROBLEMS AND SOLUTIONS 
_ Epitep sy Howarp Eves, University of Maine 


Send all communications concerning Elementary Problems and Solutions to Howard 
Eves, Mathematics Department, University of Maine, Orono, Maine. This department wel- 
comes problems believed to be new, and demanding no tools beyond those ordinarily furnished 
in the first two years of college mathematics. To facilitate their consideration, solutions should 
be submitted on separate, signed sheets, within three months after publication of problems. 


PROBLEMS FOR SOLUTION 
E 1381. Proposed by H. E. Hart, City College of New York 
Show that 
+ 22/21 + 28/31 + + = 0 


has no real roots. 


E 1382. Proposed by Ian Connell, University of Manitoba 
Let u, be the mth term of the sequence 
1, 2, 4, 5, 7, 9, 10, 12, 14, 16, 17,---, 
where one odd number is followed by two evens, then three odds, etc. Prove that 
Un = In — [(1 + V8n — 7)/2], 
where square brackets denote the integral part. 


E 1383. Proposed by E. J. Burr, University of New England, N.S.W., Aus- 
tralia 


Let P, be the problem: “Find a set of »>2 consecutive positive integers, 
the greatest of which is a divisor of the least common multiple of the remaining 
n—1.” Show that precisely one of the problems P, has a unique solution. 


E 1384. Proposed by J. H. Butchart, Arizona State College 
Construct a circle through two given points, separated by a given circle, 
which shall cut the given circle at the smallest possible angle. 


E 1385. Proposed by John Burke, Saranac Lake, N. Y. 


A homogeneous die is “loaded” at a single point and is cast ” times. The 
one spot appears m times, the two spot m2 times, - - - , the six spot mg times. Find 
the (probable) position of the load. 
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SOLUTIONS 
Generalization of a Russian Olympiad Problem 


E 1337 [1958, 708; 1959, 424]. Proposed by the Student-Faculty Colloquium, 
Carleton College 


What is the largest value of m, in terms of m, for which the following sentence 
is true? If from among the first m natural numbers any 7 are selected, among 
the remaining m—zx at least one will be a divisor of another. 

II. Solution by J. L. Selfridge, IB M Research Center, Yorktown Heights, N. Y. 
The published solution [1959, 424] contains an oversight. The proof is correct 
only for m<12. If m212 and we select n= [m/2]—1 numbers, the remaining 
m—n numbers need not contain 1 or (Rk, 2k) as asserted by the solver. This is 
easily shown by selecting 1 and all numbers exactly divisible by an odd power 
of 2. For example, when m= 12 select 1, 2, 6, 8, 10. The remaining numbers are 
odd or an odd multiple of 4, and no one is twice another. Following is believed 
to be a correct solution of the problem. 

If 1, 2,---+, [m/2] are selected from 1, 2, - - +, m, then none of those re- 
maining divides another of them. This shows that » < [m/2]. Now we show that 
for any nS [m/2]—1, some a and bd will remain unselected with b= 2*a. Suppose 
that no such a and 6 exist. Then we must select all but one from each sequence 
r, 2r, - ++, 2*r. These sequences are disjoint if r is odd, and we count those 


selected by noticing that selecting all but the largest in each sequence means 
selecting 1, 2,---, [m/2]. 


Locus of Centers of a Family of Equilateral Triangles 
E 1351 [1959, 141]. Proposed by J. F. Darling, Woodstown, New Jersey 


Find the locus of the centers of the equilateral triangles whose sides a, }, c 
pass through three fixed points A, B, C respectively. 

Solution by Helen M. Marston, Educational Testing Service, Princeton, N. J. 
Construct the three circles O., Os, O. for which the segments BC, CA, AB 
respectively subtend 120° arcs, the minor arc in each case being directed toward 
the inside of triangle ABC. From any point P, on O, draw P.AP, (line a) and 
P.BP, (line b). Angles P.P,C and P.P,C each equal 60° by construction, whence 
P;, C, P. are collinear (line c) and triangles P,P,P,. represent one set of equi- 
lateral triangles abc. Let M., Ms, M. be the midpoints of the minor arcs BC, 
CA, AB respectively. The center M of triangle P,P,P, is the intersection of 
P.M, and P,M;, which must intersect at 60°. Since the locus of M includes 
points M,, Ms, M., these three points form an equilateral triangle and the locus 
of M is the circle circumscribing it. 

If circles O/ , Of, Of are constructed similarly but with minor arcs directed 
toward the outside of triangle ABC, we obtain a second set of equilateral tri- 
angles abc, a second set of minor-arc-midpoints M/, Mg, M?, and a circle cir- 
cumscribing them for the locus of M’. 

The required locus is therefore a pair of circles. It can be shown that they 


| 
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have a common center G, which is the centroid of triangle ABC. In the limiting 
case where A, B, C are collinear, the two circles coincide. When A, B, C form 
an equilateral triangle, there is only one set of equilateral triangles abc and the 
locus of M is the single point G. 


Also solved by Leon Bankoff, W. B. Carver, J. W. Clawson, Michael Goldberg, L. D. Gold- 
stone, W. R. McEwen and Amos Nannini (jointly), D. C. B. Marsh, Beckham Martin, D. C, 
Stevens, and the proposer. 

The loci of M and M’ may be alternatively described as the circles having centers at the 
centroid of triangle ABC and passing through the isogonic centers R and R’ of triangle ABC, or 
as the circumcircles of triangles Of Of O! and 0,0,0.. 

Goldberg, Marsh, and Martin used cartesian coordinates; Carver, Clawson, and the proposer 
used conjugate coordinates. Goldstone pointed out the connection between the problem and Prob- 
lem 4301 [1950, 565]. Martin considered the more general problem of finding the locus of the 
centroids of all triangles of a given species whose sides a, b, c pass through three fixed points A, B, C 
respectively. 


The Case of the Careless Typesetter 
E 1352 [1959, 141]. Proposed by C. W. Trigg, Los Angeles City College 


In setting the type for the multiplication (abc) (bca)(cab) = 234235286, 
a>b>-c, in which the unit’s digit is 6, the remaining digits of the product be- 
came pied. Restore them to their proper order. 

Solution by W. R. Talbot, Lincoln University, Missouri. Let m=(abc) and 
let N be the desired product. To establish limits on m, we notice that if c=1, 
then the largest m, namely 981, gives VN = 159,080,922, which is too small. There- 
fore 988>n>981. But the sum of the digits of N is 35=2 mod 3; whence n 
and its permutations are congruent to 2 mod 3. Then is 986 or 983, of which 
only 983 has a product of digits ending in 6. It is now readily found that 


(983) (839) (398) = 328,245,326. 


Also solved by A. N. Aheart, Merrill Barnebey, D. A. Breault, W. E. Buker, J. M. Calloway, 
G. B. Charlesworth, P. L. Chessin, R. J. Cormier, Eleanor G. Dawley, Monte Dernham, J. E. 
Faulkner, J. F. Foley, Michael Goldberg, L. D. Goldstone, H. W. Gould, Robert Gregorac, E. E. 
Griffee, P. G. Hodge, Jr., J. H. Hodges, J. M. Howell, A. R. Hyde, F. S. Innis, Jr., Sister Kenneth 
Kolmer, W. R. McEwen, D. C. B. Marsh, Helen M. Marston, J. H. Means, Franklin Mohr, 
D. L. Muench, J. B. Muskat, D. J. Persico, J. P. Phillips, J. L. Pietenpol, C. F. Pinzka, B. E. 
Rhoades, J. A. Rodriguez, J. L. Selfridge, D. C. Stevens, W. B. Stovall, Jr., Harry Weingarten, 
Mildred West, Charles Wexler, Dale Woods, Ronald Wyllys, and the proposer. 


Eccentricity of an Elliptical Section of a Cylinder 
E 1353 [1959, 141]. Proposed by L. A. Kenna, University of Arizona 


What is the eccentricity of an ellipse formed by a plane cutting the axis of 
a right circular cylinder at an angle of 6? 


Solution by A. N. Aheart, West Virginia State College. Denoting the lengths 
of the major and minor axes by 2a and 26, we have 


= (1 — = (1 — sin? 9)? = | cosa. 


xis of 


ngths 
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Also solved by R. G. Albert, E. F. Allen, J. W. Baldwin, Leon Bankoff, Merrill Barnebey, 
H. F. Bechtell, A. P. Boblétt, D. A. Breault, C. S. Carlson, G. B. Charlesworth, P. L. Chessin, 
A. G. Clark, R. J. Cormier, C. H. Cunkle, Alfred Diarnbey, Michael Difon, C. W. Dodge, Ernest 
Enns, G. W. Erwin, Jr., J. E. Faulkner, Susan L. Friedman, Todd Gitlin, George Glauberman, 
Michael Goldberg, L. D. Goldstone, F. D. Grogan, J. D. Haggard, R. Henry, Richard Holt, A. F. 
Horadam, R. H. Hou, A. R. Hyde, G. S. Innis, Jr., M. S. Itzkowitz, J. D. E. Konhauser, W. E. 
Lawrence, W. R. McEwen, D. C. B. Marsh, Beckham Martin, C. S. Ogilvy’s freshmen calculus 
class, D. J. Persico, J. L. Pietenpol, C. F. Pinzka, B. E. Rhoades, L. A. Ringenberg, D. A. Robin- 
son, J. A. Schumaker, Ruth G, Smith, D. C. Stevens, W. B. Stovall, Jr., W. R. Talbot, C. W. Trigg, 
T. C. Wales, Charles Wexler, L. H. Williams, R. H. Wilson, Jr., Alan Wofford, Dale Woods, J. W. 
Young, and the proposer. 


Editorial Note. The eccentricity of the conic section formed by a plane cutting at an angle @ 
the axis of a right circular cone of generating angle ¢ is |cos 6| /cos ¢. 


A Property of the Coefficients of (1-+-x)?~? 
E 1354 [1959, 141]. Proposed by P. L. Chessin, University of Maryland 


If --- +a,-2x?-*, where p is a prime, then 
a,+2, a2—3, a3+4, are all multiples of p. 


Solution by C. F. Pinska, University of Cincinnati. Multiplying correspond- 
ing members of the k congruences p—1= —i (mod p), 1=2, 3, - -- , k+1, gives 


(p — 2)\(p — k — 2)! = (—1)*(R + 1)! (mod 9). 
For k<p we have (k!, p) =1 and thus 
(p — 2)!/(p — k — = (—1)*(k + 1) (mod 9), 
demonstrating that a,+(—1)***(k+1) is divisible by p. 


Also solved by A. N. Aheart, R. G, Albert, J. W. Baldwin, H. F. Bechtell, R. F. Brown and 
Anna Endelman and Joel Levy (jointly), J. M. Calloway, G. B. Charlesworth, A. E. Danese, 
Alfred Diarnbey, H. B. Emerson, J. E. Faulkner, N. J. Fine, Todd Gitlin, George Glauberman, 
Michael Goldberg, L. D. Goldstone, Joseph Hammer, Leonard Hauer, Richard Holt, J. H. Hodges, 
Irving Katz, Jack Klugerman, Sidney Kravitz, W. E. Lawrence, Robert McGuire, D. C. B. Marsh, 
Andrzej Mgkowski, Valerie Miké, Mary P. Moseley, D. L. Muench, Stewart Nagler, K. K. Norton, 
C. S. Ogilvy, D. J. Persico, Benjamin Sapolsky, L. J. Schneider, Arnold Singer, D. C. Stevens, 
H. R. Stevens, Gerald Stoller, W. R. Talbot, W. A. Veech, T. C. Wales, Charles Wexler, Dale 
Woods, Ronald Wyllys, J. W. Young, David Zeitlin, and the proposer. 


A Detail of a Larger Work 
E 1355 [1959, 141]. Proposed by A. J. Goldman, National Bureau of Standards 


If y(x) >0 and y=log [1+ (log x)/x+y/x] for all sufficiently large x, prove 
that y~(log x)/x as x. (This is a step whose details are left to the reader 
in a paper by S. Chowla and F. C. Auluck, Some properties of a function con- 
sidered by Ramanujan, J. Indian Math. Soc., vol. 4, 1940, pp. 169-173.) 


Solution by David Zeitlin, Remington Rand UNIVAC, St. Paul, Minn. From 


1+ (log x)/x+ y/x = >1+y 
we have 


ng 
rm 
the 
old- 
the 
Or 
oser 
the 
286, 
be- 
and 
= 
lere- 
ce n 
hich 
oway, 
J. 
E. E. 
nneth 
Mohr, 
B. E. 
‘arten, 


728 ADVANCED PROBLEMS AND SOLUTIONS [October 
0<y < (log x)/(x — 1), 
whence y(x)—0 as ©. Therefore 


lim xy/(log x) = lim y/(e¥ — 1) = 1. 


Also solved by D. A. Breault, R. F. Brown and Anna Endelman and Joel Levy (jointly), 
N. J. Fine, George Glauberman, Walter James, D. C. B. Marsh, Morris Morduchow, K. K. 
Norton, C. F. Pinzka, D. C. Stevens, Peter Treuenfels, and W. A. Veech. 


ADVANCED PROBLEMS AND SOLUTIONS 
EpITED By E. P. Starke, Rutgers University 


Send all communications concerning Advanced Problems and Solutions to E. P. Starke, 
Rutgers University, New Brunswick, New Jersey. All manuscripts should be typewritten 
with double spacing and margins at least one inch wide. Problems containing results believed to 
be new or extensions of old results are especially sought. Proposers of problems should also en- 
close any solutions or information that will assist the editor. In general, problems in well- 
known textbooks or results in readily accessible sources should not be proposed for this de- 
partment. 


PROBLEMS FOR SOLUTION 
4863. Proposed by W. F. Stinespring, Institute for Advanced Study 
For what constant k (if any) can the inequality 
|4+B| 4| + 


be established for all ” by m real symmetric matrices (n= 2)? Here the notation 
| A| means the positive square root of A*, and the relation S between two 
matrices means that their difference is nonpositive definite. 


4864. Proposed by D. J. Newman, Massachusetts Institute of Technology and 
AVCO Research and Development 


Let a coin be tossed repeatedly, and let p be the probability that at some 
point the number of heads exceeds twice the number of tails. Prove 


p= — 1). 
4865. Proposed by L. Lewin, Enfield, England 
Defining the inverse tangent integral of the second order by 
0 x F 
provethat 6T7%2(1) — 47i2(1/2) — 2Ti2(1/3) — = log 2. 
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4866. Proposed by F. D. Parker, University of Alaska 


Show that a semi-magic matrix A (the sums of the rows and columns are 
all equal) can be decomposed into a sum B+C such that for integral K, 


(B+ C)*¥ = BE + C¥, 
4867. Proposed by Ky Fan, Oak Ridge National Laboratory 
Let f be a measurable function such that 0<f(x) <1/2 on [0, 1]. Prove 


exp f log f(x)dx exp f log (1 — f(x))dx 


4868. Proposed by Ronald Pyke, Columbia University 


Show that for all 8 Se and for all real x the function f defined by f(0, 8B) =1 
and 


= (x #0) 


is finite and satisfies the relationship 
S(%, B) = 


SOLUTIONS 


Matrices Having a Given Multiplication Table 
3099 [1924, 455]. Proposed by S. A. Corey. 


Let H, I, J, and K be the complex quaternion units, 
a1 dB}, 
respectively, which have the multiplication table 


H H I -!I 
I H I 
J J-K JE 


Ki-J K J 


Prove that there exists a set of four matrices, involving no imaginaries, which 
have the same multiplication table. 
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Solution by D. C. B. Marsh, Colorado School of Mines. Using known cor- 
respondences of a+i+cj+dk and x+y to 


4 
—-b a-d ¢ 
and 
—d@d-c 


respectively, we can establish composite 8 by 8 matrices with real components 
which have the required multiplication table. One set follows: 


10-1 0041041 @ 1 
101001 0-1 6% 

H= = 

4) o-1 0-1 1 0-1 0 6 
@ 4 4-06 64 8a 

Linear Independence 


4797 [1958, 452; 1959, 322]. Proposed by D. J. Newman, Massachusetts 
Institute of Technology and AVCO Research and Development 


Prove that all expressions like 7\/19/4—3.\/7+8-/6/5 are irrational. More 
specifically, prove that the square roots of the square-free integers are linearly 
independent over the rationals. 

II. Editorial Note. P. T. Bateman has raised the objection to the solution 
printed in the April issue that the c; appearing in formula (3) may be all zero 
and that a proof to the contrary is not given. A similar oversight seems to exist 
in each of the other contributions to this problem, so that we have no solution 
available at this time. 

Bateman also notes that the result of the problem is contained in Theorem 2 
and its Corollary in A. S. Besicovitch, Linear independence of fractional powers of 
integers, Proceedings of the London Mathematical Society, XV (1940), pp. 3-6. 
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Reducibility of a Polynomial ° 
4815 [1958, 713]. Proposed by N. S. Mendelsohn, University of Manitoba 


Let a:,-- +, a, be distinct integers. Discuss the reducibility of the poly- 
nomial f(x) =(x—a) - + - («—a,)+1 over the field of rationals. 


Remarks by Hans Rademacher, University of Pennsylvania. This problem is 50 
years old and due to I. Schur. A solution is printed in Pélya-Szegé, Aufgaben 
und Lehrsdtze, vol. 2, VIII, pp. 346-7. Many generalizations have been dis- 
cussed in the meantime. References can be found in Brauer and Ehrlich, On the 
irreducibility of certain polynomials, Bull. Amer. Math. Soc., vol. 52, 1946, pp. 
844-856. 


Also solved by P. T. Bateman, J. H. Hodges, D. C. B. Marsh, W. F. Trench, and the proposer. 
The result is: the polynomial is irreducible except when n=2 and a; and az differ by 2, and when 
n=4 and 4, de, ds, a4 are in arithmetic progression with common difference 1. 


An Integral Equation 


4816 [1958, 779]. Proposed by M. S. Klamkin, A VCO Research and Develop- 
ment, Wilmington, Mass. 


Solve the integral equation 


where a and b are independent of x. 
Solution by Robert Weinstock, University of Notre Dame. Introducing the 
change of variable u=t—x and defining F by 


(1) F(x) = f Po(x — tdi = f (u + x)*o(—u)du, 
0 —z 
we may differentiate the integral equation and obtain 


+ x)*¢(—u)du = ab} f + + x)o(—u)du, 


that is 
3 = ab[F(x)]*-°F'(x) = ab(d/dx){(b — 1)—[F(x)]*“}, 


provided b¥1; the case b=1 is handled separately below. We introduce 
m=b/(1—b) and integrate to obtain 


(2) F(x) = [¢ — (3/om)x}->, 


where ¢ is an arbitrary constant. 
Differentiating (1) three times and using (2) we obtain 


(x) = = $(3/am)*(m + 1)(m + 2)(m + 3)[c — 


r 
| 
\ 
| 


732 ADVANCED PROBLEMS AND SOLUTIONS [October 


For the existence of the integrals involved we must have m>0 (whence 0<b<1), 
a>0, and x <amc/3. 

In case b=1, the differential equation for F reads F’(x) = (3/a) F(x), whence 
F(x) =2(a/3)*Be*/*, where B is an arbitrary constant. We then have ¢(x) 
= 4F’"' (x) = Be*/*, which satisfies the integral equation for b=1, provided a>0, 
for all x. 

If a=0, there is clearly only the trivial solution ¢(x) =0. 


Also solved by H. E. Fettis, Emil Grosswald, C. C. Yalavigi and the proposer. 


Entire Functions 
4817 [1958, 779]. Proposed by I. S. Gél and Stanley Kaplan, Cornell Univer- 
sity 
Determine all entire functions f(z) such that |f(z)| =1 whenever |z| =1. 


Solution by C. Goffman and M. Henriksen, Purdue University. Write f(z) 
=2"g(z), where 7 is a nonnegative integer, g is entire, and g(0) #0. Since | z| =1 
implies | g(z)| =1, either g is constant, or the open mapping g sends the unit 
circle onto itself. In the latter case, the Schwarz reflection principle tells us that 
if s: and z are collinear with the origin and | ziz2| =1, then 


| g(zs)g(z2) | = 1. 


(cf., e.g., Bieberbach, Conformal Mapping, New York, 1953, p. 76). It follows 
that g is bounded on a neighborhood of infinity. So, by Liouville’s theorem, g is 
a constant. 


Hence f(z) =az", where |a| =1. 

Also solved by I. N. Baker, P. T. Bateman, Harley Flanders, D. S. Greenstein, Benoit La- 
chapelle, W. F. Trench, and the proposers. 

A Property of the Miquel Circle of a Complete Quadrilateral 

4818 [1958, 779]. Proposed by Hiiseyin Demir, Zonguldak, Turkey 

Let d; be the sides of a complete quadrilateral, and A,; be the vertex on 
d;, d;. Let t; be the triangle formed by the sides other than d;, and (O,;) denote 
the circumcircle of ¢;. Denote the Simson line of a point S; of (O;) with respect 
to by Dy. 

Then prove that, if D; and d; are parallel for all 7, (1) the line S,O, passes 


through the vertex Ag, (i, p¥q, r), and (2) the points S; all lie on the Miquel 
circle (O). 


Solution by the proposer. (1) Let the projections of S; and O, on d, be denoted 
by U;,, Vor. Then, using directed angles, we have 
XV prOpAgr = X AirAgiAg (from the circle (O,)), 
= £ UiUigAg (since D; is parallel to d,), 
=f (from the circle ig). 
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Now, since O,V,, is parallel to S;U;, we get the required collinearity of S;, 
0,, Ag. 

(2) Since the line of centers 0,0, is perpendicular to the radical axis FA pq, 
where F is the Miquel point, we have successively 


x 0,0,0, AgiF Api, 
= (from the circle (0,)), 


= 
= pr (from the circle (0,)), 
= £ 0,50, (from property (1)), 
and S; lies on the Miquel circle (QO). 
Also solved by A. E. Landry. 
Sum of Five Distinct Squares 


4819 [1958, 779]. Proposed by P. T. Bateman, University of Illinois 


Suppose k& is a given nonnegative integer. Show that every sufficient large 
positive integer is a sum of five squares none of which is less than k? and no two 
of which are equal. 


Editorial Note. After having submitted an independent proof, the proposer finds stronger re- 
sults obtained by E. M. Wright, using Vinogradov’s method. [Quart. J. Math. Oxford Ser. (2), vol. 
4, 1933, pp. 37-51.] Wright showed that every large positive integer can be expressed as a sum of 
five squares whose ratios are roughly prescribed in advance. This certainly implies the present 
theorem. 

Also solved by L. Carlitz, using a result of Walfisz’. 


RECENT PUBLICATIONS 
EpITED By RIcHARD V. ANDREE, University of Oklahoma 


All books for review should be sent directly to R. V. Andree, Department of Mathematics, 
University of Oklahoma, Norman, Oklahoma, and not to any of the other editors or officers of 
the Association, 


Mathematics in Business. By Lloyd L. Lowenstein. Wiley, New York, 1958. 
xv+364 pp., $4.95. 


“This text is for use in a first course in the mathematics of business. The in- 
tention is to give the student a firm foundation for further courses in the mathe- 
matics of finance as well as accounting, business statistics, insurance, and other 
business subjects.” The book should be successful in accomplishing the stated 
aims of the author. Professor Lowenstein has a gift of writing in clear and 
simple, as well as concise, English. Apparently the students for whom this book 
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was written are those without exceptionally strong backgrounds in high school 
mathematics, and although the degree of rigor attained is superior to that in 
many similar books, the approach is not too subtle for his audience. 

Chapters I and II review arithmetic and beginning high school algebra. This 
material is not an integral part of a course for those with normal preparation. 
Logarithms and the Binomial Theorem are treated in appendices whereas pro- 
gressions are taken up when needed in the body of the text. The business mathe- 
matics chapters cover the topics: percentage; profit and loss; simple interest and 
discount; compound interest; depreciation; annuities with applications to 
amortizations, sinking funds, and bonds. The book is not intended to be a 
treatise on the theory of interest and Professor Lowenstein has achieved a splen- 
did balance between theory and practice. More attention is given to simple 
interest methods than in some books. 

Unusual features of the book are the listing of the answers to the odd- 
numbered exercises at the end of each chapter and very effective use of footnotes 
for clarification of theory. The text includes 61 pages of the usual financial 
tables. The logarithm table is only a five-place table but the logarithms of inter- 
est rates are given to seven decimal places. 

C. L. SEEBECK, JR. 
University of Alabama 


Experimental Designs in Industry. Edited by Victor Chew. Wiley, New York, 
1958. 268 pp. $6.00. 


This is a collection of some of the papers given at a symposium held Novem- 
ber 5-9, 1956, conducted by the Institute of Statistics of the University of North 
Carolina, State College Section, and sponsored by the Mathematics Division 
of the Air Force Office of Scientific Research. 

The papers are of a technical, statistical nature but use relatively elementary 
mathematics. Up-to-dateness is the keynote. Parallel treatments of experimental 
design in general, factorial designs, and designs for exploring response surfaces 
are given via theoretical discussions, examples of applications, and extensive 
bibliographies. Ample detail is given so that the reader can understand the 
differences between experimental problems in industry and agriculture (the 
origin of many of the designs and the vocabulary) and so that the special tech- 
niques discussed here can be used by the working statistician. 

I. RICHARD SAVAGE 
University of Minnesota 


A Modern Approach to Intermediate Algebra. By Henry A. Patin. Putnam’s, 
New York, 1958. 224 pp. $3.75. 


According to the preface, this text is a one-semester modern approach to 
intermediate algebra, with two major purposes. The first is “to acquaint the 
student with a logical system as it appears in mathematics.” The preface also 
states that the topics and their order of presentation are fairly traditional and 
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that the presentation, although not rigorous, is coherent. 

Thus, the major new contribution which this book purports to make is the 
development of algebra as a logical system. However, most of this is contained 
on pages 6 through 19, where the postulates for the integers, then the postulates 
for the rationals are stated and the basic rules of ordinary algebra are “proved.” 
The treatment is very terse and abstract and is begun without warning, motiva- 
tion, or any attempt to relate the laws to the student’s past experience. The 
only hint that something is afoot is the assertion, “We will now state our 
fundamental assumptions.” The next thing the reader sees is 


1.1 a+b=b+a 


Unlike certain recent texts on the collegiate level, the problems make no 
contribution toward preparing the student for the next abstract formulation, or 
in helpng him actually devise it; they are mere practice in identifying the laws. 
Here is one: “Show what laws have been used in —3—2=(—3)+(—2) 
= —1(3+2).” While a person acquainted with modern abstract algebra might 
like this presentation in chapter 1, the novice would probably find Chapters 
1 and 2 of Birkhoff and MacLane more accessible. 

Chapter 1 contains the author’s major effort at infusing modern ideas into 
his work. Thereafter, the book grows more traditional, chapter by chapter. By 
the time the one on coordinate geometry was written, modern ideas were 
abandoned completely. For example, the introduction of the concept of func- 
tion is shockingly inadequate for a text published in AD 1958: “An algebraic 
expression in x—3/2x+7/2 is called an algebraic function of x.” There is no 
definition beyond this. 

There are few imaginative or challenging problems. Most could have been 
copied directly from the easier parts of traditional texts. Since one of this 
book’s claims to the modern spirit is that excessive manipulation is avoided, the 
lists of problems contain only those involving easy manipulations. Conse- 
quently the book would not satisfy the teacher of a traditional course, who would 
consider it watered down. 

The final chapter is a four and one-half page discussion of algebra as a deduc- 
tive system. This might be clear enough for a mature reader, but is certainly 
impenetrable to most 11th-grade high school students. 

All prospective authors of modern texts should read this work, since it indi- 
cates how much thought must be given to a pioneering work, by its many 
evidences of having been poorly thought out: little or no motivation, off-the- 
cuff lists of problems, language unsuited to the level of 11th-grade students, 
pretentiously abstruse proofs, introduction of concepts such as variable and 
function without actual definitions, etc. 

Another point strongly made by this book is that the writing of a modern 
text should mot be undertaken unless the author has made contemporary think- 
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ing in all branches of elementary mathematics a part of his working knowledge. 
One final comment. This book could be used in the second semester of a 
traditional 9th-grade course. 
The type, paper and layout are superior; for this the publishers should be 
commended. 
Ep WALTERS 
William Penn High School 
York, Pennsylvania 


Calculus. By Walter Leighton. College Mathematics Series. Allyn and Bacon, 
Boston, 1958. x +373 pages. $6.95. 


The publication of this new volume in the College Mathematics Series fol- 
lows by approximately a year that of the textbook Calculus with Analytic Geom- 
etry by Johnson and Kiokemeister in the same series. The Johnson and Kioke- 
meister text pleased many instructors, but had the disadvantage of being too 
difficult for most students. It contained perhaps an excessive emphasis on rigor. 
The present text goes too far in the other direction, and the author omits proofs 
which he should include. No proofs are given for the standard limit theorems, 
or for the existence of the definite integral. The proof that lim,.» (sin x)/x=1 is 
based upon the assumption that cos x is continuous at x =0, and the author does 
not justify this assumption. The proof given for the Chain Rule is the one which 
most modern instructors try to avoid, and the lack of generality of this proof 
is pointed out by the author in a Critique. The author does not empashize that 
the derivative of a function f(x) is taken with respect to x. Although the defini- 
tion of a derivative is given on page 14, the idea that we differentiate f(x) with 
respect to the variable x is apparently not mentioned until page 58, after the 
discussion and proof of the Chain Rule and the subsequent introduction of the 
differential. In short, the able student will feel that there is material which he 
needs, but which the author has not given him. His discouragement will be 
increased when he consults the references listed at the end of Chapter 1, and 
finds them far beyond him. An excellent reference text for this course is the 
Johnson and Kiokemeister book referred to above. 

In view of the number of proofs which are omitted, it seems unnecessary to 
include an outline of “the basic algebraic theory of the representation of a ra- 
tional function in terms of partial fractions” (see pp. 183-186). 

Much of the discussion is too scanty. The chapter on Solid Analytical 
Geometry contains no adequate discussion of any surfaces except the plane, or 
any curves except the straight line, and contains not a single problem on quadric 
surfaces. The discussion of partial derivatives of functions of functions is rushed, 
and the reader is given the impression that if he wishes to learn the subject 
adequately he should consult other textbooks. Infinite series are motivated by 
Taylor’s formula with the remainder. No discussion is given of operations with 
series. The book contains a short chapter on differential equations. 
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The writing style is good, and the book will be easy for most students to 

read. The problems appear in general to be well chosen. Problem 7, page 38 is 
incorrectly stated, while most students will abandon Problem 5, page 50 before 
they see the point to it. 

Dick WIcK HALL 

Harpur College, 

State University of New York 


A First Course in Statistics. By R. Loveday. Cambridge University Press, Lon- 
don, 1958. x+121 pp. $1.75. 


This is a nonmathematical introduction to descriptive statistics. It does 
not go very far nor does it give a very full discussion of any of the topics it 
contains. It does have worked-out illustrations that might be useful for an 
untrained person operating a desk calculator. It would not be useful in a course 
in statistics having any mathematical prerequisite. In the words of the author, 
“It may ... prove helpful to students . . . who require an elementary introduc- 
tion to the subject before embarking upon a more mathematical treatment.” 

G. Mapow 
Stanford Research Institute 


Some Aspects of Multivariate Analysis. By S. N. Roy. Wiley, New York, 1957. 
viii+214 pp. $8.00. 


In the words of the author, this monograph does not attempt to cover the 
entire area of multivariate analysis but, rather, to deal with some important 
developments with which he has been associated during recent years. In par- 
ticular, he omits a general discussion of such topics as factor analysis, classifica- 
tion problems and the multivariate counterpart of analysis of variance. Some of 
this omitted material will be treated in a future monograph. 

The first two-thirds of the book deals mainly with obtaining confidence 
bounds for certain parametric functions. The testing of hypotheses in so far as 
it is developed is mainly a means to this end. The parametric functions under 
consideration are a set of natural measures of departure from the usual null 
hypothesis, there being a single such function in some cases and a set of such 
functions in more complicated cases. 

The last third of the book consists of appendixes concerning matric theory, 
quadratic forms, linear transformations, jacobians and integration theory. The 
last two pages contain fifty-five references. 

The author is to be commended for the excellence of this monograph. The 
printing is also very good. It is somewhat regrettable however that it is pub- 
lished on 83” X11” paper rather than a more standard size. 

PauL M. HuMMEL 
University of Alabama 
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Scientific Programming in Business and Industry. By Andrew Vazsonyi. Wiley, 
New York, 1958. xix+474 pp. $13.50. 


In this book the author states in the preface “what he (the reader) can ex- 
pect to get is a thorough understanding of what scientific programming is and 
what it means to the business world.” The book is addressed to the business man 
and the writer claims success in presenting the material to business men and 
confidence that this book, seriously studied, is suitable for business men. While 
the reviewer agrees that there are some business men who may read this book 
he feels that it is by far a better book for a college graduate well grounded in 
statistical theory and mathematics than it is for the average business man of his 
acquaintance. 

The book contains chapters on mathematical models, three on linear pro- 
gramming, one each on convex programming, dynamic programming and game 
theory. Several chapters on various kinds of applications are also included. 
Among these are statistical inventory control, transportation allocation, and 
three chapters on scheduling. The treatment is expository and the author goes 
to some length to point out problems not solved by the devices at hand. While 
matrices are avoided generally the author does not explain notations for sums 
nor his nonstandard use of x:,: for x1 for example. Hence the book will not be 
simple for anyone not familiar with certain concepts of linear equations and 
finite dimensional linear spaces. The author is to be commended for the care 
with which he states assumptions. 

The book should be valuable for college and high school teachers of mathe- 
matics who need to have some inkling of the applications of mathematics. It 
should be pointed out that the kind of applications given have been carried out 
mainly by industrial giants and not by small business men. Scientific program- 
ming is still in its infancy as far as both theory and practice go. 

Preston C. HAMMER 
University of Wisconsin 


Applied Mathematics for Engineers and Physicists. (2nd ed.) By Louis A. Pipes. 
McGraw-Hill, New York, 1958. xi+723 pp. $9.50. 


This well-written text is for the student of engineering or physics. It con- 
tains an amazing number of topics in a clear, concise style, thus making it a 
possible introductory source book for mathematical methods applicable to the 
engineering sciences. It does not treat topics depending on probability. It leaves 
the proof of basic mathematical theorems and the elaboration of special cases to 
mathematical analysis and advanced treatments. Much of it is at a slightly 
higher level than that of the student just out of a meager calculus course. There 
is considerable new material not found in the first edition. 

The first chapter covers infinite series up through uniform convergence on a 
slightly higher level than the beginning calculus. Theorems are well stated. Dis- 
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cussion and applications bring out clearly their significance. The use of operators 
on numerous occasions is started here with the exponential operator form for 
Taylor’s series. Each chapter closes with exercises and references. 

A short chapter giving the elementary functions of complex numbers pre- 
pares for the development in the following chapter of Fourier series and the 
Fourier integral in complex form. 

The fourth chapter treats linear algebraic equations, determinants, and 
matrices. The elementary properties of determinants are covered very rapidly. 
A numerical method due to F. Chio for evaluating large determinants is given. 
The work on matrices is far more complete than expected. Silvester’s theorem 
permits the elegant solution of a system of linear differential equations in terms 
of the exponential function of a matrix. The confluent cases, where this method 
becomes especially valuable, is left to the literature. 

Following a chapter on Newton’s and Graeffe’s method for solving equations 
comes one on linear differential equations with constant coefficients. Chapters 
7, 8, and 9 are on “Oscillations of Linear Lumped Electrical Circuits,” “Vibra- 
tions of Elastic Systems with a Finite Number of Degrees of Freedom,” and 
“The Differential Equations of the Theory of Structures.” Here is a profusion 
of applications of the earlier methods and some extensions. For example, one 
finds the analysis of a general network in matrix notation using Laplace trans- 
forms, Lagrange’s equations, stability conditions, and Rayleigh’s method for 
calculating frequencies. 

For systems having repeated components and for numerical work there is a 
treatment of finite differences and difference equations. The work on partial 
differentiation has a discussion of maxima and minima, including Lagrange 
multipliers, and the evaluation of some definite integrals. Material on gamma, 
beta, and error functions prepare for Bessel functions and Legrendre Polyno- 
mials. Legendre functions of the second kind are omitted. 

In “Vector Analysis” there is a clear discussion of curvilinear coordinates. 
Gauss’s theorem and Stokes’s theorem are used to find divergence and curl in 
curvilinear systems. The fundamental equations of hydrodynamics and of heat 
flow in solids are developed. The chapter closes with an introduction to the basic 
tensors and their applications. 

The next three chapters cover the wave equation, Laplace’s equation, and 
the equation of heat conduction. 

The first chapter on complex variables goes through Laurent expansions, use 
of residues, and the argument principle. The other carries mapping through 
Schwartz’s transformation and the applications of conjugate functions to po- 
tential and hydrodynamic problems. 

Chapter 21 gives a more extended treatment of Laplace transforms. The 
p-multiplied form is used so as “to retain the operational forms of Heavyside.” 
A heuristic treatment of the Fourier-Mellin theorem permits finding inverse 
transforms by residues. There are applications to the impulsive functions, par- 
tial differential equations, and integral equations. An appendix gives the basic 
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working theorems and a ten-page table of the transforms “most frequently oc- 
curring in the transient solution of electrical problems.” 

Some methods of treating nonlinear oscillatory problems are given in the 
last chapter. Some of the topics are “Forced Vibrations of Nonlinear Systems,” 
“The Method of Kryloff and Bogolieuboff,” and “Relaxation Oscillations.” This 
61-page chapter closes with 20 references, the most recent having a 1956 date. 

It is a very useful book. 

LAFon 
University of Oklahoma 


Polynomial Expansions of Analytic Functions. By R. P. Boas, Jr. and R. C 
Buck. Springer-Verlag, Berlin, 1958. viii+-77 pp. 19.80 DM. 


This little book differs from most items in the Ergebnisse series in that, 
rather than being a survey of a field, it is a semi-expository presentation of the 
authors’ approach to a topic. 

By a basic set { p,} of polynomials is meant a Hamel basis of the vector space 
of all polynomials. This book deals with classes of analytic functions f which 
have expansions of the form f(z) = }-a,p,(s), where the sequence {an} is defined 
in a prescribed manner. One such procedure, exploited by Whittaker and others, 
consists simply of substitution in the power series for f,so that an = > raf (0)/R!, 
where z*= )-7npa(z), the last summations all being finite since { Pa} is a basic 
set. The authors’ point of view allows other possibilities for the coefficients. The 
main idea is that of kernel expansion, which we describe briefly. 

Given {p,}, let {u,} be a sequence of functions, and define the kernels 
w) = pn(2)un(w). Let T: F-+f be a function-to-function-transformation 
defined by f(z) =1/(2mi){rK(z, w) F(w)dw, T being a suitably chosen path. Then 
f(z) = where a,=1/(27i) Jrun(w)F(w)dw. The simplest example is 
that for which p,=2" itself. Two classical choices for u, are u,=w" and up, 
=w"/n!, so that the respective kernels are (1—zw)— and e*”. With proper 
choice of F, these kernels yield the Taylor series and the Pélya representation. 

The method is applied to a wide class of basic sets which the authors call 
generalized Appell polynomials. Many of the polynomial sets in the literature 
are special cases. Particular attention is given to Appell, Sheffer, and Brenke- 
Huff polynomials. Expansion theorems are obtained involving Borel and 
Mittag-Leffler summability as well as convergence. The form of the typical 
theorem is to characterize the class of functions for which the particular type 
of polynomial expansion converges to the function in the specified sense. 

While the book is by no means light reading, it should be within the scope 
of a large segment of the mathematical public. In particular, it should make an 
excellent seminar topic. 

CASPER GOFFMAN 
Purdue University 
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Introductory Mathematical Analysis. By Edgar D. Eaves and Robert L. Wilson. 
Allyn and Bacon, Boston, 1958. xv +424 pp. $5.95. 


This book is intended to serve the needs of freshmen enrolled in a terminal 
course meeting three hours a week for one academic year. The students are ex- 
pected to have completed at least one year of high school mathematics; how- 
ever, because so much is usually forgotten before the students begin their college 
mathematics, the book opens with a chapter of 37 pages providing a rather 
thorough review of the fundamental operations of arithmetic and algebra. The 
following seven chapters are devoted to an exposition of the concepts of func- 
tion, instantaneous rates, differentiation and integration, with additional topics 
in algebra introduced when they are needed. Then follow chapters on trigo- 
nometry, simple and compound interest, the exponential and logarithmic func- 
tions and their derivatives, sequences, progressions and annuities, and a brief 
chapter on probability. 

The text is well written in a style that is both lucid and thorough. The 
authors amply demonstrate their teaching skill by anticipating students’ diffi- 
culties and by the care with which definitions and theorems are formulated. In 
the reviewer’s estimation the authors can be commended for their good judg- 
ment in the selection of applications likely to interest students at this level and 
for achieving a standard of rigor considerably higher than in most textbooks of 
this category. The exercises, however, are nearly all quite simple, adequate for 
the run-of-the-mill student, but not sufficiently challenging to the abler ones. 

J. M. 
Queens College 


Calculus for Electronics. By A. E. Richmond. McGraw-Hill, New York, 1958. 
viii+407 pp. $6.00. 


This book is designed for use as a mathematics text in technical institutes 
and training-in-industry programs. As the title indicates, it is primarily con- 
cerned with applications. Very few of the fundamental concepts of the calculus 
are described in terms sufficiently precise to be understandable. Many of its 
mathematical assertions are too vague to be called statements; a goodly number 
are false if taken literally. In summary, it is not a mathematics book. 

If it is possible at all for a student to learn the basic techniques of the calculus 
from this book, he can do so only by working his way through the large number 
of illustrative examples and exercises provided by the author. His fear of the 
calculus may be allayed by the fact that most of these examples are selected 
from the subject matter of electronics. Can a student learn the calculus well 
enough to apply it only by learning how to do type-problems? This reviewer 
doubts it. 


M. HENRIKSEN 
Purdue University 
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Studies in the Mathematical Theory of Inventory and Production. By Kenneth J. 
Arrow, Samuel Karlin, Herbert Scarf. Stanford University Press, Stanford, 
1958. x +340 pp. $8.75. 


This book is a collection of papers exceptionally well inter-related. In part I, 
after giving a historical background, the nature and structure of the inventory 
problem are discussed. Part I ends with extensive summaries of the remaining 
fourteen chapters. Part II considers the finding of optimal policies in deter- 
ministic inventory processes. Part III discusses the same for stochastic inven- 
tory processes. Part IV concerns operating characteristics of inventory policies. 

The mathematical tools used range from statistics through the calculus of 
variations, including differential equations and programming. 

As the authors state in Chapter 2, they are interested only in existence 
theorems, etc. as long as they “are useful tools toward arriving at effective solu- 
tions,” that is “effectively computable solutions.” 

This book is a welcome addition to the literature of applications of mathe- 
matics in business. 

ALBERT NEWHOUSE 
University of Houston 


An Introduction to the Theory of Integration. By A. C. Zaanen. Interscience, 
New York, 1958. ix+254 pp. $7.25. 


The author is to be commended for succeeding in the difficult task of writing 
a good book, of readable length, on modern integration. One of its most welcome 
and noteworthy features, which adds greatly to its value and usefulness, is its 
inclusion of certain related topics from classical and functional analysis, in 
which the theory has applications. 

In developing the subject, the author tends to strike a middle position be- 
tween the completely measure-theoretic approach and the more abstract treat- 
ment in which the integral is introduced as a linear functional. He treats the 
basic construction and extension problem in the following fashion. The existence 
of “elementary” integrals of the Daniell type is established by the methods of 
measure theory. The problem of extending an arbitrary elementary integral, 
defined on a linear lattice of real valued functions, is then reduced by the use of 
ordinate sets, to the problem of extending measures on semi-rings to the g-rings 
they generate; and this latter problem is solved by the classical technique of 
outer measures. As a matter of personal taste, this reviewer would have pre- 
ferred a treatment independent of ordinate sets. 

After proving the principal properties of the extended integral and discussing 
the various relations between integrals derived abstractly from linear functionals 
and those arising more concretely from measures, the author proves Fubini’s 
theorem in the form given by Stone. He then gives some applications of the 
theorem by establishing some results on integration by parts, the gamma func- 
tion, and fractional integration. 
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Numerous other topics from analysis and functional analysis, including 
chapters on ergodic theory, unitary transformations on Le, the Radon-Nikodym 
theorem, and Lebesgue’s theorem on the differentiation of an indefinite integral, 
to mention a few, are also in the book. In addition, the book contains an excel- 
lent collection of exercises. - 

R. A. Kunze 
Massachusetts Institute of Technology 


Zahlwort und Ziffer, Bd. Il. By Karl Menninger. Vandenhoeck and Ruprecht, 
Gottingen, 1958. 314 pp. $6.00 (D.M. 24.80) Paperbound. 


This is the second volume comprising the second part of the author’s original 
treatise, A Cultural History of Number. The material of this book can be found 
scattered in other books and journals. The present treatment is a thorough, 
scholarly, and very readable coherent logical-historical development. The his- 
tory of the development of number symbols and computation, both as things 
and as written symbolic manipulations, becomes an exciting drama, which is 
made all the more vidid by hundreds of diagrams, photographs, and drawings. 

The treatise begins with the use of fingers, both to represent numbers up to 
100,000, and to do simple operations on numbers. This is followed by a chapter 
on the way common people represented number by sticks, ropes, and other 
material things, leading to the use of letter symbols and finally unique numerical 
symbols. The way in which abaci and counting boards paved the way to the 
various algorisms performed with written symbols completes the story for the 
Western world. A final chapter treats the development of number language and 
number writing in China and Japan. 

Although written in German, it is simple, yet elegant in style, and can be 
read easily by a second-year student of the language. Volumes I and II, bound 
together, are now available for $9.50 (D.M. 40.00). 

Howarp F. FEHR 
Teachers College, Columbia University 


BRIEF MENTION 


Planning of Experiments. By D. R. Cox. Wiley, New York, 1958. vii+308 pp. $7.50. 


While this book is not, in any sense, a mathematical or statistical book, it certainly 
provides worthwhile reading for the mathematician and statistician, as well as being an 
interesting reference for nonmathematically minded colleagues, who are always with 
us. 


Theories of Figures of Celestial Bodies. By Wenceslas S. Jardetzky. Interscience, New 
York, 1958. xi+186 pp. $6.50. 


A mathematical treatment of the shape of a rotating body. 


Theoretical Mechanics: An Introduction to Mathematical Physics. By J. S. Ames and 
F. D. Murnaghan. Dover, New York, 1958. ix+462 pp. $2.00. 


A republication of the 1929 volume by Ames and Murnaghan. 


l, 

y 

g 

l- 

of 

e 

l- 

ig 

1e 

ts 

in 

e- 

t- 

1e 

ce 

of 

il, 

of 
gs 

of 

e- 

ils 

i's 


744 RECENT PUBLICATIONS [October 


An Introduction to Fourier Methods and the LaPlace Transformation. By Philip Franklin. 
Dover, New York, 1958. x +289 pp. $1.75. 


A corrected edition of the 1949 book Fourier Methods by the same author. 


Introduction to Bessel Functions. By Frank Bowman. Dover, New York, 1958. x+135 pp. 
$1.35. 


A reprinting of Bowman's book published in the late 1930's. 

The Theory of the Potential. By William Duncan MacMillan. Dover, New York, 1958. 
xi+469 pp. $2.25. 

A welcome republication of MacMillan’s 1930 volume. 

The Conduction of Heat in Solids. By H. S. Carslaw and J. C. Jaeger. Oxford University 

Press, New York 16, N. Y., 1959. viii+510 pp. $13.45. 

This revised edition of the well-known standard by Carslaw and Jaeger will be wel- 
comed by all. Call it to the attention of your engineering and physics faculty. 
Dynamics and Nonlinear Mechanics. By E. Leimanis and N. Minorsky. Volume II, 

Surveys in Applied Mathematics. Wiley, New York, 1958. xii+206 pp. $7.75. 


This book is divided into two portions, the first of which deals with recent advances 
in the dynamics of rigid bodies, celestial mechanics, and mathematical exterior ballistics, 
by Leimanis. The second portion by Minorsky deals with the theory of oscillations, in- 
cluding relaxation theory. Both portions conclude with up-to-date bibliographies. 


Dictionary of Astronomy and Astronautics. By A. Spitz and F. Gaynor. Philosophical 
Library, New York, 1959. vi+439 pp. $6.00. 


The understandable definitions in this volume are not limited to a few lines, but often 
run to half or even a full page discussion. 
Vector Analysis with an Introduction to Tensor Analysis. By A. P. Wills. Dover, New 
York, 1958. xxxii+285 pp. $1.75. 
A republication of Wills’ 1931 text. 
The Fundamental Principles of Quantum Mechanics with Elementary Applications. By 
Edwin C. Kemble. Dover, New York, 1958. xviii+611 pp. $2.95. 
A reprinting of Kemble’s 1937 expansion of an earlier work. 


Quantum Electrodynamics. Julian Schwinger, Ed. Dover, New York, 1958. xvii+424 pp- 
pp. $2.45. 


A collection of papers which appeared in various journals during the first half of the 
present century. Well selected, but unfortunately some of the print is too small for com- 
fortable reading. 


Statics and the Dynamics of a Particle. By William Duncan MacMillan. Dover, New 
York, 1958. xviii+430 pp. $2.00. 


A reprint of MacMillan’s 1927 book, which was originally designed as a text book in 
mechanics. 


A Compendium of Mathematics and Physics. By Dorothy S. Meyler and O. G. Sutton. 
Van Nostrand, Princeton, N. J., 1958. x +384 pp. $5.00. 


A collection of formulae from elementary mathematics, physics, and statistics. 
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Linear Groups with an Exposition of the Galois Field Theory. By Leonard Eugene Dick 
son. Dover, New York, 1959. xvi+312 pp. $1.95. 


Since Dickson’s book was written in 1900, Finite or Galois Field Theory has under- 
gone extensive development. Many of the proofs given here could be materially shortened 
and the modern approach would be somewhat different. In spite of the work by Wedder- 
burn, Chevalley, Dieudonné, and Artin, this work of Dickson’s is still eminently read- 
able and very much worthwhile. It could well be required reading for every beginning 
graduate student. 


The Fourier Integral and Certain of Its Applications. By Norbert oe Dover, New 
York, 1959. xi+201 pp. $1.50. 


This reprint of the Cambridge University Press 1933 publication will be of particular 
interest to those who wish to study Cybernetics or Time Series by the same author, since 
much of the related mathematics is given here. 


College Mathematics for Freshmen (2nd ed.). By P. K. Smith and H. F. Schroeder. Van 
Nostrand, Princeton, N. J., 1959. x +314 pp. $4.25. 


This revision of the authors’ 1948 work retains much of the flavor of the original 
book. The title is somewhat misleading, since the content is clearly remedial mathe- 
matics. 


Lectures on the Theory of Elliptic Functions. By Harris Hancock. Dover, New York, 
1958. xxiii+498 pp. $2.55. 


A reprint of Hancock’s 1909 book by the same title. 
Elliptic Integrals. By Harris Hancock. Dover, New York, 1958. 104 pp. $1.25. 
A reprint of Hancock’s 1916 monograph. 


The Foundations of Euclidean Geometry. By Henry George Forder. Dover, New York, 
1958. xii+349 pp. $2.00. 


The current interest in high school geometry has probably spurred the reprinting of 
this 1927 work. 


The Elements of Non-Euclidean Geometry. By D. M. Y. Sommerville. Dover, New York, 
1958. xvi+274 pp. $1.50. 


A reprint of Sommerville’s 1914 lectures. 


NEWS AND NOTICES 


By Lioyp J. Montzinoo, JRr., University of Buffalo 


Readers are invited to contribute to the general interest of this department by sending 
news items to L. J. Montzingo, Jr., Mathematical Association of America, University of 
Buffalo, Buffalo 14, New York. Items must be submitted at least two months before publica- 
tion can take place. 


IBM COMBINATORIAL PROBLEMS INSTITUTE 


An institute on Combinatorial Problems was sponsored by IBM Corporation at its 
Lamb Estate Research Center, near Ossining, New York, from July 6 to August 14, 
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1959. The institute was devoted to studies and seminars in mathematical programming, 
combinatorial design problems, game theory, and switching theory. The following ac- 
cepted IBM’s invitation to attend: Professor A. W. Tucker, Chairman, Princeton Uni- 
versity, Dr. G. B. Dantzig, The RAND Corporation, Professor M. M. Flood, University 
of Michigan, Dr. D. R. Fulkerson, The RAND Corporation, Professor R. E. Gomory, 
Princeton University, Dr. A. J. Hoffman, General Electric Company, Dr. E. F. Moore, 
Bell Telephone Laboratories, Mr. J. Riordan, Bell Telephone Laboratories, Professor 
H. J. Ryser, Ohio State University, Professor C. B. Tompkins, UCLA. 

The results of the Institute will be presented at a symposium, to be held at Princeton 
University in the fall of 1959, jointly sponsored by the Society for Industrial and Applied 
Mathematics and IBM Corporation. 


NSF CONTINUES PROGRAMS 


The National Science Foundation announces that it will continue to support the 
following programs during the summer of 1960 and the academic year 1960-61: 


Undergraduate Research Participation 
Research Participation for Teacher Training 
Science Training for Secondary School Students 


Colleges, universities, and other educational institutions desiring to consider participa- 
tion in these programs may receive full information by addressing: Special Projects in 
Science Education, Scientific Personnel and Education, NSF, Washington 25, D. C. 


PRELIMINARY ACTUARIAL EXAMINATION PRIZE AWARDS 


The winners of the prize awards offered by the Society of Actuaries to the nine under- 
graduates ranking highest on the score of Part 2 of the 1959 Preliminary Actuarial 
Examination are as follows: 

First prize of $200: Harold M. Stark, California Institute of Technology. 

Additional prizes of $100 each: Stephen L. Alder, Harvard University; William G. 
Brown, University of Toronto; Richard M. Dudley, Harvard University; David D. 
Grossman, Harvard University; Bertrand I. Halperin, Harvard University; Ralph E. 
Miller, Harvard University; Steven A. Orszag, Forest Hills (N. Y.) High School; Barry 
Wolk, University of Manitoba. 

The Society of Actuaries has authorized a similar set of nine prizes for the 1960 ex- 
amination on Part 2. 

The 1960 Preliminary Actuarial Examinations will be administered by the Society 
of Actuaries at centers throughout the United States and Canada on May 11, 1960 and 
on November 16, 1960. The closing date for the May examinations is April 1, 1960, and 
for the November examinations, it is October 1, 1960. Further information concerning 
these examinations may be obtained from The Society of Actuaries, 208 South LaSalle 
Street, Chicago 4, Illinois. 


AFFILIATE MEMBERSHIP IN IRE COMPUTER GROUP 


Members of the Mathematical Association of America are cordially invited to 
affiliate with the Professional Group on Electronic Computers of The Institute of Radio 
Engineers. The Professional Group on Electronic Computers is an association of IRE 
members and affiliates with professional interest in the fields of analog and digital elec- 
tronic computers. The group is concerned with the advancement of the electronic com- 
puter field and serves to aid in promoting close cooperation and exchange of technical 
information among its members. 

The fee is $6.50 per year and includes the quarterly Transactions of this professional 
group. Application forms may be obtained from The Institute of Radio Engineers, 1 
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East 79th Street, New York 21, New York. Write attention of Mr. L. G. Cumming, 
Technical Secretary. 
NEW JOURNAL 


The Department of Commerce, through its Office of Technical Services (OTS), 
has started a new journal, TECHNICAL TRANSLATIONS. This journal appears twice a 
month and contains abstracts of foreign scientific publications (including mathematics) 
for which translations (in many cases from the Russian) are available through OTS. 
Further details may be obtained from Office of Technical Services, U. S. Department of 
Commerce, Washington 25, D. C. 


PERSONAL ITEMS 


Fresno State College: Professor Roy Dubisch will be on leave during the academic year 
1959-60, and will be a Lecturer at the University of California, Berkeley; Professor 
G. D. Alkire will be Department Chairman (acting) during the absence of Professor 
Dubisch; Dr. W. L. Allen, Litton Industries, Santa Monica, California, Professor D. L. 
Boyer, New Mexico College of Agriculture and Mechanic Arts, and Dr. W. D. James, 
Stanford University, have been appointed Assistant Professors. 

Seattle University: Dr. B. R. Toskey has been promoted to Assistant Professor; Dr. 
S. A. Husain, Purdue University, has been appointed Assistant Professor. 


Associate Professor B. H. Arnold, Oregon State College, has been promoted to Pro- 
fessor. 

Mr. P. L. Chessin, Westinghouse Electric Corporation, Baltimore, Maryland, has 
accepted a position as Applied Science Representative with the International Business 
Machines Corporation, Washington, D. C. 

Mr. D. M. Estes, Baylor University, has accepted a position as Staff Assistant with 
the Southwestern Bell Telephone Company, Dallas, Texas. 

Dr. T. S. Ferguson, University of California, Los Angeles, is on sabbatical leave and 
will spend the academic year at Princeton University. 

Mr. F. W. Gibson, Douglas Aircraft Corporation, Santa Monica, California, has been 
appointed Graduate Research Mathematician with the Management Sciences Research 
Project, University of California, Los Angeles. 

Dr. Seymour Ginsburg, National Cash Register Company, Hawthorne, California, 
has accepted a position as a member of the technical staff in the Information Processing 
Research Department of the Research Laboratories of the Hughes Aircraft Company. 

Mr. R. B. Hall, Arizona State College, has accepted a position as Research Engineer 
with North American Aviation, Downey, California. 

Associate Professor W. B. Houston, Jr., Morehouse College, has been appointed As- 
sistant Professor at Carleton College. 

Assistant Professor T. R. Jenkins, Washington State College, has accepted a position 
as Research Scientist with the Lockheed Aircraft Corporation, Missiles and Space Divi- 
sion, Sunnyvale, California. 

Mr. J. C. Jones, University of Kentucky, has been appointed Teacher of Mathematics 
at the Cincinnati Public Schools. 

Mr. P. J. Knopp, Harvard University, has been appointed Instructor at Spring Hill 
College. 

Professor Stephen Kulik, Utah State University, has been appointed Professor at 
Long Beach State College. 

Dr. Milton Lees, Institute for Advanced Study, has been appointed Visiting As- 
sistant Mathematician at Brookhaven National Laboratories, Upton, Long Island, New 
York. 

Assistant Professor Kenneth Loewen, Tabor College, has accepted a Graduate School 
Fellowship at The Pennsylvania State University for next year. 
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Mr. L. I. Lowell, Hughes Aircraft Company, Tucson, Arizona, has accepted a posi- 
tion as Research Engineer at the Missile Division of North American Aviation, Inc., 


Downey, California. 


Dr. C. A. McCarthy, Yale University, has been appointed C. L. E. Moore Instructor 


in Mathematics at Massachusetts Institute of Technology. 


Dr. Amin Muwafi, University of Florida, has been appointed Supervisor at the 
Ministry of Education, Amman, Jordan. 


Mr. W. V. Neisius, Systematics, Inc., New York, has been promoted from General 


Manager to Vice-President. 


Mr. C. C. Nielsen, Utah State Agricultural College, has accepted a position as 
Mathematician with the Hercules Powder Company, Magna, Utah. 
Mr. Hugh Noland, Montana State College, has been appointed Assistant Professor at 


Chico State College. 


Mr. G. W. Reitwiesner, Aberdeen Proving Ground, Maryland, has joined the Applied 
Mathematics Division of the National Bureau of Standards, Washington, D. C. 

Dr. J. E. Wilkins, Jr., Nuclear Development Corporation of America, White Plains 
New York, has been promoted to Manager, Research and Development Operations. 


Associate Professor Fred Brafman, University of Oklahoma, died February 4, 1959. 
He was a member of the Association for seven years. 

Associate Professor Myrtle C. Brown, North Texas State College, died March 8, 
1959. She was a member of the Association for thirty-six years. 

Professor Emeritus E. R. Smith, Florida State University, died March 21, 1959. 
He was a member of the Association for thirty-six years. 


THE MATHEMATICAL ASSOCIATION OF AMERICA 


Official Reports and Communications 


Professor H. M. Gehman, Secretary Treasurer, announces that the following 121 per- 
sons have been elected to membership by the Board of Governors on applications duly 


certified. 


Emit_ AMELOTTI, M.S. (Illinois) 
Chairman of Dept., Villanova 
University 

Capt. EpwarD ANLIAN, M.S. (S.U. of 
Iowa) Instr., United States Air 
Force Acade 

Paut ANTHONY, M.A.B.(New York) 
Professor, New York City Com- 
munity College 

Woodland Junior High 
School, East Mead iow, New York 
Ricuarp A. BAKER, Student, Port- 


.(North Texas 
S.C, Mathematician, National 
Bureau of Standards 

C. Beck, M.S. (S.U. of Iowa) 
Instr., Wisconsin State College 

St. Stanislai < Cleve- 
land, Ohi 


Mrs. Lots Browsaav, Student, Wis- 
consin State College 
BLUMSTEIN, M.A. (Buffalo) 


Research Physicist, Cornell Aero- 
nautical Lab. 

Epwarp E. BosMan, M.A. (Columbia) 
Head of Dept., Barrett School, 
Cresco, Pennsylvania 

HELEN M. Bownen, M.A. (Wisconsin) 
Teacher, Bay Ride High School, 
Brooklyn, New York 

ALFRED G. Boyp in E.E. (North 
Dakota) Manager, Marketing 
Planning, General Electric 

M. Pui_srick Brincess, A.M. (Har- 
vard) Head of Dept., Roxbury 
Latin School, Massachusetts 

ALEXIS GUILBAULT, 
M.S.(Montreal) _Instr., Prevost 
High School, Fall River, Massa- 
chusetts; Grad. Student, Co- 
lumbia University 

BROTHER JosepH HeEIsLer, B.S. (St. 
Edward’s) Grad. Student, Uni- 
versity of Notre Dame 

Joun Student, McGill Uni- 
vers 

Brown, Student, Los 


Angeles State Colle; 
F, BROWNELL, B.S. 
Ed. (Northern Illinois) 


Asst., Mankato State College 
Rosert L. CAMPBELL, M.Ed. (Louisi- 
ana S.U.) Instr., University of 


Hawaii 
WALTER J. CARPENTER, M.A. (North 
Carolina) ofessor, 


North Georgia College 

JAMEs W. CEDERBERG, Student, Uni- 
versity of Kansas 

Joun_ C. Cuitwoop, Jr., M.A. 
(Texas) Asst. Professor, Uni- 
versity of Kansas City 

AmBRosE R. CLARKE, Ed.M.(St. 
Lawrence) Head of Dept., 
tavia High School, New York 

Wattace A. M.S. (Wisconsin) 
Asst. Professor, Michigan College 
of Mining and Technology 

Martna J. Coiquitt, A.B. (Mercer) 
Teacher, A. L. Miller Junior 
High School, Macon, Georgia 
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Rosert E, Comey, M.S.(N.Y.S.T.C.,, 
Albany) Instr., Hudson Valley 
Technical Institute, Troy, New 
York 

RoseErt M. Cotton, Student, Swarth- 
more College 

Jose B. Counts, Actuary (Facultad 
de Ciencias Economicas) IBM 
World Trade Corp., Buenos 
Aires, Argentina 

Noe, E. Curr, M.S.(Kentucky) 
Asso. Engineer, The Martin Co. 

Caro. Cummincs, Student, Texas 
Technological College 

Artey T. M.A.(North 
Carolina) Professor and Head 
of Dept., Salem College 

Evan L. Davis, Jr., Student, Uni- 
versity of Colorado 

KENNETH J. Davis, Student, Wofford 


College 

Ernest J. Ecxert, M.A. (Southern 
State Colle: 

Don LS Estes, “Student, Baylor Uni- 


Cuanizs ts FEpeERsSPIEL, M.A. (Mich- 
igan) 1509 St. Mary's Street, 
Raleigh, North Carolina 

Joun A. FiGcueroa, Student, Los 
Angeles State College 

Donap E, FitzGeRaip, A.M. (Yale) 
Instr., Long Beach City College 

G. FLAHERTY, A.B, (Salve 


Regina) Instr., Salve Regina 
Colle: 
A. Geracuty, Ph.D. 


(Notre Dame) Jr. Research 
Instr., University of Notre Dame 

Mrs. ALVA . GotpsmiTH, B.A. 
(Smith) Teacher, Charles E. 
Gorton High School, Yonkers, 
New York 

Mrs. RicHarp H. GRAMANN, B.S. 
(Illinois) Teacher, Arlington 
Township High School, Illinois 

Harvey W. GREENE, M.A. (Missouri) 
Mathematician, Jersey Produc- 
tion Research Co. 

JoHN Ps HA.ey, Student, University 
of Tu 

J. Harrison Hancock, M.S. (Brown) 
Mathematician, Naval Research 


Harry R. Hannu, A.B.A, (Pennsyl- 
vania) Research Manager, The 
Philadelphia Inquirer, Pennsyl- 


vania 

Harvey A. HARTENSTEIN, Student, 
Long Beach State College 

LeonaRD J. Hasster, B.S.(Yale) 
— Pomfret School, Connec- 
ticut 

Jane M. Hit, M.A.(Ohio S.U.) 
Teacher, Washington, D.C. Pub- 
lic Schools 

Nick Hotonyak, Jr., Ph.D. (Illinois) 
Physicist, General Electric, Syra- 
cuse, New York 

Joun H. Huntzincer, M.S. 
Head of Dept., Northeast 
School, Philadelphia, 


vania 
Juri Katviste, B.S. in E.E.(Wash- 
ington) Research Engineer, Boe- 
ing Airplane Co. 
ARTHUR KLEIN, Ed.M. (Boston T.C.) 
Master, Boston Latin School 
Kurtz, Student, South 


C, 
Dakota School of Mines and 
Technology 

Joun P. LASCHENSKI Col- 
lege of the Holy 

MartHa Student, Wisconsin 
State College 


Howarp Levi, Ph.D.(Columbia) 
Asso. Professor, Columbia Uni- 
versity 

F. A. Limouzin, Professor Emeritus, 

le Normales, France 

Mary E, Locan, Student, University 
of Kentucky 

Robert C. Lyngss, B.A. (Oxford) In- 
—— of Schools, London, Eng- 


lan 

LAVAL MATHIEU, Student, University 
of Montreal 

MortHer M. CarRMELITA MCINTYRE, 
Ed.M. (St.Louis) Instr., Spring- 
field Junior College 

LEONARD McKee, B.S. (Cali- 
fornia) Senior Research Chem- 
ist, Shell Chemical Corp. 

Kart Dipl.-Phys. (Han- 
sische) Teacher, Heeresoffizier- 
schule II, Hamburg, Germany 

Joun W. Mitsom, B.A. (Pennsylvania 
S.U.) Grad. Student, South- 
western Louisiana Institute 

Davip B. Moser, Student, Emory 
University 

Preccy J. MuLttns, Student, Univer- 
sity of Alabama 

Francis X. A.B. (Boston 
Coll.) Lt. (jg), Naval Amphibi- 
ous Base, Virginia 

Sanpra G. NEss, Student, Emory Uni- 


versity 

NaTHAN O. Nites, M.S. (St. Law- 
rence) Asso. Professor, Naval 
Academy 

HELEN PatsNner, B.A.(Hunter) Lec- 
turer, Barnard College; Grad. 
Student, Columbia University 

STEPHEN D. Pastor, Student, Oberlin 


ge 

Rust K. N. Pater, M.A. (Columbia) 
Research Asst., Columbia Uni- 
versity 

Georce B. Peck, M.S. (Illinois) Asst. 
Professor, Arizona State Univer- 
sity 

Pencuuck, B.Ed. (Alberta) 
Vice Principal, Ft. Assinboine 
School, Alberta, Canada 

RicHARD W. PIKE, Student, Wiscon- 
sin State College 

Howarp L. Prouse, M.S.(Wayne 
S.T.C.) Supervisor, Math. Dept., 
Mankato State College 

NorMaAN J. PULLMAN, M.A. (Harvard) 
Grad. Asst., Syracuse University 

CuaRLEs W. Ray, B.A. (San Francisco 
S.C.) Grad. Student, San Fran- 
cisco State College 

H. RICHARDSON, Student, 
Chico State College 

H. Roserts, B.E. (Vander- 


bilt) Grad. Student, Emory 
University 

Paut ROETHEL, Student, Wisconsin 
State College 


James H. Ror, Student, Eastern 
Kentucky State College 

Epcar A. Rose, Ed.M. (Rochester) 
Teacher, Monroe High School, 
Rochester, New York 

Mrs. Dorotuy S. RuTLEDGE, B.A, 
Birmingham-Southern) Grad. 
tudent, Emory 

E. Ryuvus, B. alifornia, 
Berkeley) Librarian, 
of California 

Joun Student, University 


sas 

Joun F. H. Scutvep, M.A. (Colum- 
bia) Head of Dept., Cato Merid- 
ian Cen’ School, New York 
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Davi I. Scunerper, Student, Oberlin 


HILBert SCHULTz, Student, Wisconsin 
State College 

NATHAN T. SEELY, Jr., M.A.(Penn- 
— Acting Chairman 

pt., Agricultural and Techni- 

cal College of North Carolina 

Cart J.Sinxe, Ph.D.(Purdue) Asst. 
Professor, Calvin College 

Sister CATHERINE JOSEPHINE, Ph.D. 
(Boston) Asso. Professor, Em- 
manuel College 

SKELTON, Student, Emory Uni- 


rsity 

Ww. SLINKMAN, M.S. 
S.C., Pittsburg) Instr., Univer- 
sity of Colorado 

. Stoss, B.S.(Pomona) Re- 

search Asst., University of Cali- 
fornia, Berkeley 

Cuartes J. Smitu, A.B. (Washington) 
Grad. Student, Washington Uni- 
versity 

GEORGE M. SNEED, Student, Emory 
University 

BERNARD SOHMER, Ph.D. (New York) 
~~ , College of the City of New 


FRANK W. Starks, III, Student, J. M. 
Atherton High School, Louisville, 
Kentucky 

F, STEExE, Asst. Professor, 
Heidelberg College 

Rosert M. Stone, M.Ed. (Howard 

Prof oward 


Payne) Asst. essor, Hi 
ollege 

Epwarp E B.A.(S.C. of 
Washington) Math, Analyst, 


General Electric Co. 
JaMEs J. SWISHER, M. 
rado) Instr., College of the 


Sequoias 
DonaLp G, TEMPLETON, Student, Uni- 
versity of Chattanooga 
Mrs. J uvanita S. Totson, M.S. (How- 
Chairman of Dept., Roose- 
=a High School, Washington, 


F. Tynpatt, B.S, (Franklin 
& Marshall) Grad. Student, 
Brown University 

Harotp R. VANDENBURGH, B.A. 
(North Texas S.C.) Mathema- 
tician, Project Matterhorn, Prince- 
ton University 

Rosert G. VAN Meter, A.M. (Duke) 
Asst. Professor, Geneva College 

Joun R. Veatcu, Student, Whitman 


ollege 
M. WAGNER, M.S. (Iowa 
S.C.) Engineer-Computer, North 


American Aviation 

Norman H. Watton, Student, San 
Jose State College 

Utysses V. Warp, M.A.(Wayne 


— Instr., Howard Univer- 
J. Warne, M.S. (New York) 
Asst., University of Tennessee 


JouN WEISSMAN, M.A.(Temple) In- 
str., Rutgers University 

Marvin J. Winer, Student, Univer- 
sity of Buffalo 

Burton M. Woopwarp, Student, 
University of Florida 

S. Wuncu, Ph.D. 

$80. Professor, 

FREDERIC Wyatt, B.A.(Union) 

. Henry & Co, 

Lituian ZARLING, M.A. (Minnesota) 
Asso. Professor, Wisconsin State 
College 
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THE APRIL MEETING OF THE KENTUCKY SECTION 


The annual meeting of the Kentucky Section of the Mathematical Association of 
America was held at Centre College, Danville, Kentucky, on April 25, 1959. Professor 
W. J. Robinson, Chairman of the Section, presided at both the morning and afternoon 
session. Fifty-three persons attended the meeting including forty-three members of the 
Association. 

At the business meeting the following officers were elected for the coming year: 
Chairman, Professor W. C. Royster, University of Kentucky; Secretary-Treasurer, Pro- 
fessor V. F. Cowling, University of Kentucky; Traveling Lecturer, Professor T. J. 
Pignani, University of Kentucky. 

The following papers were presented: 


1. The numerical solution of the Van der Pol equation, by Professor W. Krogdahl, University 
of Kentucky, introduced by the Secretary. 

Limit cycle solutions of the Van der Pol equation were obtained by numerical integration for 
integral values of the parameter from 1 to 10. The characteristics of these solutions were compared 


with the characteristics of solutions for large parameter values as determined by asymptotic 
formulae. 


2. Singularities of analytic functions, by Mr. J. P. King, University of Kentucky. 
The author uses the transformation z=t+i? to obtain new criteria that the point z=1 bea 
singular point for f(z) = anz*. 


3. Some remarks on grading, by Professor A. G. Anderson, Western Kentucky State College. 

A scheme whereby much of the arbitrariness in grading can be eliminated is the subject of the 
paper. The method is not designed for experienced teachers or advanced courses but should enable 
a new department member to conform to an institution’s grading policy while at the same time 
incorporating ideas of his own. 


4. A class of odd typically real functions in an ellipse, by Professor W. C. Royster, University 
of Kentucky. 

Let denote the class of odd functions (*)f(z) =2+ 2n41(Z), Tn(z) =cos cos™z), 
with the properties that (a) f(z) is regular in the ellipse with foci +1 and semi-axes a>b>0, (b) 
f(z) is real on the real axis and (c) if z lies in E then Im {z}-Im {f(z)} >Oand Re {z}-Re {f(z)} >0. 


It is shown that if f(z) given by (*) belongs to the class T» then deny1R™ +d2n1R-™| <2, n=1, 
2, 


5. A problem in musical composition, by Professor W. H. Spragens, University of Louisville. 
In this paper the author studies problems of intervals as related to the twelve tone scale. 


6. Discrete analytic functions, by Mr. W. T. Sledd, University of Kentucky. 
An expository talk. 


7. Report on the Visiting Lecturer Program, by Professor S. E. Pence, University of Kentucky. 

8. Remarks on MacIntyre’s method of solving a Pfaffian, by Professor T. J. Pignani, University 
of Kentucky. 

An expository talk concerned with MaclIntyre’s (Proc. Edinburgh Math. Soc. vol. 4, 1934- 
1936) method of solving the exact differential equation P(x, y)dx+Q(x, y)dy=0. 

9. Some problems in constructive function theory, by Professor George Piranian, University of 
Michigan. (Invited one-hour address). 

V. F. Cow ine, Secretary 
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THE APRIL MEETING OF THE TEXAS SECTION 


The annual spring meeting of the Texas Section of the Mathematical Association 
of America was held April 17-18, 1959 at the University of Texas, Austin, Texas. There 
were 150 in attendance, including 100 members of the Association. 

The following officers were elected for the coming year: Chairman, Professor W. T. 
Guy, Jr., University of Texas; Vice-Chairman, Professor P. R. Culwell, San Antonio 
College; Secretary-Treasurer, Professor C. R. Sherer, Texas Christian University. The 
next meeting will be held at San Antonio College, San Antonio, Texas. 

Professor H. J. Ettlinger gave a detailed report of the High School Mathematics 
Contest which was given on a trial basis. Through Professor Ettlinger’s expert leader- 
ship, the results were very satisfactory and the section decided to continue the contest 
next year. 

The following program was presented: 


1. Diagonal forms of odd prime degree over finite fields, by Professor J. F. Gray, St. Mary’s 
University of San Antonio, introduced by the Secretary. 

Consider the problem of representing zero by a non-trivial linear combination of pth powers 
in a finite field k. By extending a theorem of D. J. Lewis, it is proved that p pth powers suffice 
when p is an odd prime. If, in addition, p>5, then p—1 pth powers are shown to suffice and, for 
p27, p—2 pth powers are sufficient. In general, where p and fo are odd primes and p> po, there 
exists an integer such that over a finite field k, ++ + = 0(aiEk) 
has a nontrivial solution in k. The values (3)=0 and #(5)=1 are best possible. That po—t(po) 
> logs for an infinity of primes would follow from an undoubted, but unverified, conjec- 
ture of Dickson’s which implies the existence of an infinity of prime pairs (p, g) with g=2p+1. 


2. Some mathematical aspects of the problem of stability, by Professor W. S. McCulley, Agricul- 
tural and Mechanical College of Texas. 

Consideration of an electro-mechanical system described by a second order linear ordinary 
differential equation with constant coefficients indicated that negative real parts for the roots of 
auxiliary equation represent stability. For feedback systems, information about stability may be 
obtained by Routh-Hurwitz criterion, Bode plots, Nyquist criterion, root locus method, describing 
functions, or closed loop pole-zero location. 


3. The accumulation of round-off error in the numerical integration of ordinary differential 
equations, by Professor Peter Henrici, University of California, Los Angeles. (By invitation.) 

The paper studies the dependence of the accumulated round-off error (r.o.e.) in the numerical 
integration of an ordinary differential equation by a general multi-step method in the sense of 
Dahlquist (Math. Scand., vol. 4, 1956, pp. 33-53) on the local r.o.e. Under the assumption that 
the local r.o.e. are independent random variables with known distributions, a formula for the 
variance of the accumulated r.o.e. is derived. The results are tested by numerical experiments in 
which a given differential equation is solved by the same method under varying initial conditions. 


4. An overall look at extended analytic geometry, by Professor R. S. Underwood, Texas Tech- 
nological College. 

When X and Y are particular single-valued functions of variables, a locus on the X Y-plane 
(occasionally vacuous) is implicitly associated with any equation in the variables. Extended 
analytic geometry provides ways of finding these loci rapidly by use of algebra, plane analytic 
geometry, and in some cases, calculus. The system yields solutions, or proves inconsistence, for 
various Diophantine and simultaneous equations, including all linear-versus-quadratic pairs. Trick 
results appear, as well as properties of equations that are independent of plotting rules. Criteria 
for quadric surfaces shorter than the conventional ones are developed. 
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5. Extended quadrics of Darboux, by Professor Dale Maness, Howard Payne College. 

The equation of the most general non-composite quadric hypersurface having contact of 
second order at a general point with an analytic hypersurface in n-space is found. The classical 
tangents of Daroux are found to be specializations of cubic (m—2)-dimensional cones to be known 
as the “cones of Darboux.” The classical directions of Darboux are generalized for hypersurfaces 
and are to be known as the “extended directions of Darboux.” A one parameter family of hyper- 
quadrics to be known as the “extended quadrics of Darboux” is defined and specializes to give 
the Darboux family of quadric surfaces of three space. 


6. The unitary equivalence of 3X3 matrices, by Professor Carl Pearcy, Rice Institute, intro- 
duced by the Secretary. ’ 

Consider the collection A"A™ - ++ A%, where m, m2, m; is any finite sequence of non- 
negative integers, and A is any m Xn complex matrix. Specht showed that the traces of the above 
collection form a complete set of unitary invariants for m Xm complex matrices. Murnaghan showed 
that for n =2 a subset of three of the above traces suffices, and the author has shown that for n=3 
a subset of eleven traces suffices. 


7. Numerical integration with weight functions X and X*, by Professor R. E. Greenwood, Uni- 
versity of Texas. 

The use of the correspondence Pm x*f(x)dx~ (Xin) is advantageous in evaluating first 
moments of means (k=1) and second moments (k=2) for certain empirically determined func- 
tions. The set of needed abscissas {x;..},i=1, - + + , m, may be called a set of generalized theorem- 
of-the-mean points. It is shown that for some values of m one or more of these points may fall out- 
side of the interval (a, b) thereby precluding effective use of the numerical approximation corre- 
spondence. Other values of n yield usable results. This work is a portion of an article by Greenwood, 


Carnahan and Nolley which is scheduled to appear in a forthcoming issue of Mathematical Tables 
and Other Aids to Computation. 


8. Inversion of Jacobian matrices, by Professor Louis Brand, University of Houston. 

In a transformation of coordinates x*=fi(y!,-+-+, y"), i=1,---+, m, the m? derivatives 
dy*/axi may be computed from the fact that the Jacobian matrices (dx*/dy’) and (dy*/dx’) are 
reciprocals when written so that 7 is a row index. The inversion of (dx‘/dy/) is particularly simple 
when the transformation is from one orthogonal system to another. Then (dx*/dy/) can be expressed 
as the product MD of an orthogonal matrix M and a diagonal matrix D, whose elements are the 
lengths of its column vectors; and (dy*/dx/)=(MD)-!=D™ M’ where M’ is the transpose of M. 
The calculation may also be made by observing that the vectors ém:=(dx!/dy', +++, dx"/dy*), 
=(dyi/ax!, +++, dy!/dx*) satisfy e;-e! and hence form reciprocal sets. Thus =gitemr; or 
since gi* =0(j#K), gii=1/g;; for orthogonal transformations, = gliem =€mj/ With 
summation suspended. 


9. Sub-structure vs. super-structure in presenting real algebra on the secondary and elementary 
collegiate level, by Professor J. F. Gray, St. Mary’s University of San Antonio. 

In the current concern to present the real number system, both on the secondary and on the 
elementary collegiate level, from a truly mathematical point of view, two approaches stand out. 
The first proceeds from structure to substructure, beginning with a set of postulates, assumptions, 
or properties of the real number system essentially equivalent to the axioms for an ordered field and 
later descending to study the various sub-systems of interest. The alternative approach proceeds 
from structure to super-structure, ascending from a consideration of the natural numbers to the 


rational and real fields. Reasons and supporting experiences are cited to strongly suggest that the 
first approach is far more desirable. 


10. Equilibrium position of a certain flexible cable, by Professor Guy Johnson, Rice Institute. 
11. Concerning matrix norms, by Professor Ben Fitzpatrick, Jr., University of Texas. 
12. Peculiar sort of power series, by Professor G. R. MacLane, Rice Institute. 
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13. Coverings, by Mr. Terence Reed, Rice Institute, introduced by the Secretary. 


14. Advanced calculus for undergraduate majors, by Dr. Fred Bamforth, Texas Christian Uni- 
versity, introduced by the Secretary. 


C. R. SHERER, Secretary 


THE JUNE MEETING OF THE PACIFIC NORTHWEST SECTION 


The twelfth annual meeting of the Pacific Northwest Section of the Mathematical 
Association of America was held at the University of Oregon, Eugene, Oregon on June 
19, 1959 in conjunction with the 558th meeting of the American Mathematical Society 
and jointly with the Pacific Northwest Section of the Society for Industrial and Applied 
Mathematics. Professor Kenneth Bush, chairman of the Section, presided over the 
meetings. One hundred and eleven persons were in attendance, including 69 members of 
the Association. 

At the Business Meeting the following officers were elected: Chairman, Professor 
A. E. Livingston, University of Washington; Vice-Chairman, Professor T. G. Ostrom, 
University of Montana; Secretary-Treasurer, Professor K. S. Ghent, University of 
Oregon. 

Professor Allendoerfer reported briefly on plans of the Mathematical Association of 
America. 

The program was as follows: 


1. Set theory and numbers, by Professor J. L. Kelley, University of California, Berkeley (by 
invitation). 

The purpose of the paper is to construct the real number system within intuitive set theory. 
The void set is 0, and for each set x we let x’ be the set whose only members are x and the members 
of x. (i.e., yYEx’ iff y¥Ex or y=x). An inductive class A is a class such that OCA and if xGA then 
x'GA. A set n is a natural number iff m belongs to each inductive class. The Peano postulates are 
easily proved. Using the properties of the natural numbers as constructed, addition and multiplica- 
tion are defined as cardinal addition and multiplication. The development is simple and intuitively 
natural. 


2. Boundary value problems in ordinary differential equations, by Professor René De Vogelaere, 
University of California, Berkeley, introduced by the Secretary. 

A general method of numerical analysis, called the discretization method is defined using as 
illustration the problem of integrating a function over an interval. The definition, inspired by work 
of Kantorovich, Kelley, Lemaitre and Crout, contains as special cases a large number of the 
methods used for the solution of differential equations, ordinary and partial, as well as integral 
equations; in particular, difference, Galerkin and least-square methods. The compactness of pres- 
entation that this definition permits is illustrated for the solution of boundary value problems of 
differential equations. Advantage of restricting the functions involved to non-equidistant points 
over the prevailing practice (difference methods) is stressed. 


3. An elementary approach to Lebesgue integration, by Professor H. H. Schaefer, State College 
of Washington. 

This is a brief elementary exposition of a basically well-known method to obtain the concept 
of Lebesgue integral for, say, real functions on a compact interval JC R*. Here “elementary” 
means that no use is made of the notions of measure and measurable function; familiarity with the 
Riemann integral and with the completion of a metric space is assumed. The (didactical) stress is 
laid upon the fact that with this material, along with the basic tools of calculus and using a concept 
equivalent to (Lebesgue) measure 0, the most important properties of the L-integral can be de- 
duced including the convergence theorems of B. Levi, Fatou, and Lebesgue. The method extends 
easily to the Lebesgue-Stieltjes integral and shows rather clearly the relation to the corresponding 
Riemann integrals. 
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4. Some mathematical aspects of the theory of life contingencies, by Professor S. A. Jennings, 
University of British Columbia. 

The theory of life contingencies offers a variety of interesting and important applications of 
elementary calculus which many beginning students might find more meaningful than the con- 
ventional discussions of velocity, acceleration, centroids and moments of inertia. To illustrate this 
contention, the speaker discussed how, starting with ¢p.z, the probability that a person aged x will 
survive to age x+#, the force of mortality may be defined as a logarithmic derivative, and Gom- 
pertz’ and Makeham’s laws of mortality derived by solving a differential equation. Of more interest 
from the point of view of calculus applications, however, are problems involving integration. The 
speaker mentioned several examples of this type. 


5. Mathematics goes to market, by Professor R. E. Gaskel!, Oregon State College. 

Interest of the government in mathematics and mathematicians developed from military 
necessity between 1940-1945, and this has carried over to industry since then. There are, as a 
result, growing numbers of mathematicians employed in industry, in fact the growth could be 
described as explosive. The dramatic effect of high-speed computing machinery is easily visible, 
but fewer people—mathematicians as well as non-mathematicians—realize the no less dramatic 
effect of new developments in mathematics itself on industry. A need for more mathematicians is 
becoming evident, but beyond this more mathematics will be required of all of us—producers, 


consumers, and even those few who insist that they are merely observers. 


K. S. GHENT, Secretary 


CALENDAR OF FUTURE MEETINGS 
Forty-third Annual Meeting, Conrad Hilton Hotel, Chicago, Illinois, January 28-30 


1960 


The following is a list of the Sections of the Association with dates of future meetings 


so far as they have been reported to the Associate Secretary: 


ALLEGHENY Mountain, Grove City College, 
Grove City, Pennsylvania, April 30, 1960. 

ILurnots, Illinois Wesleyan University, Bloom- 
ington, May 13-14, 1960. 

INDIANA 

Iowa, State University of Iowa, Iowa City, 
October 16, 1959. 

Kansas, Kansas State College of Pittsburg, 
April 30, 1960. 

Kentucky, University of Kentucky, Lexington, 
April, 1960. 

LoulIsIANA-MIssISSIPPI 

MARYLAND-DistricT OF COLUMBIA-VIRGINIA, 
American University, Washington, D. C., 
December 5, 1959. 

METROPOLITAN NEw YORK 

MIcHiGAN, University of Michigan, Ann Arbor, 
March 26, 1960. 

MINNESOTA 

Missouri, Central Missouri State College, 
Warrensburg, April 30, 1960. 

NEBRASKA, University of Nebraska, Lincoln, 
April 23, 1960. 

NEw JERSEY, Princeton University, November 
7, 1959. 


NorTHEASTERN, Boston College, Chestnut Hill 
Massachusetts, November 28, 1959. 

NORTHERN CALIFORNIA, University of Cali- 
fornia, Berkeley, January 16, 1960. 

Oxnt1o, Kent State University, May 7, 1960. 

OxvaHoma, Oklahoma City University, October 
23, 1959. 

Paciric NORTHWEST, State University of Mon- 
tana, Missoula, June 17, 1960. 

PHILADELPHIA, University of Delaware, New- 
ark, November 28, 1959. 

Rocky Mountain, United States Air Force 
Academy, Colorado Springs, May 6-7, 1960. 

SOUTHEASTERN, University of South Carolina, 
Columbia, April 1-2, 1960. 

SOUTHERN Ca.irornia, Los Angeles State Col- 
lege, March 12, 1960. 

SOUTHWESTERN, Air Force Missile Develop- 
ment Center, Holloman Air Force Base, 
New Mexico, April, 1960. 

Texas, San Antonio College, April, 1960. 

Upper NEw York State, University of Roch- 
ester, May 7, 1960. 

Wisconsin, Mount Mary College, Milwaukee, 
May 7, 1960. 
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Research 
Scientists and Engineers 
with MS, PhD or ScD 


The Scientific Research Staff Invites 


SCIENTISTS AND ENGINEERS 
TO INITIATE ORIGINAL RESEARCH 


“The Scientific Research Staff of Republic Aviation Is Engaged 
In Performing Theoretical and Experimental Research In the Physical 
Sciences Vital to The Growth of Aeronautics and Astronautics. Qualified 
Individuals are Offered Generous Support In Carrying Out Investigations 
to Demonstrate the Validity of Ideas Leading To Significant 
Advances in The-State-Of-The-Art.” 


Members of Republic’s Scientific Research Staff have been carrying on independent 
investigation in a progressive research environment since the formation of the group 
three years ago. Each individual is encouraged to pursue areas of research in which 
he feels he may make the greatest contributions. 

The ability of this environment to aid in bringing theoretical concepts into the 
realm of feasible engineering has been amply demonstrated. An example is the Plasma 
Pinch Engine conceived by members of this staff. Originally backed by Republic 
funds, it is now receiving supplemental support under government contracts. Among 
other research in advanced stages are programs on lifting fans and new methods of 
structural analysis. 

TODAY THE RESEARCH STAFF 
IS BEING AUGMENTED — 


to emphasize existing areas 

and to open new ones. 

IF YOU FEEL YOU CAN : 

CONTRIBUTE TO THE Supporting Republic’s expanding effort in research and develop- 
WORK OF THIS GROUP, ne R h Center, including 7 modern taboratories 
Republic is ready to discuss will be compteted this yee. @ne hour from New York City, yet 
your interests with you. situated in the center of Long Island, this laboratory is ideally 
Salaries are high, commensurate located for both working and living. 

with talent and creativity. 


Immediate Opportunities Exist for Scientists with Appropriate 
Research Experience and Advanced Degrees in the Following: 


Methods of 
Mathematical Analysis 
Perturbation Procedures Asymptotic Expansions 


Write in confidence directly to Dr. Theodore Theodorsen, 
Director of Scientific Research 


STEP AVIATION 
Farmingdale, Long Island, New York 


der 

1gs, 

of 

‘on- 

this 

will 

om- 

rest 

The 

ary 

isa 

| be 

ble, 

atic 

is is 

ers, 

-30) 

ngs 

Lill, 

vali- 

[on- 

ew- 

orce 

960. 

ina, 

lop- 

ase, 

och- 

& 


Senior Research Scientists 


Appointments are now available on the professional staff of the Naval 
Weapons Laboratory for research scientists with Ph.D. or equivalent ex- 
perience. Openings exist in pure and applied mathematics, theoretical 
physics, and celestial mechanics. Starting salaries range from $7510 
to $11,595 per annum. Naval Weapons Laboratory programs include 
diverse challenging research areas. The most advanced computing equip- 
ment and capable junior scientists are available for assistance. Excellent 
working atmosphere, pleasant community, economical housing and 


recreational facilities. 


For further information, write to the Director, Computation & Exterior 


Ballistics Laboratory. 


NWL 


U. S. Naval Weapons Laboratory 
Department of the Navy 
Dahlgren, Virginia 


Expanding the Frontiers of Space Technology in 


@ The understanding aca solution. of 


complex orbit and trajectory problems 
associated: with space navigation place of 
heavy emphasis on and 
analysis. Lockheed Missiles and Space 
Division is erigaged in several areas of 


In orbit and penta: eaveliite orbit the 
effort is in two ite orbit pre- 
dictions where eed is endeavoring 
to eaiculate near-earth orbits so precisely 
that the gravitational potential of the 


earth and hence its shape and mass dis- anal, 


“tribution may be deduced. This involves 


sophisticated statistical techniques for 


dling tracking information as well 
tical and 


as mathema' 
Second: Lockheed 


investigating trajectories in restricted 
three-body problems in the development 
navigation and vehicle 
Lockheed mathematicians regard is 
bs re and highly intricate science as 
value, rather than an ini 
a problem, which is the present 
approach. 
ENGINEERS and SCIENTISTS 
If you. are experienced in mathemati 
ulus of variations; numeri- 
cal pins Bor ordinary and partial differ- 
ential equations; statistics; or error 
ysis, we invite you to join us in the 
accomplishment of one of the moms Cred 
lenging efforts in man’s history. 
Research and 


J-79, 962 W. 


hLockhheed MISSILES AND SPACE DIVISION 


SUNNY VALE, PALO ALTO, VAN NUYS, 


SANTA CRUZ, SANTA MARIA, CALIF- 


CAPE CANAVERAL. FLORIDA ALAMOGORDO. NEW MEXICO 
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IBM. APPLIED SCIENCE... 


a professional career 
that didn’t exist ten-years ago offers 


BROAD NEW OPPORTUNITIES 


Ten years ago, the Applied Science Representative was almost 
unheard of. Today, at IBM, he has a “frontier’’ career with all 
the associated challenges, excitement, and opportunities to make 
rapid progress based on personal capabilities. 


The IBM Applied Science Representative is a professional 
consultant to the top echelons of business, industry, and govern- 
ment management. He deals with the mathematical formulation 
of scientific and management problems and proves the feasibility 
of solutions with electronic data processing systems. 


Successful men and women candidates will be given comprehensive 
training to equip them to direct their education and talents 
toward a successful career in Applied Science. 


You can become an IBM Applied Scientist, if: 


e you have an M.S. or Ph.D. in engineering, mathematics, or one of the 
physical sciences . . . or a B.S. in engineering, mathematics, or one of the 
physical sciences and a Master’s degree in business administration. 

e you have a B.S. in mathematics or one of the physical sciences, plus 
experience in computer programming. 

e you are interested in exercising a free hand in assignments requiring 
fresh thinking and new, creative approaches. 

e you are able to communicate ideas to technical and nontechnical people. 
e you are interested in joining the leader in the fast-expanding electronic- 
systems field. 


For further information, contact your nearest IBM branch office, or: 
IN THE EAST IN THE MIDWEST 
Mr. W. F. McClelland 
Manager of Applied Science 
IBM Corporation I 
Eastern Region 


Region 
425 Park Avenue 618 South Michigan Avenue _ Wilshire Boulevard 
New York 22, New York Chicago 5, Illinois Los Angeles 5, California 


Applied Science Representatives are located in 
over 100 IBM branch offices throughout the country. 


DATA PROCESSING DIVISION 


fem 
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MATHEMATICIANS WITH ADVANCED DEGREES 


Applied Mathematics Laboratory, 
David Taylor Model Basin has or- 
ganized an Operations Research Di- 
highs to explore the application of 
eed computers to the solution 
of mg ay in war gaming, damage 
assessment, and tracking. 
Other research positions available at 
the Applied Mathematics Laboratory 


grees. A variety of high speed digital 
computers are available on station. 
A LARC computer system is sched- 
uled for delivery in Spring 1960. 
You can become associated with 
some of the outstanding mathema- 
ticians in the field. 

Applications are invited from mathe- 
maticiang with experience and ad- 


» for applicants with MS or PhD.de- vanced * rs 


Write to: Personnel Officer (3 
DAVID LOR MODEL BASIN 
Washington 7, D.C. 
Positions are for duty in Washington, D.C., and 
are in the Career Civil Service. 


Recently published 


—the first in a new series entitled 
UNIVERSITY TEXTS IN THE MATHEMATICAL SCIENCES— 


PROBABILITY AND STATISTICAL 
INFERENCE FOR ENGINEERS: 
A First Course 


by CYRUS DERMAN and MORTON KLEIN, Department of Industrial 
and Management Engineering, Columbia University 

This is a compact introduction to probability theory and statistical infer- 
ence, directed toward engineering students who have had a course in integral 
and differential calculus. The volume has developed from a course taught by 
the authors in the School of Engineering at Columbia University, and pro- 
vides the student with the fundamentals required for professional application 
and further study in areas strongly dependent on probability theory and 
statistical inference. 
July 1959 156 pages illustrated $3.75 


OXFORD UNIVERSITY PRESS, 417 Fifth Avenue, New York 16 
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Announcing 


INTRODUCTION TO 
LAPLACE TRANSFORM 


By Dio L. Holl, Clair G. Maple, and Bernard Vinograde 


This introduction to the Laplace transform provides adequate theoretical 
discussions and applications to problems in the physical sciences and engi- 
neering for undergraduate students. It may be used by students after only a 
year of calculus and a course in differential equations. It establishes the es- 
sential theorems and methods in the first three chapters, and then concen- 
trates on applications in the next two chapters. The sixth chapter contains an 
extension of the theory to finding transforms of functions possessing infinite 
discontinuities. 174 pages, illustrated, price $4.25 


APPLETON-CENTURY-CROFTS, INC. 
35 WEST 32ND STREET, NEW YORK I, N. Y. 


You Haw You Will 
d Texts that 
A reflect Af 
pprove ppreve 
recommendations 
introducti for teaching 
the Foundations & 
Fundamental Concepts e 
by JOHN T. MOORE 
by HOWARD EVES Advance readers specifically commend 
od thematics 
& CARROLL V. NEWSOM this new = 
363 pp. $6.75 tial chapters on SETS and the use of 
our math majors aa toe wen TIONS; the “well conceived” treatment 
and has proved highly satisfactory. We RINEHART of plane analytics de-emphasizing conic 
found students ENJOY this Company sections; the “excellent treatment” of 
i SEQUENCES; COMPLEX NUMBER 
book as well as profiting from it 232 Madison ran i a yea 
I can unreservedly recommend it as New York 16 SYST ‘very ingenious 
excellent for any math major.” interesting PROBLEMS. 
Prof. John G. Kemeny Examination copies will be 
Dartmouth College available in early December 
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A NEW REPRINT 
American Mathematical Monthly 


Now available 

Volumes 21-22, 1914-1915 

Single volumes, paper bound .................. $16.00 
Previously reprinted 

Volumes 1-20, 1894-1913 
$335.00 
300.00 
Single volumes, paper bound ................. 15.00 


Reprinted with the permission of the original publisher 
Please address orders and inquiries to 


JOHNSON REPRINT CORPORATION 
111 Fifth Avenue, New York 3, New York 


Investigate the widest variety 
of openings in recent years 


at the Knolls Atomic Power Laboratory 


For the first time in recent years, current 
-penings here extend into disciplines 
generally considered outside of the tra- 
ditional nuclear areas. As a result, ex- 
cellent opportunities exist today for men 
interested in entering the nuclear field 
for the first time, as well as for recent 
graduates and, of course, experienced 
nuclear engineers and scientists. If 
you’ve been thinking of exploring pro- 
fessional opportunities at KAPL, we sug- 
gest you make your initial inquiry today. 


U.S. Citizenship Required 


Current OPENINGS: 


Advanced engineering mathematics (PhD) 
Advanced | numerieal analysis (PhD) 


. Math I statistics & design of experiments (PhD) 
programming 


Mathematical analysis & computer 
Cc 


Experimental reactor physics 
Theoretical reactor physics 
Nuclear analysis 

Mass spectrometry 


Forward your resume in confidence, 
including salary requirement. 
Please also state your job interests. 
Address: Mr. A. J. Scipione, Dept. 6-MJ. 


GENERAL @® ELECTRIC 
Schenectady, N.Y. 
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Distinguished additions to 
HIB Harper’s Mathematics Series 
Charles A. Hutchinson, Editor 


VECTOR ANALYSIS 


with applications to geometry and physics 


Manuel Schwartz, University of Louisville 
Simon Green, Assumption University of Windsor, Canada 
W. A. Rutledge, University of Tulsa 


A text for beginning students, a reference for advanced students, and a book that 
will be extremely useful to engineers, VECTOR ANALYsiIs is distinguished by a pro- 
fusion of examples and exercises and by a clear, thorough consideration of the 
subject as both a mathematical discipline and a language of physics—the subject 
matter of vector analysis is illumined so that the student can handle problems 
effectively and the book is so designed that the reader can use it on a self-study 
basis. The numerous examples range from simple to complex problems; nearly 
all the problems in the exercises have answers; and there is a wide range of ap- 
plications: statics, kinetics, dynamics, potential theory, electricity and magnetism 
—some elasticity, hydrostatics, and hydrodynamics—techniques in analytic and 
differential geometry. 


556 pp., $7.50 


ADVANCED COMPLEX CALCULUS 


Kenneth S. Miller, New York University 


For a first course in complex variable theory (designed to follow a course in real 
variables such as is presented in the same author's ADVANCED REAL CALCULUS), 
this text is for upper class mathematics majors and beginning graduate students as 
well as for graduate engineering students. “Natural” extensions from the real to 
the complex are carefully pointed out so that complex variable theory appears as a 
)) logical continuatien of real variable theory. The first five chapters are in one-to-one 
correspondence with the first five chapters of the abovementioned text, and intro- 
duce the reader to the basic fundamentals concerning functions of a complex 
variable. The remaining four chapters, departing from this parallel development, 
take up the calculus of residues, important theorems, conformal mapping (in- 
cluding a discussion of the elliptic integral), and a discussion of Laplace integrals. 
235 pp., $5.75 


H|B 


Publishers Since 1817 


49 East 33d Street 
wg | Harper & Brothers New York 16, New York 


TIAA life i insurance . 
the best jeieition for your money 


Teachers Insurance and Annuity Association is a nonprofit insurance company 
founded by Carnegie organizations in 1918. Its primary function is to provide low- 
cost insurance and annuities for employees of colleges, universities, and other 
nonprofit educational and research organizations. 


A $20,000 TIAA 10-Year-Term policy costs only $82.40 per year issued to a man 
aged 34.* This is just one example of the many low-cost TIAA plans available. 
TIAA insurance costs are low because no agents are employed, no commissions 
are paid, and there are few occupational hazards in academic employment. 

A compact reference booklet, the Life Insurance Guide, describes the different 
TIAA policies and is available to help you plan an adequate life insurance pro- 
gram. To get your copy and an illustration of a low-cost TIAA policy issued at 
your age, just return this coupon. 


*$132.20 annual less $49.80 cash dividend at end of . Future 
$49. paid year. dividend amounts cannot 


Employing Institution 


! 730 Third Avenue, New York 17, New York 
| Please send me a Life Insurance Guide and an illustration of low-cost protection at my age. ‘ 
| Name. Date of Birth } 
| Address | 
Ages of pendents 
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PLANE TRIGONOMETRY 


By Gorpon P. Futter, Texas Technological College. Second Edition. 281 pages, $4.75 (With 
Tables) 

This greatly strengthened new edition emphasizes analytic trigonometry and brings it into proper 

balance with numerical trig. It includes clear, simple discussions and explanaions; illustration of 

each new topic with problems worked in detail; practical problems taken from physics, surveying, 

and aviation; a complete and adequate treatment of logarithms; a discussion of trigonometric equa- 

tions and inverse functions; a full chapter devoted to complex numbers. 


BOOKS BY ROSS R. MIDDLEMISS, Washington University 
COLLEGE ALGEBRA 
344 pages, $4.50 


This excellent text contains a complete coverage of topics usually taught in a standard course. The 
course is made more valuable and stimulating by the greater emphasis on reasoning and clear 


thinking; this method combats the student’s tendency toward mechanical operations unaccompanied 
by real thought. 


ALGEBRA FOR COLLEGE STUDENTS 
394 pages, $4.75 


A new reatment of the author’s College Algebra, this text is designed for the less advanced students. 
For a slower, more detailed study, the fundamental material—through quadratic equations with 
one unknown—has been expanded. The lessons have been shortened and geared in treatment to a 
somewhat less mature student with a background of only one year of high school algebra. Emphasis 


is upon a real understanding. 
ANALYTIC GEOMETRY. Second Edition. 
310 pages, $4.50 


Written especially for students preparing for calculus, this text emphasizes graphs and the elemen- 
tary properties of polynomials, rational fractional functions, exponential and logarithmic functions, 
and trigonometric functions. In this second edition, new problems, new figures, new illustrative 
examples, and sometimes a new point of view are introduced. 


DIFFERENTIAL AND INTEGRAL CALCULUS. Second Edition. 
497 pages, $5.95 


This important text presents the topics customarily covered in the first course in the calculus. 
Special attention has been given to topics which are inherently difficult, and in many places the 
explanation has been expanded for greater clarity. A chapter on solid geometry has been added. 
Many chapters have been entirely rewritten. 


Send for copies on approval 


McGraw-Hill Book Company, Inc. 
330 West 42nd Street : New York 36, N. Y. 
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texts in the new mathematics series 


SESE 
QAAAAN 


under the direction of CARL B. ALLENDOERFER 
INTRODUCTION TO MATHEMATICAL STATISTICS 


by ROBERT V. HOGG and ALLEN T. CRAIG, 
both, University of Iowa 


ready November 1959 


Designed for advanced students in mathematics, this book presents the 
fundamental concepts of statistics and probability in a reasonable and nat- 
ural order, showing interrelationships among them wherever possible. This 
work contains modern theoretical developments not found in books written 
at this mathematical level. An answer pamphlet will be available gratis. 


MODERN MATHEMATICS: AN INTRODUCTION 
by SAMUEL I. ALTWERGER, The New School for Social Research 
ready February 1960 


Designed as a basic integrated text in mathematics for liberal arts students, 
the contents are logically arranged on an axiomatic basis, permitting the 
ooo gene of topics ranging from elementary arithmetic (number theory) 
through portions of the calculus. This mature presentation contains man 

contributions from modern mathematics. A solutions manual will be avail- 


able gratis. 
COLLEGE ALGEBRA 
by WILLIAM H. DURFEE, Mount Holyoke College 
ready Spring 1960 
theory and solution of 
ORDINARY DIFFERENTIAL EQUATIONS 


by DONALD GREENSPAN, Purdue University 
ready Spring 1960 


The Macmillan 


60 FIFTH AVENUE, NEW YORK 11, N.Y. 


GEORGE BANTA COMPANY, INC., MENASHA, WISCONSIN, U.S.A. 
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